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1  INTRODUCTION 

In  attempting  to  describo  the  motion  of  an  earth  satellite  the  basio 
problem  is  to  express  analytically,  and  with  sufficient  accuracy,  tho  deviation 
of  the  motion  from  the  simple  path  that  would  be  described  if  tho  only  force 
acting  were  due  to  a  uniform  gravitating  sphere.  Three  of  the  actual  perturbing 
forces  arise  from  the  asphericity  of  the  earth,  from  attraction  by  other  bodies  - 
notably  tho  sun  and  moon  -  and  from  the  pressure  of  solar  radiation.  The 
effects  of  these  three  sources  of  perturbation  aro  studied  in  this  Report, 

Other  sources  of  perturbation,  not  studied  here,  include  atmospheric  drag  and 
small  forces  due  to  electromagnetic  drag  and  Einstein's  modification  of  Newton's 
law  of  gravitation. 

Since  the  first  satellites  wore  launched  in  1957  a  spate  of  papers  on 
satellite  perturbations  has  appeared.  The  effects  of  the  earth's  asphericity, 
of  luni-solar  attraction  and  of  solar  radiation  have  customarily  been  treated 
separately  and  there  is  a  particularly  large  literature  on  the  first  of  these. 

The  object  of  the  present  paper  is  the  development  of  a  general  analysis  from 
which  results  applicable  to  all  three  sources  of  perturbations  may  be  obtained. 

Most  writers  of  papers  on  satellite  perturbations  assume  their  readers 
have  a  wide  background  knowledge  -  in  particular  an  understanding  of  the 
Legendre  method  of  expanding  the  earth's  potential.  In  the  hope  of  making  this 
Report  more  comprehensible  to  those  who  lack  this  background,  and  to  make  the 
account  as  complete  as  possible,  the  opportunity  has  been  taken  to  give  some 
of  the  fundamental  material,  mainly  in  tho  Appendices  and  sometimes  rather 
unconventionally, 

A 

Merson  lias  investigated  the  motion  of  a  satellite  of  a  parent  body  whose 
gravitational  field  is  axi-syrrmctric,  by  expanding  the  potential  at  a  general 
exterior  point  as  a  series  of  terms  involving  Legendre  polynomials  and  negative 
powers  of  the  radial  distance  to  tho  point  from  an  origin  on  the  axis  of  - 
symmetry.  Taking  the  earth  to  be  a  parent  body  of  this  type  he  has  obtained 
the  main  effects  of  the  first  six  parameters  of  asphericity  -  tho  J  coefficients  - 
on  the  orbital  elements  of  an  earth  satellite, 

2 

Cook  has  considered  tho  perturbations  of  an  earth  satellite  duo  to  tho 
attraction  of  a  distant  body  (the  sun  or  moon),  Ko  assumes  the  distant  body 
does  not  move  during  one  complete  revolution  of  the  satellite  and,  like  Merson, 
describes  tho  perturbations  by  giving  the  changes  in  the  orbital  elements  of 
tho  satellite  over  such  a  revolution. 


Now  the  disturbing  field  arising  from  a  stationary  body,  considered  as  a 
point  mass,  is  axi-aymme trio.  Hence  the  effect  on  satellite  motion  may  be 
obtained  by  generalizing  Merson's  theory,  using  a  series  of  terms,  in  the 
Legendre  expansion  of  potential,  with  positive  powers  of  the  radial  distance, 
Cook's  results  being  for  the  .  irtt  term  in  the  series.  Such  a  generalization 
is  made  in  this  Report,  the  potential  being  expressed  as  a  doubly  infinite 
series,  powers  of  the  radial  distance  running  from  -co  to  ■*<».  The  effect  of 
the  general  term  of  the  series  is  studied.  When  results  are  applied  to 
satellites  of  the  earth,  asphericity,  lunar  and  solar  perturbations  must  be 
dealt  with  separately  since  there  is  a  different  axi3  of  symmetry  associated 
with  each  (the  earth's  polar  axis  and  axes  pointing  towards  the  moon  and  sun 
respectively).  It  is  a  straightforward  matter,  however,  to  transform  results 
associated  with  one  axis  into  the  appropriate  form  relative  to  another  axis 
and  it  is  natural  to  regard  the  earth's  polar  axis  as  the  fundamental  reference. 

Apart  from  a  consideration  of  secular  and  long-period  perturbations  in 
section  10,  the  Report  follows  the  J.ferson-Cook  approach  of  studying  perturbations 
over  a  complete  revolution  of  a  satellite.  The  treatment  owes  much  to  the  work 

3  4 

of  Groves  and  Kaula  .  These  authors  give  very  general  accounts  of  satellite 
motion  in  the  earth's  field,  with  axial  symmetry  no  longer  assumed  and  the 
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short-period  effect  of  every  Legendre  term  included.  In  a  subsequent  paper 
Kaula  gives  an  equally  general  development  of  the  lunar  and  solar  disturbing 
functions.  The  very  generality  of  these  papers,  however,  involves  two 
disadvantages:  the  analysis  is  formidable  from  the  point  of  view  of  the 
average  reader,  and  there  is  a  dearth  of  concrete  results  at  the  end;  algebraic 
expressions  are  given  which  cover  all  first-order  terms  (of  the  perturbations 
of  the  elements),  but  it  is  very  difficult  to  write  down  specific  terms  frcm 
these  expressions. 

In  this  Report,  therefore,  the  aim  has  been  to  make  the  general  formulae 
as  simple  as  possible  and  such  that  the  generation  of  particular  expressions 
is  always  an  easy  and  obvious  process.  To  the  former  end,  and  despite  the 
extension  to  positive  powers  of  the  radial  distance  already  mentioned,  the 
formulae  have  been  kept  less  general  than  those  of  Groves  and  Kaula  in  that 
short-period  terms  are  averaged  out  as  soon  as  possible,  while  the  restriction 
to  axially  symmetric  terms  is  itself  an  enormous  simplification.  The  important 
short-period  terms  corresponding  to  Jg  only  -  i,  e,  to  the  dominant  oblateness 
term  in  the  case  of  the  earth  -  are  listed  separately  in  section  10.1  of  the 
paper.  To  assist  in  the  derivation  of  particular  results,  each  of  the  five 
basic  functions  defined  in  the  Report  is  tabulated  as  far  as  the  fourth 
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harmonic  in  the  Legendre  expansion.  The  functional  expressions  for  higher 

harmonics  are  best  obtained  by  the  use  of  recurrence  formulae;  a  list  of  the 

6 

required  formulae  is  given,  generalizing  some  work  of  Merson  • 

Prom  the  formulae  for  the  perturbations  of  on  earth  satellite  after 
one  complete  revolution,  expressions  are  derived  for  secular  and  long-period 
variations  valid  over  many  revolutions.  Formulae  for  luni-solar  perturbations, 
relevant  to  a  not-too-distant  satellite,  are  given  for  three  situations:  for 
a  complete  revolution  of  the  satellite,  the  sun  or  moon  being  considered 
stationary;  for  the  average  over  a  revolution  of  the  sun  or  moon;  and  for 
further  averages  with  respect  to  revolutions  of  the  perigee  of  the  satellite 
and  of  it3  line  of  nodes  relative  to  the  line  of  nodes  of  the  sun  or  moon. 

2  PREVIOUS  WORK 

The  considerable  literature  dealing  with  satellite  perturbations  has 

been  remarked  upon  and  some  of  tho  papers  (Refs.  1-6)  mentioned.  It  may  be  of 

interest  to  give  a  brief  survey  of  some  of  the  more  important  contributions 

to  this  literature  and  we  start  with  those  which  tackle  the  asphericity 

perturbations.  Fairly  detailed  descriptions  of,  and  comparisons  between,  the 

7  8 

main  theories  are  provided  in  useful  papers  by  A.  H.  Cook'  and  Kaula  . 

9 

First,  historically,  was  a  paper  by  King-Hele  .  His  results  are  limited 

mainly  by  the  neglect  of  higher  pavers  of  the  eccentricity  and  by  the 

assumption  that  the  existence  could  be  postulated  of  a  constantly  inclined 

"orbital  plane"  in  which  the  satellite  would  always  lie.  The  theory  of 

10 

Brenner  and  Latta  provides  a  natural  extension  of  King-Hele's  methods,  the 
characteristic  of  this  method  being  its  reliance  on  first  principles  rather 
than  on  the  traditions  of  celestial  mechanics, 

^arson's  theory  has  provided  the  basis  for  all  orbital  analysis  so  far 

conducted  at  the  R.A.E.  He  carries  to  its  logical  conclusion  the  idea  of 

introducing  non-osculating  elements  which,  as  far  as  Jg  is  concerned,  are  as 

smooth  as  possible  over  one  revolution  of  the  satellite.  The  starting  point 

is  the  set  of  six  Lagrangian  planetary  equations  for  osculating  elements,  as 

it  is  also  with  the  theories  of  several  other  authors  including 

12 

Zhongolovitch  (see  Ref.  1l)  and  Kozai  ,  Kozai’s  theory  is  the  one  used  by 
the  Smithsonian  A3trophysical  Observatory  to  analyse  Baker-Nunn  camera 
observations. 

Kozai’s  paper  is  to  be  found  in  a  noteworthy  number  of  the  Astronomical 
Journal  which  also  gives  the  theories  of  Brouwer”*  ^  and  Garfinkel^.  These 
authors  have  more  advanced  starting  points  which  involve  the  methods  of 

r-  i 
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Hamilton,  Jaoobi,  Delaunay  and  Von-Zeipel.  An  intermediate  orbit  ia  required  - 
the  baaio  Kepler  ellipse  in  the  case  of  Brouwer,  but  a  pseudo-ellipse,  much 
closer  to  the  true  orbit,  in  the  oase  of  Garfinkel,  Brouwer's  theory  is  used 
by  the  National  Aeronautics  and  Space  Administration  (NASA)  as  one  basis  for 
the  derivation  of  orbits, 

a  r  4/ 

The  classical  method  of  Hansen  has  boen  modified  by  Musen  to  give  a 
theory  which  is  part  analytical,  part  numerical  and  is  intended  to  facilitate 
the  practical  representation  of  satellite  motion  to  very  high  accuracy, 

17 

Musen ' s  theory,  also,  has  been  used  by  NASA  and  formulae  developed  by  Fisher 
allow  it  to  be  compared  directly  with  Brouwer's, 


No  survey  would  be  complete  without  mention  of  the  series  of  papers  by 

jO 

Vinti  of  the  U.S.  Bureau  of  Standards,  In  the  first10  of  the  three  main  papers 

of  this  series  Vinti  introduced  spheroidal  co-ordinates  in  place  of  the  usual 

spherical  co-ordinates.  By  this  means  he  obtained  separability  of  the 

Hamil  ton-Jacobi  equation  for  motion  in  a  field  which  accounts  for  more  than 

99  per  cent  of  the  earth's  asphericity.  In  terms  of  the  J  coefficients 

(see  section  4.1)  the  field  is  such  that  J 2n  =  (-l)n+1  Jg  arxi  J2n+1  =  and 

hence  can  represent  the  dominant  J 2  coefficient  completely,  residual  deviations 

being  due  to  the  values  of  +  Jg,  Jg  -  J^f  etc,  Vinti  originally  had  a 

second  free  parameter  in  the  potential  so  that  J y  say,  could  be  represented  by 

a  choice  of  origin  with  *  0  and  ^n+1  =  ^2»  but  this  idea  was  no* 

pursued  beyond  the  first  paper.  In  the  second  main  paper1'*  the  formal  integrals 

from  the  first  are  evaluated,  two  possible  sets  of  orbital  elements  defined 

and  the  exact  solution  for  the  intermediate  orbit  developed.  In  the  third 
20 

paper  Von  Zeipel's  method  is  used  to  obtain  the  perturbations,  relative  to 
tho  intermediate  orbit,  due  to  J.  +  ard  J,. 

W  J 


On  passing  to  the  subject  of  luni-solar  perturbations  one  finds  far  fewer 
papers.  The  main  reasons  for  this  are  that  for  a  close  earth  satellite  the 
perturbations  are  much  smaller  than  the  asphericity  perturbations  and  that  the 
subject  has  less  obvious  novelty  and  so,  perhaps,  less  intellectual  attraction. 
Perturbations  of  the  lunar  orbit  about  the  earth  due  to  the  sun  and  other 
disturbing  bodies  have,  of  course,  been  studied  for  many  years  and  the  theories 
of  Hansen,  Delaunay,  Hill  and  Brown  are  described  in  such  standard  texts  as 
Ref s,21,  22  ard  23.  Although  the  problem  of  luni-solar  perturbations  of  an 
earth  satellite  is  superficially  the  same,  it  is  actually  different  enough  for 
most  workers  to  have  attacked  it  independently  of  the  traditional  theories. 
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The  paper  of  Cook  ,  alroady  mentioned,  considers  only  the  main  term  in 

tho  Legendre  expansion  of  tho  disturbing  function  and  takos  the  disturbing 

body  to  bo  fixed.  Allan^,  using  vector  methods,  considers  a  further  term 

25  25 

in  the  expansion  -  tho  parallactic  term  -  and  Smith  yet  anothor.  Smith 
also  allows  for  the  motion  of  the  disturbing  body  and  for  long-period  inter¬ 
action  with  tho  secular  effects  of  the  earth’s  asphoricity. 

26 

Kozai  develops  tho  disturbing  function  for  tho  effect  of  a  disturbing 
body  at  a  general  position  in  its  orbit,  working  with  the  main  term  only. 

Muson,  Bailie  and  Upton  '  include  the  parallactic  torm  in  a  similar  development. 
Rci,27  includes  sliort-poriod  terms  und  requires  13  pages  for  the  expression 
of  the  disturbing  function  in  toms  of  orbital  parameters.  With  secular  and 
long-period  terms  only  it  still  requires  6?  P ages;  of  these,  2?  pages  cover 
the  main  torm  in  tho  Legendre  expansion  (Types  1  to  5  on  pp.  25-27  of  Ref,  27) 
and  agree  with  Kozai’ 3  development  (pp,  8-9  of  Ref, 26)  apart  from  a  sign  error 
in  one  of  Kozai’ s  terms  (line  13  on  page  9  of  Ref. 26). 

(After  the  rest  of  the  present  paper  was  written,  the  author  confirmed 

that  the  error  in  the  disturbing  function  just  referred  to  i3  indeed  in  Ref, 26 

and  not  in  Ref.  27.  This  was  done  by  comparison  with  a  third  paper,  recently 

42 

written  by  Myrna  Lewis  ,  Her  expression  for  tho  disturbing  function  - 
equation  (l8)  of  Ref. 42  -  itself  contains  a  copying  error,  which  is  not 
transmitted  to  the  rest  of  the  paper,  in  that  the  middle  term  of  its  third  line 
should  be  multiplied  by  2.  Ref.  42  goes  beyond  Ref 3, 26  and  27  by  giving  a  full 
development  of  secular  and  long-period  perturbations  to  the  orbital  elemonts. 
Ref.  27  gives  no  such  development  and  Ref.  26  gives  only  the  secular  toms  for 
Q  and  w). 

3  LAGRANGE’S  PLANETARY  SQUAT! QMS 

Wo  suppose  tliat  a  satellite  is  moving  in  a  gravitational  field,  the  main 
effect  of  which  is  an  inverse-square-law  attraction  towards  a  fixed  point,  0, 
Then  if  the  field  is  conservative,  the  potential  may  bo  written  |i/r  +  U  where  r 
is  the  distance  from  0  and  ^  is  a  constant.  The  first  term  is  from  the  central 
force  and  the  second  torm,  U,  from  the  rest  of  tho  field.  This  term  is  called 
the  disturbing  function  and  we  shall  assume  it  to  bo  a  function  of  position 
but  not  of  time.  In  the  absence  of  perturbing  force,  the  satellite  would 
describe  a  fixed  Kepler  conic  section  orbit  -  wo  suppose  an  ellipse.  The 
cllipso  into  which  the  actual  orbit  would  "freeze’’  if,  at  time  t,  tho  perturbing 
force  wore  instantaneously  removed  is  called  the  osculating  ellipse  at  time  t. 

It  can  be  specified  by  six  elements;  if  tho  gravitational  field  is  that  of  the 
earth  (centre  0)  these  are  normally  taken  as  the  following  set:  send -major 


axis  a,  ecoentricity  e,  inclination  i,  right  ascension  o’  the  node  fl,  argument 
of  perigee  u  and  moan  anomaly  at  the  epoch  (when  t  =  0)  x* 

It  is  important  to  note,  to  avoid  confusion,  that  \  -  as  with  the  other 
five  elements  -  is  defined  for  the  current  time  t  and  varies  -  like  the  others 
duo  to  the  perturbing  force;  the  definition,  however,  involves  an  extension  of 
the  instantaneous  "frozen"  orbit  back  to  t  =  0,  Use  of  this  element  has  an 
inherent  drawback.  For  if  Ax  is  the  change  in  %  between  two  times  t^  and  t2, 
Ax  will  depond  on  two  extrapolations  back  to  t  a  0.  But  for  an  arbitrary 
displacement  of  the  time  origin  (say  to  a  system  defined  by  t'  =  t  +  const), 
the  extrapolations  would  be  back  to  a  different  epoch  (t'  =  0)  and  so  A%  would 
have  a  different  value.  This  is  clearly  mo3t  undesirable. 

To  avoid  this  difficulty  we  replace  the  element  X  by  another  element, 

usually  denoted  by  x*  but  here  by  <r  for  ease  of  notation,  such  that,  if  n  is 

2  3 

the  mean  motion  given  by  n  a  =  p, 

t 

(T+Jndt  =X  +  nt  »  0 ) 

o 

where  both  sides  of  this  equation  give  the  mean  anomaly,  M,  in  the  osculating 

orbit.  The  element  <r  itself  suffers  from  the  disadvantage  that  it  is  not,  as 

is  x>  immediately  defined,  for  a  given  t,  from  the  osculating  orbit.  By  its 

23 

use,  however,  it  can  be  shown  that  the  equations  for  the  variations  of  the 
elements  arc  easier  to  integrate. 

23 

These  equations,  Lagrange's  planetary  equations  ,  may  now  be  stated  m 
the  following  form:- 
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In  these  equations  the  disturbing  function  U  is  assumed  to  be  a  function  of 
a,  e,  i,  0,  GJ  and  M,  where 

t 

M  =  o'  +  J  n  dt  ,  (8) 

o 

so  that  the  dependence  of  n  on  a  does  not  enter  into  ~  # 

oa 

Once  the  equations  have  been  set  up  there  is  no  difficulty  in  picking  out 
the  main  first-order  terms  of  the  solution*  For  we  simply  set  a,  e,  i,  Q  and  u> 
constant  on  the  right  hand  sides,  after  the  partial  derivatives  have  been 
formed.  To  eliminate  short-period  terms  we  then  integrate  over  a  complete 
revolution  of  M,  i.e,  over  one  orbit,  (Short-period  terms  are  included  in  the 
very  general  treatments  of  Groves^  and  Kaula^j  for  a  satellite  in  the 
gravitational  field  of  the  earth,  however,  the  only  significant  short-period 
terms  arise  from  -  defined  in  the  next  section  -  and  these  are  given 
explicitly  in  section  10,1), 

4  THE  DISTURBING  FUNCTION 

4* 1  Form  assumed 


The  earth's  gravitational  potential  may  be  expressed  as 


12 


£*U 


oo  n 

'['  *  X!  E®"  ^  (8ln  @)  <Cn,m  003  3X  *  V  Sln  n,X)] 


n=1  TftsO 


This  only  differs  from  the  form  and  notation  recornnendod  by  the  International 

pQ 

Astronomical  Union  in  that  U,  in  equation  (9),  is  used  just  for  the  disturbing 

function  and  not  for  the  entire  potential.  The  derivation  of  the  expression 

is  outlined  in  Appendix  A.  The  notation  Is  as  follows:  C  and  S  are 
rr  n,m  n,m 

numerical  coefficients,  r  is  the  distance  from  the  earth's  centre  of  mass, 

R  is  the  mean  oquatorial  radius  of  the  earth,  (3  i3  latitude,  X  is  longitude 
measured  east  and  (sin  3)  is  the  associated  Legendre  function  defined  in 
Appendix  B.  The  expansion  is  valid  for  points  outside  a  sphere  containing  the 
total  mass  of  the  earth;  neglecting  the  mass  of  the  atmosphere,  the  radius  of 
the  smallest  such  sphere  is  only  slightly  in  excess  of  R. 

An  important  special  case  arises  if  the  assumption  is  made  that  the 

gravitational  fie-ld  is  symmetrical  about  the  earth' 3  polar  axis.  With  the 

28 

same  deviation  from  the  recommendations  of  the  I.A.U.  as  before,  the 
potential  may  then  bo  written 


OO 


£)  Pn  « 


(10) 


where  is  a  numerical  coefficient  and  P^  (sin  3)  =  P°  (sin  3)  is  Legendre’s 

polynomial.  Equation  (9)  clearly  reduces  to  equation  (lO)  on  taking 

C  =  -J  and  all  other  C  and  S  coefficients  zero, 
n,o  n 

The  assumption  of  polar  symmetry  gives  a  very  good  approximation  to  the 
actual  field  of  the  earth.  It  leads  to  considerable  simplification  in 
integrating  Lagrange's  equations.  One  reason  for  this  is  that  the  three 
'orientating'  elements  in  the  set  of  orbital  elements  defined  in  section  3  - 
that  is  to  soy,  i,  fi  and  u  -  are  defined  relative  to  the  assumed  axis  of 
symmetry  (the  polar  axis)  and  the  plane  perpendicular  to  it  (the  equator). 

In  this  paper,  though  we  are  considering  perturbations  associated  with 
arbitrary  central  fields,  the  standard  application  is  to  the  field  of  the 
earth  and  it  i3  convenient  to  use  the  notation  which  refers  to  the  earth. 

The  restriction  to  axial  symmetry  is  assumed  so  that  the  potential,  outside 
a  sphere  containing  the  central  mass,  is  described  by  equation  (lO),  However, 
though  the  generalisation  to  equation  (9)  is  not  made,  we  wish  to  generalize 
equation  (lO)  another  way.  Gravitating  mass  is  assumed  to  exist  outside  a 


13 


sphorc  of  radius  ^  as  well  as  inside  a  sphere  of  radius  r^, 
two  spheres  is  at  0,  a  point  on  the  axis  of  symnetry  but  not 
centre  of  mass  of  either  region. 

Then  the  potential  may  be  taken  to  be  given  by 

£  +  u  -  &[i  -  pn(?J  pe  <sinp>] 

n=-oo 

where  C=. n  if  n  »  0,  and  t  =  -n-1  if  n  <  0;  R  is  now  an  arbitrary  fixed  value  of 

r.  This  expansion  is  valid  at  all  points  betv/een  the  two  spheres,  i.e,  for 

which  r .  <  r  <  r0.  The  terms  with  n  >  0  arise  from  the  interior  mass  and  the 
1  2 

terms  with  n  <  0  from  the  exterior  mass,  A  derivation  of  the  expansion, 
extended  to  the  case  of  axial  non-symmetry,  is  given  in  Appendix  A. 

Throe  of  the  J  coefficients  may  be  eliminated  very  simply.  Since  Pq  (sinf3) 
is  constant  (+l)  is  associated  with  a  constant  contribution  to  U  which  may 
bo  disregarded.  Jq  will  be  zero  so  long  as  the  constant  p  is  equal  to  GJL , 
whore  G  is  the  universal  gravitational  constant  and  Ik  is  the  total  interior 
mass;  thus  J  may  bo  eliminated  by  replacing  p,  if  necessary,  by  p  s  GJL , 

Finally,  may  bo  made  zero  by  taking  0  at  the  centre  of  the  interior  mass. 

However,  it  is  convenient  to  retain  all  terms  in  the  expansion  during  the 
general  development. 

Thus  the  form  assumed  for  the  disturbing  function  is  given  by 

00 

u  =  yv  ,  (12) 

n=-oo 

where 

Un  =  -  p  Jn  Rn  r"""1  Pe  (sin  3)  .  (13) 

4«2  Development  in  terms  of  elements  and  mean  anomaly 

Equation  (13)  gives  Un  as  a  function  of  r  and  3  (X  is  missing  because  of 
the  assumed  axial  symmetry)  at  all  points  of  a  certain  region  in  space  and 
hence  at  all  points  of  any  orbit  that  lies  entirely  within  this  region.  In 
order  to  integrate  Lagrange's  equations,  (2)  to  (7)>  it  is  necessary  to 
replace  this  function  of  r  and  3  by  a  function  of  the  elements  a,  e,  i,  Q,  u 
and  of  the  mean  anomaly  M  for  the  particular  orbit.  This  we  can  do  by 
expressing  r  and  3  as  functions  of  a,  e,  i,  fl,  w  and  M  by  the  usual  relations 
for  the  Kepler  ellipse.  Two  stages  are  required, 

m  i 

j 


The  centre  of  the 
necessarily  at  the 


(11) 


14 


( 


First  we  replace  Py  (sin  0)  by  a  function  of  i,  w  and  the  true  anomaly  v. 
3  4 

We  follow  Groves  with  minor  modifications  to  the  notation  and  considerable 
simplification  due  to  the  assumod  axial  symmetry. 

If  y  is  the  angle  measured  from  the  north  apex  of  the  osculating  orbit 
back  along  the  orbit  to  the  satellite  we  have 

Y  n  j  -  (w  +  v)  (14) 


by  definition  and 


sin  0  =  sin  i  co3  y 


(15) 


by  elementary  spherical  trigonometry.  V.'e  now  apply  the  addition  theorem  for 
Legendre  polynomials  which  is  given  by  equation  (Cl)  of  Appendix  C,  replacing 
n,  m,  <f> ,  01,  0  and  X  by  l,  k,  x/2  -  0,  %/2,  i  and  y  respectively.  This  gives 

l 

P4  (sin  p)  =  y\  (°)  *£  (cos  i)  00s  kr  ,  (16) 

k=0 

where  u^  -  1  if  1:  =  0,  and  =  2  if  k  *  Q. 

Now,  from  equation  (B2)  of  Appendix  B, 


k  (=03  i) 

P.  (eo3  i)  s  sin"  i - r  , 

d  (cos  i) 


(17) 


where  the  second  (k'th  derivative)  factor  is  a  polynomial  in  cos  i.  This  poly¬ 
nomial  does  not  vanish  when  cos  i  =  1  (unless  k  >  C  when  it  disappears  completely) 
and  in  fact 


(cos  i) 
d  (cos  i)k 


U+k) 


2*  (C- k) r.  k : 


when  i  =  0  . 


Hence  it  may  bo  normalized  -  in  a  certain  useful  sense  -  and  we  define,  for 

k  <  6, 


ok  / ,  ,  '  |  ,  |  d^P  1  (cos  1) 
(i'l  _  2  C^-k )  ■  kl  _ _ _ 

' '  "  (u+k) 1  .  ,  .vk 

v  7  d  (cos  1) 


(18) 


vVe  now  combine  equations  (13),  (16),  (l7)  and  ( 1 S) ,  writing 


( 
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t 

\  -  £u*  -  <19> 

kaO 

whoro 

=  -ji  Jn  Rn  r~n-i  ^  ^  |0)  sinki  Ak  (i)  coa  kf  .  (20) 

2  k* 

The  function  A  (i),  which  is  1  when  i  =  0.  is  conveniently  expressed  as 

Kj 

a  function  of  cos  i  and  f,  where 


f  =  sin^i  , 

Tabic  1  givos  the  particular  functional  expressions,  for  C  and  k  up  to 
C  s  4,  k  s  4*  A  recurrence  relation  is  given  in  section  6,2  and  derived  in 
Appendix  D, 

In  the  second,  and  more  complicated,  stage  of  the  replacement  of  r  and  0 

by  a,  c,  i,  fi,  cj  and  U  we  have  to  deal  with  the  factors  r~  ”  and  cos  ky  which 

occur  in  the  expression  of  U  .  by  equation  (20),  It  is  usual^*^  at  this 

nx  22 

point  to  appeal  to  Hansen’ 3  X  function  ,  by  means  of  which  one  can  write 

OO 

r  n  1  cos  ky  =  ^  a  n  ^  Xgn  (e)  cos  [sM  +  k  (w-|x)]  (21 ) 

S=*oo 

(remembering  that  y  =  -  u  -  v).  Expression  in  terms  of  a,  e,  i,  0,  w  and 

M  is  clearly  complete  if  this  function  is  used. 

In  the  present  application,  however,  tills  approach  is  unnecessary,  though 
the  form  of  equation  (21 )  helps  to  clarify  the  following  argument.  We  have 
resolved  to  eliminate  short-period  terms  in  the  solution  of  the  Lagrangian 
equations,  by  integrating  them  (with  respect  to  time)  over  a  complete  revolution 
of  the  satellite.  But  dt  is  proportional  to  cLM,  by  equation  (8),  Then  considera¬ 
tion  of  Lagrange's  equations  (2)  to  (7),  in  conjunction  with  the  expression  for  U 
in  terms  of  M  using  equation  (21 ),  shows  that  it  is  legitimate,  and  certainly 
simpler,  just  to  integrate  U  with  respect  to  M  before  substituting  in  the 
equations. 

Now  if  (av)  denotes  the  average  valuo  of  over  a  complete 
revolution, 
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Unk  (av)  * 


2% 

k  f  « 


11 J- R"  Ve<0> 

2k  k 


2n 


—  sinki  Ak(i)  J  r”n~^  cos  k  (v+cj-gft)  dM 


....  (22) 


by  equation  (20).  But  for  the  Kepler  ellipse, 


•=  (1  -e2)“"  dv 


dM 


(from,  for  example,  Chapter  2  of  Ref. 23).  Hence,  using  the  basic  relation 

2 

p/r  =  1+e  cos  v,  where  p  =  a(l-e  ), 

2%  . 


,  -n+i  /• 

/  r”n~  cos  k  (v+cj-p^)  dM  =  — ^ - r— r  /  (l+e  cos  v)n-1  cos  k  (v+olw|tc)  dv 

J  ac(l-e2)^  J 

°  °  ....  (23) 


Kaula4  carries  out  this  last  integration  for  n  >  1,  but  the  integration  is 
equally  valid  if  n  <  0.  Using  the  binomial  theorem  for  any  index,  positive  or 
negative,  ard  the  relation 


s  o”s 

cos  v  =  2 


for  s  a  positive  integer,  we  get 
2 % 


O 

^(q)  cos  ^s”2q^  v 

q=0 


J  (l+e  cos  v)n  ^  cos  k  (v+co-^ft)  dv 

°  2%  on  s 

=  j  (f)  ^,(q)  cos  ^s_2q^  v  003  k  dv  * 


3=0 


q=o 


....  (2k) 


where  the  series  in  s 
zero  or  negative  by 


terminates  if  n  >  0,  and  is  defined  for  n  positive, 


n-1 


=  (n-1 )(n-2)  ....  (n-s)/sl 


17 


Now,  for  any  x, 


s 


^T'  cos  (s-2q)  v  cos  x  =  ^  [cos  [x+(s-2q)  vj  +  cos  |x-(s-2q)  vj] 

q=0  q=0 


>  .) 

q=0 


q  1  cos  [  x-(s-2q)  v} 


(25) 


since 


&  -  (A)  • 


Also 


oo  3 


s=0  q=0 


0<q<s  oo  oo 

■  £  -IE 

q,S  q=0  s=q 


From  equations  (2/*.),  (25)  and  (26),  with  x  =  k(v+a)-^x), 


2k 

j  (l+e  cos  v)n  ^  cos  k  (v-hj— gx)  dv 


2k 


oo  oo 


(26) 


-  / 1  !(■:')©  (ge)3  cos  [k  (u-gx)  +  (k-s+2q)  vj  dv  • 


q=0  s=q 


....  (27) 


But,  of  the  double  series  on  the  right  hand  side,  terms  for  which 
k-s+2q  *  0  integrate  to  zero.  Hence,  since  s  =  k+2q  satisfies  s  >  q, 
equations  (23)  and  (27)  give 


2k 


-n-1 


cos  k  (v+w— gft)  dM 


2k  p~n+‘*  /  n-1  \  /k+2q\ 

TfebWl » ) 

4=0 


(■ge)k+2<1  cos  k  (u-^x)  . 


•  •  •  • 


(28) 


m 
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Nov/  define 


00 


»*<•>  ■  Efelq)(k;2<i)^2q 

q=0 


(29) 


Then  from  equations  (22),  (28)  and  (29), 


U  ,  (av)  s  •  n  J  ( — 
nkx  '  n  Vp 


\  #°> 


22k  k] 


•(1-e2)2 


sinki  Ak(i)  6|  S^e)  cos  k  (w-£k)  • 


....  (30) 


The  series  given  by  equation  (2°)  terminates  if  n  >  1,  and  S^e)  is  then  a 


polynomial  in  e  ,  But  for  this  case  k  <  n  (since  k  <  £  and  £  =  n).  For  k  =  n  the 
polynomial  vanishes  identically  but  otherwise  it  ha3  a  non-zero  leading  (constant) 


term 


,  viz  (“)■  Hence  it  may  be  nomadized  and  we  define,  for  n  >  1, 


Bn«  =  SJc)!(nk)  if  ° 


<  k  <  n-1 


and 


y  (3D 


Bn(c)  s  0  (though  this  is  actually  arbitrary). 


n 


A  formula  for  snjc(e)>  v/hon  n  >  1»  in  terms  of  (l//l~e2)  i3  given  by  equation 

(El  3)  of  Appendix  E, 

For  n  <  0  the  series  is  infinite.  Fortunately  it  may  be  summed  in  terms 
of  the  Bk(e)  polynouaal  defined  for  n  >  1.  In  Appendix  E  it  is  shown  that 


S 


me 


M  = 


(32) 


for  n  <  0  and  0  <  k  <  £  (where  £  =  -n-1  if  n  <  0  aid  £  -  n  if  n  =  0). 

Since  it  is  desirable  to  retain  a  single  formula  for  U  .  (av)  we  seek  to 

°  nk 

combine  equations  (31 )  and  (32)  before  returning  to  equation  (30).  This  can 

n  2  — 

bo  done  by  taking  the  factors  (R/p)n  and  (l-e  )2,  from  equation  (30),  with 
Snic(e)  and  writing 


,n 


rY'(  i-¥8*W  =  ("i;1)  (f)"  ^-e2)6  b*(o)  < 


(33) 


whore 


b  =  p,  h  s  n,  e  s  ^  if  n  >  1 


and 


b  s  a,  h  =  -n+1,  e  =  0  if  n  <  0  • 


The  function  B  will  now  always  be  referred  to  with  lower  suffix  h  instead  of  n. 

1c 

Particular  expressions  for  B^(e),  up  to  h  =  5,  k  s  5,  arc  given  in  Table  2. 

A  recurrence  formula  is  given  in  section  6.2. 


Before  giving  the  final  expression  for  U^^av),  it  is  convenient  to 
combine  together  all  the  numerical  factors,  each  a  function  of  n  and  k  only. 

V/e  define 

cx  / n-l\ 

Cn  "  -  ,21c  „  (kj 

for  all  integral  n.  The  coefficient  is  given  in  Table  3  for  n  running  from 
-4  to  +4  and,  therefore,  as  far  as  k  =  4*  It  is  zero  if  C- k  is  odd. 


Equations  (30),  (33)  and  (34)  now  give 


Unk  <av) 


=  n 


J  a 
n 


-1 


r  BkJe) 


(i-°2)e  2? 


(i) 


sinki 


Ck  cos  k 
n 


(u~Jk) 


(35) 


and  this  is  the  final  result  in  the  development  of  U,  It  is  observed  that 
Unk(av)  is  a  function  of  a,  c,  i  and  w  only;  Q  and  M  do  not  appear. 

5  INTEGRATION  OF  THE  PLANETARY  EQUATIONS 

Short-period  terns  having  boon  eliminated  from  the  disturbing  function, 
correct  first-order  perturbations  to  the  elements  can  only  be  obtained  at 
intervals  of  one  revolution.  The  integration  of  Lagrange’s  equations  (2)  to  (7) 
over  a  complete  revolution,  from  node  to  node  or  from  perigee  to  perigee,  is 
now  quite  trivial  on  taking  U^av)  from  equations  (12),  (19)  and  (33)*  If 
A£,  in  the  notation  of  Merson^,  indicates  the  change  in  a  typical  element  X> 
over  one  revolution, 


2%  d£ 
n  dt 


(av) 


t 


whore  n,  hero,  i3  the  mean  motion. 
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Tho  (first-ordor)  oxprossion  for  each  can  bo  written  down  as  a  double 
sum  of  toms  with  indox  n  running  from  -oo  to  -k>o  and  k  from  0  to  £•  For 
simplicity,  however,  we  write 

n=-^o 

and  then,  as  in  equation  (19)> 

e 

=  ^*nk  * 
k=0 

As  a  further  simplification  we  introduce,  as  additional  £  .  ,  the  quantities 
and  given  by  equations  (44-)  and  (45)  which  follow. 

When  substituting  into  equations  (2)  to  (7),  it  must  be  remembered  that  b 
is  a  function  of  a  and,  when  n  >  1,  also  of  e. 

The  expressions  for  the  now  easily  be  written  down  a3  follows 


(36) 


(37) 


ank 

enk 

ink 

nnk 

Wnk 

V 

*nk 

Pnk 


0  ,  (38) 
k  ek-1  (l-e2)£+i?  sinki  A^(i)  Bk(e)  sin  k  (u-^x)  ,  (39) 
-  e(l-e2)-1  cot  i  e^  ,  (bO) 


Ck  ek 
n 

0-e2;6-i 

.  k-2. 
sin  i 

Dk(i)  Bk(e)  cos 
</  n 

k  (u~Jx)  , 

(41 ) 

~  n  . 

nk 

cos  i  + 

*nk  > 

(42) 

-  (l-e 

9  A- 

*  pnk  ■ 

t 

(43) 

Ck  ek- 
n 

2  (i-<=2F 

-e  k. 

sin  i 

Ak(i)  EK(e)  cos 
C  '  nv  ' 

k  (u-±n)  , 

(44) 

2(n+1 ) 

< <1- 

2Ne  .  k.  ,k/.x  _k/  N 
-e  )  sin  l  A, (i)  B,  (e)  < 
v  n 

SOS  k  (w--grc)  , 

i  (45) 
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k  k 

whore  D.^i)  and  E  (e)  are  defined  by 
v  n 

D^(i)  =  k  cos  i  A*(i)  +  sin  i  A^(i)  ,  (46) 

Ejj(e)  =  [k  +(2n-k-l)e2]  b£(s)  +  e(l-e2)  ^  B*(e)  if  n  >  1  (47) 

and 

E^(e)  =  k  B*(e)  ♦  o  B^(e)  if  n  <  0  .  (W) 

u’o  observe  that  E&O)  -  E^(o)  =  k.  Thus  in  particular  D?(0)  =  E°(0)  a  0 

v  n  v  n^ 

and  this  is  why  these  functions  can  not  be  normalised  in  the  way  that  A. (i) 
and  B^(e)  were.  The  function  D^(i)  is  given  in  Table  4  as  far  as  C  =  4,  k  a  4> 
and  Ek(e)  in  Table  5  for  n  running  from  -4  to  +4  and  a3  far  as  k  s  4. 

6  EXPLICIT  EXPRESSIONS  FOR  THE  ^ 

6. 1  The  special  cases  n  =  -1.  0  and  +1 

These  special  cases  have  been  mentioned  briefly  in  section  4.1  when  the 

final  expression  for  the  disturbing  potential  was  introduced.  It  may  now  be 

seen  from  equations  (38)  to  (45)  that  the  K,n  have  their  right  values  for  these 

cases.  Thus  all  £  are  zero  except  that 

n 


For  n  =  -1,  U_^ 

For  n  =  0,  U  = 
'  o 

full  potential  is 


Thus  the  perturbation  may  be  combined  with  the  central  force  itself,  using  - 
instead  of  p  -  an  effective  p  =  p(l-JQ).  An  orbit  in  the  field  of  force  may 
be  described  either  by  constant  ’barred*  elements  based  on  n  or  else  by  the 
normal  osculating  elements  based  on  p.  In  the  lattoi  case  we  must  tolerate 
variations  in  the  elements,  due  to  the  perturbation. 


cr  =  p  =  -2  , 

o  ro 


\i  J 


-1 


R 


.  This  is  constant  and  so  may  be  disregarded. 


^  J 


j  i.e,  with  no  other  disturbing  tern  present,  the 


p(l-J  ) 
*+U  =  — 21 

r  o  r 
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O.b 


Although  the  normal  osculating  elements  aro  not  all  constant  it  is  clear  - 
from  the  fact  that  the  orbit  really  is  an  ellipse  -  that  the  changes  in 
a,  e,  i,  f)  and  w  over  a  complete  revolution  must  be  zero,  Por  <x,  however,  the 
situation  is  different,  We  note  that  (varying)  a  is  related  to  (constant)  a 
by  the  equation 


»(M)  ■  =(M)  . 

since  both  sides  give  the  square  of  the  velocity.  Thus  the  osculating  mean  motion, 
n,  varies  during  a  revolution.  Since  the  la3t  equation  leads  to 


(49) 


it  follows  fairly  easily  that  the  mean  value  of  n,  the  'mean  mean  motion', 

exceeds  n  (where  n  =  p)  by  2  Jq  n,  to  the  first  order.  Then  equation  (8) 

leads  at  once  to  the  secular  change  given  by  cr  =  -2, 

0  |i  R  sin  P 

The  final  special  case  is  for  index  n  =  +1,  Then  U.  =  -  - 2 - 

r 

and,  with  no  other  disturbing  term  present,  the  full  potential,  p/r  +  U^,  is 
to  first  order  the  same  as 


P  (sin  P) 

q 


J L 


(r2  +  £  R2  +  2r  ^  R  sin  P)_i 


Hence,  using  equation  (Al)  of  Appendix  A,  the  potential  is  essentially  the  same 

as  for  an  inverse-square-lav/  force  directed  towards  the  point  r  =  -J^  R, 

P  =  instead  of  towards  the  origin.  It  follows  that  there  is  no  change,  to 

first  order,  in  any  of  the  osculating  elements  over  a  complete  revolution. 

If  the  special  point  r  =  -J^  R,  P  =  -gx  is  used  as  a  new  origin,  the  corresponding 

new  value  of  becoiries  zero.  That  the  new  origin  is  in  fact  the  centroid  of 

interior  mass  is  easily  seen  by  considering  equation  (A7)  of  Appendix  A  with 

the  first  moment,  /  r'  cos  0  '  =  C°,  equal  to  zero, 

(matter) 


6,2  Recurrence  relations 

The  functions  A,  B,  C,  D  and  E,  used  in  equations  (38)  to  (45)»  become 
so  complicated  for  numerically  large  values  of  the  index  n  that  it  is  very 
laborious  -  and  indeed  dangerous  -  to  ’work  out  their  expansions  from  first 
principles  each  time.  Now  the  functions  are  based  on  Legendre  functions,  for 
which  there  exist  well-known  recurrence  relations.  Hence  it  is  logical  to 
look  for  such  relations  for  the  functions  A,  B,  C,  D  and  E, 
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Recurrence  relations  arc  developed  in  Appendix  D,  The  complete  sot  of 
results  is  given  here.  It  should  be  noted  that  all  but  two  of  the  relations 
are  for  the  upper  suffix  (k)  fixed, 

A^(i)  For  given  k  >  0  and  for  oach  t  >  k  +  2, 


24-1 

C+k 


cos  i  A^_1 (i) 


C- k-1 

v+k 


kc-z^ 


(50) 


The  starting  values  are  given  by 

A^(i)  =  1  and  A^+1(i)  =  cos  i  . 

B^(e)  For  given  k  >  0  and  for  each  h  >  k  +  3, 

Bh<°>  ■  i^rBh-/e)  -  (1-°2K-2(e)  •  (51) 


The  starting  values  are  given  by 


BT  . (e)  =  1  and 

k+1 


k+2  (C )  ■  1  » 


in  addition  B^(e),  for  k  >  1,  may  be  (arbitrarily)  taken  as  zero, 
k  /  \ 

(i)  n  >  0.  For  given  k  >  0  and  for  each  n  >  k  +  3> 


pk  _  (n-1 )  (n-2)  (n+k-1 )  pk 

n  (n-k)  (n-k-1 )  (n-k-2)  n-2  * 


(52) 


also,  for  k  >  1 , 


Pk  (k+l)(2k+l)  \  pk-1 

k+2  =  «c2  vA1 


(53) 


with  as  definod  immediately  after  equation  (l6). 

The  starting  values  are  given  by 

Cg  =  j  and  cj^  =  0  for  k  >  0;  in  addition,  C°  =  -1  and  =  0  if  k  +0  , 


(ii)  n  <  0,  For  given  k  >  0  and  for  each  n  with  -n  >  k  +  3, 


n 


_  ("?*k-2)  (-n-fk-1 )  (-n+k) 

(-n)  (-n-1 )  (-n-k-1 )  n+2 


(54) 


2k. 


also,  for  k  *  1, 


'-k-1  =  "  2k(k+l)  “k_1  -k 


Tho  starting  values  arc  givon  by 


(55) 


c-k-2  -  0 


o.i  =  -1  . 


^(i) 


There  i3  no  throe-term  recurrence  relation  when  k  >  1. 


However  the  case  k  =  0  is  the  most  important  one  and  for  this  case,  for 
each  C  >  2, 


D°(i) 


2C-1  .  / .  \ 

—  003  1  De-i(l) 


The  starting  values  are  given  by 


D“(i)  =  0 


D°(i)  =  -f  , 


where  f  =  sin^i,  It  follows  immediately  that  the  leading  term  of  D°(i)  is 

V 

C(C+l)f  cos  i  if  0  is  even  and  £(C+1 )f  if  0  is  odd. 

Ei>> 

Again,  there  i3  no  three-term  recurrence  relation  unless  k  =  0;  for  this 
case,  and  for  each  n  >  2, 


E°(o) 


22lls°  (e) 

n-1  n-1 N  ' 


-  ^7  En-2<‘>  • 


The  starting  values  arc  given  by 


E  (e)  =  0 

o 


E°(e)  =  c‘ 


It  follows  that  the  leading  term  of  *^n(e)  io  ^  n(n+l)e  , 

For  n  <  0,  E°(e)  is  given  at  once  -  see  equation  (223)  of  Appendix  E 


E°(e) 

n 


=  2°  (e)  . 

-n-1 


(58) 
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6. 3  f or  n  >  2 

As  follows  from  section  4.1  the  terms  Ufi,  with  n  >  2,  of  the  disturbing 

function  relate  at  once  to  the  problem  of  satellite  motion  about  an  axi-syrrmetrio 

body.  In  the  case  of  the  earth,  the  coefficient  J_  is  of  order  10”^  and  J,, 

-6  ^  3 

J,  etc  are  of  order  10~  . 

4 

■j 

Merson  has  given  e  .  i  ,  ft  ,  w  and  cr  for  n  =  2  to  6  inclusive:  i.e. 
he  has  given  synthetic  expressions  for  the  £n  from  which  the  can  at  once 
be  picked  out.  However,  Ref.  1  gives  no  general  procedure  for  writing  down 


Such  a  general  procedure  is  now  available  from  the  recurrence  relations 
of  section  6.2.  These  relations  and  the  basic  equations  (38)  to  (45)  have  been 
used  to  obtain  the  complete  set  of  non-zero  £  ,  ,  extended  up  to  n  =  9,  which 

luC 

follows.  For  convenience,  i^,  and  arc  not  given,  since  they  may  at 
once  be  obtained  from  e^,  and  p^,  using  equations  (40),  (^2)  and  (43)* 

It  is  obsex-ved  that  non-zero  £  ,  can  only  occur  when  n-k  is  even  and 

0  <  k  <  n-2.  In  every  caso,  of  course,  a^  is  zero;  e^  -  and  hence  i^  - 

is  zero  for  k  =  0  (no  secular  terms  occur  for  the  eccentricity  or  inclination). 


P20 


S-S-i  °31 


p31 


S-sJt  % 


-  ^  (l-c2)  (l  f)  sin  i  cos  o, 

i 


cot  i  sin  w  , 


e~ 1  ( 1  +4  e2 )^1  -  ^  f  ^  sin  i  sin  u 
I2e  (l-e2)2  ^  f^  sin  i 


i  sin  u  . 


> 


iff,. 

-  ^  +  I  o2)  (i  -  5  f  +  ^  f2)  . 
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■}f  «  (1-e2)  f  sin  2co 

e2  ^1  -  ^  cos  i  cos  2  co  , 

^(1+2  e2)f  (1-^f)cos  Z' 


cos  tto 


-«  ■  -¥•"<«’>*  '('-I') 


C03  2  CO 


=  ^  (l-e2)  ^1  e2y  ^1  f  +  f2^  sin  i  cos  u  , 

.  .15.  e('1  +  le2W1.?ifti2Sf2\ootiiJin„  , 


=  -^e_1  (1  +  tLe2*fe4)(1-|f  +  f  fV"i  1»«  , 


=  -45  e(l-e2p '^1+^e2^1-If+|!-f2^ 


sm  i  sin  co  . 


rgp-  e2  (l-e2)  f  A  -  ^  sin  i  cos  3w  , 


=  -  sin  i  cos  i 


sin  3co  , 


-  4“p-  e(l+2e2)  f  ^1  - 1-  f ^  sin  i  sin  3co  , 

-  ^p-  (i_e2)2  f  -  g  sin  i  sin  }u 


n  =  6  =  -4^ri+5e2  +  :^e^>l  A-^f  +  44f2 


-ir^56 +-re^^“Jf+rf  jcosi  > 
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’62 


n62 


'62 


P62 


64 


n 


64 


* 


521  e2. 

32 


^  e(l-e2)  fl  +  j  s2^  f  ^1-3  f +  ^  f2^  sin  2o 

2if2\ 

■  ^(nfs2  +  2s^(i-3f^f2)=os2w  , 

=  Z|Z1  O2  (i-»2)*  ^1  +  i  e2^  f  ^1-3  f  f2^  003  2 u 


A  ♦  ^  s2)  ^1-6  f  +  2^  ffc  ]  cos  i  cos  2  w  t 


»3  (I-2)  f2 


64  =  512 


sin  4w  , 

\  ,v/  / 

■si? e4  f  ^  ■  20 f)  003  1  003  ltU  ’ 
15Zle2(l+Zs2)f2(1-^f)cos^  , 


p6,  ■  ^  '2  M f)  " 


n  =  7  e 


71 


Q 


71 


*■ 


71 


e 


Q 


73 


73 


*73 

p73 


=  “  ^  (l-e2)  (j  +  f  e2  +  8e4)  (1_4^f  +^f2"^f  f3)  3in  1  003  w> 

=  ^e(l+fe^fe^)(l^f  +  ^f2-^^coti3inw  , 


P71  = 


125.  - 

16  e 


a,.S2f*.i2,>W 


sin  l  sxn  w- 


105  e(l-e2)2  +  Z  e2  +  le^  -ZZ  f  +  51  f2  f3^  sln  1  sln  “  • 

■  -  W e2^-32)  (-l-=2H'-rf^f2) 3in  1 303  3" 

=  ^e3(1  +  8  e2)(1-ff  f  +  ^Lf2)sin  1  003  1  sin3“  > 

=  ^e(1*lte2-  e4)f  (1-4Lf+Mif2)3ini  3inJ"  * 

=  e3(l-e2)5  +  |  e2^  f  ^1  -  f  +  -^  f2^  sin  i  sin  lu  . 


jm 

if' 
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J 


(’-S') 


sin  i  cos  5“  , 


l&X  e5  f  (j  -|o  f)  sin  1  003  1  sin  5u  ’ 

l&X  ®3  0  +  I  g2)  f2  0  -12  f)  Sin  1  Sin  5w 

^  e^(l-t'2)2  f2  sin  1  Ein  5“ 


n  =  8  nc 


-  ^(-i1  ^2*f  0  - 18  f2-^ f3t^ ft)  ■ 

-  S£  <'-=¥  (i  *  §.  sM2i  e\  If  ^  (,  .  ,»  ,.22  f2-  ft  f3 

*^)* 

-^e(’-e2)  (l*f=24=4)f  (l-f  f^f2-^f5)si„2U  , 


e(,-e2)  (l  *f  e2,i  /)  f  (l  -f  f  ^  ^  ^  ^ 


p64  ~ 


"^e2  (1  +  3  e2  +  ^e4)(1-11  f  +  ^|2f2-181  f3)  cos  1  003  2“  » 

-  W'M* 

cos  2  u  , 

-  2Ws2"ii  (1*¥e2*!o e4) f2  (1_  ¥f*  3^f2)c°s  . 

-  ^  «4(-=2)i  o *y <z  o- -  f f2) 


29 


86 


86 


'86  = 


p86  = 


-^7i8e5(l"e2)f3(1’l4f)sin6u  ' 
"  327§e6f2  (1-yf)cos  ico36"  , 

-g§e4(-l«2 


32768  ”  ^  +  2  e  )  f  ^  "  14  f)  003  6  “  * 


cos  o 


IL=  .3  e 


91 


n 


91 


_  .  Ill 

32 


0*r 


e2+ile4  +  ii/ 

6  +  8  6  +  64 


) 


'91 


x  ^1-33  f  ♦  f2  -  f3  *  f4^  cot  i  sin  u  , 


A -11  f*^f2-4£  sin  i  Sin  u 


91 


16  ‘  *  128  ‘  *)  3i: 

,2  «' 


(in 


sin  u 


'93 


fi 


93 


'93 


x  (^f  f+i?f2-irf3) sinicos3u  * 

_2^  3/  i  2  3  4\ 

128  e  y  +4  e  +16  e  ) 

x  f+^f2-|^f3)sinic0s  i  sin; 

.^e(1+^e2  +  ge4+|e6)f 

x  O-f  f+jiFf2-lrf3)sini sin3w  » 


30 


93 


'95 


95 


p95 


95 


'97 


97 


'97 


97 


><  ('- £if*^r2-irfi)3ini  sin  3u  • 

^.Ve2)(wle2)r2(l-|^i^)^ 

-^=5(,^e2)fHf*^f2)sini 


in  i  cos  5^  , 
i  cos  i  sin  5“  , 
sin  i  sin  5w  , 


1024 


(l-e2)"*  ^1  +  |-e2jf2  ^1  ~  \  f  +  ff  f2^  sin  1  sin  5w 


i  cos  7  u  , 


ifaieW) 

”  327b8  e7  f 2  ^1  - sin  i  oos  i  sin  7  u  , 

-  e7(l-e2)-’  f3  (l- jg  sin  1  sln  7“ 


in  i  sin  7^  , 


6.4  for  n  <  -2 

The  recurrence  relations  of  section  6.2  and  the  basic  equations  (38)  to 
(45)  have  been  used  to  obtain  the  non-zero  extended  for  negative  n  as 

far  as  n  =  -6.  Application  of  the  results  will  be  considered  in  section  8  for 
n  *;  -3  and  in  section  9  for  n  =  -2.  The  rernarks  of  section  6.3,  concerning  the 
listing  of  results,  apply. 


n  =  -2 


'-21 


3  2  »— 

“  2  (l“c  '  sin  003  “  >  0_o1 

3—1  2  — 

r?  e”  (l-e  )2  sin  i  sin  u  ,  P  21] 


3  2  — ~ 

^  e(l-e  )  2cot  i  sin  u> 
-3e  sin  i  sin  u  . 


013 


31 


e  v  =  ^-e(l-e2)2  f  sin  2  u  ,  Q  e2(l-e2)  2  cos  i  cos  2uf 


1^22  =  ^  (1-e2)2  f  cos  2u  , 


P-32  =  "  2 


15  2, 

=  -  “  e  f  cos  2  w  . 


=  ^  (l-e2)2  ^1  +  ^  e2^  sin  i  cos  u>  , 

=  -  15-  e(l-e2)"^  ^1  +  £  e2^  (j  -  ^  ^  cot  i  sin  u 
=  -  ^  e”1  (l-e2)^  ^1+2  e2^  (jl  -  J  sin  i  sin  u 
=  e  ^1  +  4  e2^  ^  "  4  f  ^  sin  ■*■  3^n  u  • 


sin  i  sin  u  . 


=  e2(l-e2)2  f  sin  i  cos  }  u  , 

=  -  e^(l-e2)~2  cos  i  sin  i  sin  3  w  > 

=  -  *j?p-  e(l-e2)2  f  sin  i  sin  3  u  , 
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e  f  sin  i  sin  }  u  , 


n-  ?  =5  n«3o  =  ^  ^“e2)  2  (1+5«2+  8^«4)(l“  4f)cos  i  , 

^-50  =  -?0-2)i(l^«2)(l-5f*?'2)  » 

P-50  -  3(u5e^e^)(l-5^ff2)  . 

®-52  =  "  <l^'e(1'e2)^  (j  +2®2)  f  (1"^f)  8in  2“  » 

n-52  "  "  ^  e2(1_e2)_¥  ^  +  j  ^  “  3  f)  008  1  008  2u  * 
♦-52  =  -  ^  (1-e2)2  (1+e2)  f  ^1 cos  2u  , 

P-52  =  ^®2  (Uje2)f  (l'Jf)  003  2“  • 


I 


32 


'-54 


n 


2205  3/„  „2  .  , 

=  ~  ■'./  e  (1-e  ) 2  f  sin  4u  , 


-54 


128 

220 


if  2  _ ^ 

^20  e  Al-e  )  2  f  cos  i  cos  4<*>  , 


’Ley.  =  -  e2(l-e2)2  f2  cos  4u  , 


-54 

P-54 


64 


i  2 

e  +  f  cos  4“  . 


n  = 


—6  e 


-61 


1  cos  u 


n 


-61 


'-61 


P-61  = 

e-63  = 
n-63  " 
*-63  = 
P-63  = 

e-65  = 


0 


-65 


'-65 


P-65  * 


^  e(1-e2)-i  (l  *  f  =2  ♦  |  e1*)  (,  .  f  f  .  ^  f2)  cot  i  tin  .  , 

Dfp  e_1(l-t;2)2  ^1  +  p-  e2  + -p  "  2  f +  ^’  f2)  Sin  1  Ein  u  > 

-  e  ^  +  f  e2  +  f  (^1  -  f  +  ^r-  f ' ^  sin  i  sin  u  . 

-  e‘"(l-e2)2  ^1  +  0  e2^  f  ^1  -  f ^  sin  i  co3  3co  , 

~|5.  c3(l-e2)"^  ^1  + 1  e2^j  ^1  -  4p  cos  i  sin  i  sin  3«  , 

e(l-e2p  ^1  +  |  e2^  f  ^1  - 1  fj  sin  i  sin  3u  , 

-  ^-p  e3  ^  +  3  ^  “  0  sin  1  sin  3w  . 

-  26^'  c^(1~°2)2  f2  sin  i  c°s  5“  , 

72765  5  2  — 1 

2048  e',(l“e  )  2  f  cos  i  sin  i  sin  5“  > 
e3(l-e2)'-J  f2  sin  i  sin  5«  , 

-  ~^Q2k  ®5  f2  Sin  1  Sin  5“  * 


THE  ORBITAL  PERIOD 


The  period  of  an  orbit  is  perhaps  its  most  fundamental  characteristic.  A 
good  estimate  of  it  can  be  made  by  use  of  equipment  no  more  complicated  than 
a  clock.  It  is  desirable,  therefore,  to  consider  the  effect  on  the  orbital 
period  of  the  perturbations  considered  in  this  paper. 
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For  an  unperturbed  orbit  there  is  a  unique  period  given  by 


T 


y 


where  n,  the  mean  motion,  relates  to  a,  the  semi-major  axis,  by 

2  3 

n  a  =  p 


(59) 


For  an  orbit  which  is  changing  all  the  time  the  period  mu3t  be  carefully 
defined;  different  definitions  lead  to  different  quantities  associated  with 
the  same  orbit,  though  not  so  very  different  if  the  perturbations  are  small. 

The  two  most  useful  definitions  lead  to  the  anomalistic  period,  the  time 
between  successive  passages  through  perigee,  and  the  nodal  (or  draconic)  period, 
the  time  between  successive  passages  through  the  ascending  node.  We  are  here 
mainly  interested  in  the  latter,  T^,  but  a  forme la  for  the  former  is  given  at 
the  end  of  the  section. 

Following  equation  (59),  one  would  like  to  relate  to  n^,  the  value 
of  at  the  node;  the  latter  does  not  possess  a  first  order  variation  over  a 
conplete  orbit,  since  Aa  =  0.  We  write 

*  -  f(i-E,“©B*») '  (60) 

n=-«o 

where  b  is  as  defined  in  section  4.2  and  expressions  for  the  Tn  are  to  be 
evaluated, 

A  good  treatment  to  follow  is  that  of  Person  ,  who  observes  that,  since 
the  kinetic  energy  per  unit  mass  of  the  satellite  is  p  (l/t  -  l/2a)  and  the 
potential  energy  is  the  negative  of  the  potential  p/r  +  U,  the  equation  of 
energy  gives 

-  ^  -  U  =  const,  (6l) 

even  though  both  a  and  U  vary,  i/e  define  a'  from 


and  then  a1  is  constant  to  any  -  not  merely  first  -  order.  If,  now,  n’  is 
2  3 

defined  from  n'  a'  =  p,  n'  is  also  completely  constant  and  is  given  by 
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n'  =  n^1  +  +  0(U2)  .  (63) 

In  addition  to  n*  and  a’,  Mcrson  introduces  an  auxiliary  element  9  of 
which  tho  change,  A9,  over  a  complete  revolution  is  useful*  We  do  not  need 
9  since  we  alroady  have  p  in  section  5  and  A9  =  -  Ap* 

The  procedure  followed  is  to  write  down  the  first-order  cliange  over  one 
revolution  -  node  to  node  -  of  each  side  of  equation  (8);  thus 

T 

ro 

AM  =  Ac  +  /  n  dt  *  (64) 

0 

Introducing  the  constant  n1  by  means  of  equation  (63),  we  have 

T 

/n  ,  Tr 

AM  =  AC  +  n*  Tn  -  /  ISaJL  at  .  (65) 

o 

Now  the  change  in  M  over  a  nodal  revolution  can  be  related  to  the  changes, 
Ae  and  Aw,  in  e  and  w,  using  eccentric  and  truo  anomaly  as  intermediate 
variables  and  starting  from  the  relation  Av  +  Aw  =  2%,  Taking  Merson's  first- 
order  terms,  we  have 


AM  = 


251 


Aw  + 


sin  w  Ae  , 


•  •••  (66) 


where 


JL  _ 


=  1  +  e  cos  w 


Q 


The  integral  in  equation  (65)  can  be  expressed  in  terms  of  the  p^  of 


section  5;  these  p^  arose  on  integrating  equation  (7)  and  in  fact 


Q  , 


_2_  £U 

na  da 


dt  =  Ap  =  2-x 


(67) 


n=-oc 


where  n  denotes  mean  motion  on  the  L.H.S.  and  an  indicial  suffix  on  the  R.H. S. 
But  is  proportional  to  a_n~”'  and  so  we  must  have 


3  na  U 


o 


dt 


(68) 


n=-oo 


I 
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with  the  sane  re; nark  about  the  dual  uoe  of  n.  We  note  that  on  the  R.H.S.  of 

the  above  equation  pn/n+l  gives  o/o  when  n  =  -1.  Since  U_^  =  -  p  J_,|/R  we  L 

must  conventionally  take  p  /0  =  -2. 

By  now,  from  equation  (65),  we  can  express  n'T^  in  terms  of  the  en,  un, 
o"n  and  p^;  the  p^  can  be  eliminated  in  favour  of  Qn  since,  by  equations  (42) 
and  (43), 

p  =  c*  +  (l-e  )2  (cj  +  n  cos  i)  •  (69) 

rn  n  '  '  x  n  n  '  '  ' 


To  introduce  the  T  ,  which  we  are  trying  to  evaluate,  we  observe  that  from 
equation  (63),  on  substituting  for  from  equation  (l3)> 


n 


*  _ 


00  i  t  nn 

/  o  3a  J  R  > 

•"('-Ztf-'*10', 


n=-oo 


Hence,  frou.  equation  (60), 


n*  =  2k  (  1 


[l  '  EJ"©n(Tn  +  ^^(0))3  *  (?0) 


n=-oo 


Thus,  combining  equations  (65),  (66),  (68)  and  (70),  and  picking  out  the 
term3  involving  2k  J  (R/b)n  for  each  n,  we  have 


T  =  <r 
n 


.  -  jjStt  •  <«’>“  (t?  ■, 

-  "T  (’  *  'p)  sin  “  %  ‘  '  (7,) 

where  pn  is  given  by  equation  (69). 

It  is  remarked  that  for  n  =  -1  and  n  =  0  results  are  obtained  which  agree 
with  those  of  section  6,1.  'When  n  =  -1,  taking  p  ^/O  =  -2  as  earlier  and 
Pq(0)  =  1,  we  get  T  =  0;  so,  neglecting  other  disturbing  terms, 

T^  =  2k/d^,  though  we  note  that  n^  *  n'.  When  n  =  0,  <Tq  =  PQ  =  -2  and  so 


T  Sif,  -  J  (1  -*N|  ; 
n  L  0  \  rn/J 


but  from  equation  (49)  of  section  6.1,  a  =  a^  [l  +  Jq  (2a/r^  -  l)j  and  so  if 
n2  a5  =  p  =  p  (1  -  Jq),  n  =  [l  +  Jq  (1  -  3&/rn)J ;  hence  TQ  =  2tyn  as 
expected;  we  note  that  n  =  n*  (l  +  J  )  so  that  *  n*  *  n# 


<L 

£' 
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For  completeness  we  give  the  formula  for  tho  anomalistic  period,  T  ,  which 
corresponds  to  equation  (71 )•  If,  similar  to  equation  (60), 

T“  *  ’ 


n=-eo 


then 


where  r^  is  the  perigee  distance  and  (3^  is  its  latitude, 

8  APPLICATION  TO  LUNI-SOLAR  (GRAVITATIONAL)  Pr^TURBATION  THEORY 


8, 1  Introductory  remarks 

The  (disturbing)  gravitational  field,  produced  by  a  distant  body  -  such  as 
the  sun  or  moon  -  considerod  as  a  point  mass,  is  axi-symmetric,  the  axis  of 
symmetry  being  simply  the  line  from  the  main  centre  of  force  -  the  earth  -  to 
the  distant  body.  Hence  it  is  to  be  expected  that  a  first-order  account  of 
the  effects  of  the  disturbing  field  will  be  given  by  the  formulae  of  this  paper; 
it  is  only  necessary  to  interpret  the  values  of  the  J  coefficients  appropriately 
and  this  is  done  in  section  8,2, 


A  complication  arises  in  that  the  axis  of  symmetry  of  the  gravitational 
field  has  hitherto  been  assumed  to  lie  in  the  direction  of  the  north  pole,  7 /hen 
it  lies  in  another  direction  -  towards  the  distant  body  -  the  formulae  are  only 
immediately  applicable  if  the  ’orientating’  orbital  elements  -  i,  0  and  u  -  are 
assumed  to  relate  to  this  other  direction  (and  a  plane  perpendicular  to  it) 
instead  of  to  the  polar  axis  (and  equator).  The  formulae  must  be  modified  if 
i,  Q  and  u>  are  to  be  used  in  the  normal  way  and  this  modification  -  essentially 
a  rotation  of  co-ordinates  -  is  considered  in  section  8,3. 


In  this  application  of  the  basic  formulae  of  the  paper  it  is  first  assumed 
that  the  gravitational  field  does  not  vary,  i,  e,  that  the  disturbing  body  is 
stationary.  If  account  is  to  be  paid  to  tho  motion  of  the  disturbing  body  in 
its  own  orbit  relative  to  the  main  centre  of  force,  the  results  for  a  stationary 
body  must  be  averaged  with  respect  to  the  body’s  mean  anomaly.  This  subject  is 
discussed  in  section  8,4,  Yet  another  axis  is  suggested,  one  that  is  perpendi¬ 
cular  to  the  orbit  of  the  disturbing  body,  but  the  average  field  is  not  symmetric 
with  respect  to  this  axis  unless  the  orbit  of  the  disturbing  body  is  circular. 
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Consideration  of  long-term  secular  perturbations,  derived  by  a  further 
averaging  process,  is  postponed  until  section  10.3. 

8,2  Interpretation  the  J  coefficients 

The  disturbing  body  is  itself  a  centre  of  force  and  wo  suppose  that  the 
strength  of  its  field  i3  given  by  the  parameter  p^,  taking  the  local  potential 
to  be  p^/r^  where  r^  is  the  distance  from  the  body.  However  this  is  not  the 
disturbing  potential  at  a  satellite  in  orbit  about  the  main  centre  of  force, 
since  the  main  centre  is  itself  attracted  towards  the  disturbing  centre. 

Taking  r  and  R  to  bo  the  distances  of  the  satellite  and  the  disturbing 
body,  respectively,  from  the  main  centre  of  force  and  introducing  vector 
rotation  based  on  Pig.  2,  we  have 


the  accelox-ations  of  the  satellite  and  main  centre,  respectively,  towards  the 
disturbing  body  are  then 


^d%L 


and 


It  may  be  seen  at  once,  on  varying  r,  that  these  accelerations  correspond  to 
potential  terms 


and 


p^  r  cos  <p 


respectively,  where  <p  is  the  angle  between  r  and  R,  Hence  the  disturbing 
potential,  U,  is  given  by 


U 


M _ r  cos  (£ 

Vrd  R2 


) 


How  from  Appendix  A, 


00 


_L 


•I? 


+1  Pg(cos  (*>) 


'l 


n 


V003 


e=o 


n=-1 


so  that 


where  the  summation  covers  n  =  -1  and  -3,  -4,  •••«  “°°* 


(72) 


m 
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C 


Hence  tho  disturbing  potential  may  be  represented  by  equation  (ll)  with  the 
interpretation:  Jn  =  0  if  n  =  -2  and  0,  1,  . ...,  oo;  Jn  =  -  if  n  =  -1  and 

-3,  -4,  mm,  -oo;  3  is  the  distance  from  the  centre  of  force  (earth)  to  the 
disturbing  body. 

Tho  term  with  n  =  -1  is  constant  and  may  be  disregarded  as  stated  in 
section  6.1,  The  terms  with  n  =  -3,  -4  etc  are  associated  with  Legendre  poly¬ 
nomials  of  order  2,  3  etc  in  the  normal  way. 

8, 3  Perturbations  caused  by  a  stationary  disturbing  body 

dith  the  interpretation  of  the  Jn  just  made,  tho  perturbations  of  the 
elements  are  given  at  once  by  equations  (36)  to  (45)  if  i,  0  and  w  are  taken  as 
olements  relative  to  the  pseudc-ouuator ial  plane  for  which  the  'north  pole1  is 
in  the  direction  of  the  disturbing  body3  Since  v;c  prefer  to  have  tho  perturbations 
of  the  usual  elements  -  with  i,  Q  and  w  relative  to  the  true  equatorial  plane  - 
it  is  convenient  to  proceed  by  first  adding  dashes  to  i,  0  and  cj  when  used  as 
above,  so  that  the  perturbations  given  by  equations  (38)  to  (45)  are  to  a,  e, 
i',  O',  w'  and  O’;  they  involve  i'  and  o'  in  their  expression. 

To  convert  to  formulae  involving  tho  usual  elements,  relations  between  the 

dashed  and  undashed  elements  are  required.  They  may  be  obtained  by  means  of  tho 

axis  transformations  of  Appendix  ?.  To  give  the  results  here  wc  make  use  of  a 

system  of  axes  Oxyz  defined  as  follows:  Ox  is  towards  the  (true)  ascending  node 

of  tho  satellite,  Qy  is  towards  its  apex  and  0z  is  normal  to  tho  orbital  plane. 

Since  fixed  axes  are  required,  this  definition  must  be  regarded  as  strictly  valid 

2 

at  one  instant  only,  due  to  the  variation  of  i  ard  0,  With  the  notation  of  Cook  , 
wo  take  the  direction  cosines  of  the  disturbing  body  to  be  A,  B  and  C  in  the 
Oxyz  system. 

Now,  although  i'  is  at  once  defined  as  the  angle  between  the  pseudo-equator 
and  the  orbital  plane  (the  xy-plane),  O'  is  not  defined,  since  its  origin  is 
arbitrary,  Tho  transformations  involve  a  (instead  of  O'),  whore  a  is  the  angle 
in  the  xy-plane  from  Ox  (the  direction  of  the  true  node)  towards  tho  pseudo-node. 
Then  i'  ard  w*  may  be  elin.ine.tod  in  favour  of  A,  B  and  0  by  means  of  the 
relations 


w'  =  a)  -  a  , 


(73) 


which  is  obvious,  ard 

(-sin  i'  sin  a,  sin  i'  cos  a,  co3  i')  =  (A,  3,  C)  ,  (74) 


proved  in  Appendix  F, 
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Appendix  P  also  gives  the  required  expressions  for  the  perturbations 
Ai,  AQ  and  Acj  in  terms  of  Ai‘,  AfJ*  and  Aw1,  The  expressions  are: 

Ai  =  ■-.■fl-ry  Ai1  +  A  Aft*  ,  (75) 

sin  i*  9  '  ' 

sin  i  AO  =  -  -r~r  Ai*  +  B  AO'  (76) 

sin  i'  '  ' 

and 

Aw  +  An  cos  i  =  Aw*  +  An*  cos  i1  .  (77) 

We  now  have  all  the  relations  necessary  for  the  conversion  of  the  dashed 
versions  of  equations  (39)  to  (45);  equation  (38)  is,  of  course,  unaffected. 
Prom  equation  (39)  v/e  have,  since  e  =  0, 

enk  =  k  Cn  ek"1  (1-®2)2  Bfc(e)  sinki'  sink  (u'-£x)  .  (78) 

Now  Ak(i')  is  a  function  of  cos  i'  and  f*  (s  sin2i*)  so  that  from  equation  (74) 

V 

we  may  write 

A*(i>)  =  A*  (A2  +  B2,  C)  , 

and  similarly  D^(i')  =  (A2  +  B2,  C),  where  the  new  function  A*  (of  two 

O  V  v 

arguments)  is  a  polynomial.  Also,  from  equation  (73), 

cos  k  (w,--§x)  =  cos  k  (w-^x)  cos  ka  +  sin  k  (w-^x)  sin  ka  (79) 


and 


sink  (w’-^x)  =  sink  (w-^x)  cos  ka  -  cos  k  (w-^x)  sin  ka  ,  (80) 


where,  by  trigonometrical  expansion. 


cos  ka  = 


E  «-"’&) 

q(C*2q<k) 


k-2q_,  .  2q 

cos  sin  a 


(81) 


and 
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sin  ka  s  ^  (-l)8  ^2s+1^  cosk~2a"1a  a  •  (®2) 

s(0<2fl+1<k) 

On  substituting  from  equations  (79)  to  (82)  in  equation  (78)  ard  using  equation 
(74),  we  get 

*  k  <£  e*'1  (l-e2)4  ^(e)  Ak  (A2  +  B2,  0) 


j^sin  k  (u-4n)  ^  (-l)q  (^  A2q  Bk"2q 


+  COS 


k  Mu)  ^  (-1)3  (2ak+1) 


2s+1  k-2s-1 

A  B 


For  i  and  0  we  have  first,  from  equations  (40)  and  (41 )» 

4*n  .  -  k  ek0-e2)-i  *.)  °°3.1'  |(--).  sinki'  sin  k  (,.*) 

jm  i’  n  v  '  rr  '  f' 


^nk  s  °n  ek(1”e2)*^  B^(e)  -  sinki'  cos  k  (o)1-^)  . 

Applying  equations  (75),  (76)  and  (79)  to  (82),  these  lead  to 
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i*  =  -  <  eVe2)**  <(•) 


n 


ptofcJgfel  j"k  a£(a2  +  Bs(  c)  £  (_i;q ^  A2q  Bk'2q+1  c 

s  ^ 

j'k  a^(A2  +  B2,  C)  (-1)8  ^28  ^  A2 8+1  Bk"*28  C 

s 

-  d£(A2  .  B2,  0)  £  (-l)q  (kq)  A2q+1  Bk'2qJj 


cos  k  (w~ 

+  ~7 


• • • •  (84) 


and 


fl  .  sin  i 
nk 


<£  ek  (1-e2)  i  Bk(e) 


P3in2k  fk  Ak(A2  +  B2,  0)  V  (-1  )q  (kq^  A2q>1  Bk~Zq  C 

•-  *  q 

-  l£(A2  +  B2,  0)  ^  ("I)®  (2L)  A2”*1  S’1"29] 

S 

c  (u-4fl0  f,  .k/42  „2  /  >\8  /  k  ^  ,2s+2J 

■~^2  |k  Ae(A  +  B,C)^(-0  (^2s+iJA  B 


cos  k  ( w-i 


2s+2„k~2s-1, 


s 

T — » 


+  Dk(A2  +  B2,  0)^(-l)qV2q 

q 


k  N  A2q  Bk-2qt1  JJ  . 

....  (85) 


For  the  remaining  two  elements,  u>  and  c r,  we  avoid  highly  cumbersome 
expressions  by  using  -  as  in  section  5  -  the  auxiliary  quantities  \jr  and  p. 
From  equation  (77), 

A\jr  =  At* 

and  so  the  usual  analysis,  starting  from  equations  (44)  and  (45),  leads  to 


St 

J 

r 
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(“nk  +  nnk  008  1  *>  *nk  *  <£  ^  ^e>  ^  0) 


sk-2q 


[cos  k  (w-^rc)  ^  (-l)q  (j>q^  ^  B* 

-  sin  k  (w-|*)  (-1)S  ^+1)  a2&+1  Bk“2s”[j 


••••  (86) 


and 


(°nk  +  *„k  »)  Pnk 


=  2(n+1 )  C*  ek  B^(e)  Ak (A2  +  B2,  C) 


k-2q 


[cos  k  (w-2*)  ^  (-0q  k2q 

q 

-  sin  k  («-£*)  ^  (-1)3  ^k+1^  A23+1  B^23-1]  . 

s 

....  (87) 


For  use  in  the  next  section,  we  also  give  the  formula  for  the  potential 
component  U^,  for  n  <  -3.  From  equation  (35),  writing  =  -  p^/p, 


Uric(av)  =  -  Pa  cn  a_1"n  R"  ®h(e)  ®k  a«(a2  +  b2'  c) 


,k-2q 


X  [cos  k  (m-Jk)  Y  (_1)q  ^2q)  k2q  ^ 

q 

,ink(^)^(-l)a(2^)  A2-1 


-  sxn 


(88) 


V/hen  n  =  -3,  equations  (83)  to  (86)  lead  at  once  to  the  formulae  of 

2 

Cook  ,  on  multiplying  by  the  appropriate  factor  from  equation  (36),  viz 

2  2  2  O 

-2x  (a/&)  p^/p,  and  replacing  hy  Cook's  constant  K  and  a.J  by  p/n“,  where 

n  is  the  mean  motion.  Ref, 2  does  not  give  the  formula  for  A<r,  but  this  may 
be  obtained  from  equation  (87).  For  completeness  we  list  all  six  formulae 
(A  =  A  ^),  as  follows:- 
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Aa  =  0 


(89) 


Ae  =  - 


15*  Hde  (l-e2)2 

73~2 

rt  n 


|^A B  cos  2  w  -  ^  (A2  -  B2)  sin  2 <TJ : 

••••  (90) 


3*  u,  (1-e2)”2  2  2 

Ai  =  ■  ■  ■  ■; — w '  [A(2  +  3e  +  5e  003  2u) 

2  R  n  2 

+  5  B  e  sin  2w]  ,  (91 ) 


3*  |id  (l-e2)  2  2 

AO  =  - r — 2 - ”  [5  A  e  sin  2  w 

2R;  n  sin  i  2  2  1 

+  B(2  +  3e  ~5®  cos  2w)J  , 

....  (92) 


Aw+  cos  i  AQ 


3* 

2  R 


3 


2 

n 


[5  (2  AB  sin  2  w  +  (A2  -  B2)cos  2  wj 
-{2-  3(A2+B2)]] 


•  •  •  • 


(93) 


and 


Act  +  (l-e2)2  (Aw  +  cos  i  An)  =  — r — r  [(2+3  e2)[2-3(A2  +  B2)J 


R3n2 


-  15  e2  [2AB  sin  2w  +  (A2-B2)  cos  2 &$] 


m  m  •  m 


(94) 


The  author  has  written  a  Pegasus  computer  sub-routine  based  on 
formulae  (89)  to  (94).  One  point  which  is  worth  noting  is  that  the  sub-routine 
takes  these  formulae  exactly  as  they  stand  and  does  not  include  the  'second- 
order'  terms  (n  =  -4)  for  Aw  which  Cook  adds  to  equation  (93).  These  terms 
are,  in  fact,  no  more  important  for  the  element  w  than  for  any  other  element* 
They  appear  more  important  because  of  a  factor  e  in  the  denominator,  occurring 
for  Aw  only;  but  this  merely  reflects  the  indeterminacy  of  w  when  e  is  small. 

If  the  elements  e  and  w  are  replaced  by  elements  e  sin  w  aivi  e  cos  w,  then 
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•second  order'  terms  for  the  new  elements  are  of  the  same  order  of  magnitude  as 
similar  terms  for  the  renaining  elements.  This  must  bo  true  since 


and 


A(e  sin  u>)  =  sin  u  Ae  +  e  cos  w  Aw 


A(e  cos  u)  =  cos  u  Ae  -  e  sin  w  Au 


(95) 


Indeed,  if  e  is  small  and  u  large,  all  perturbations  to  e  and  u  must  be 
incorporated  indirectly,  instead  of  by  direct  addition  to  e  and  u.  Using 
e  sin  u  and  c  cos  u  as  intermediate  elements,  the  procedure  is  (i)  to  evaluate 
perturbations  to  these  elements,  given  Ae  and  Au,  (ii)  to  add  the  perturbations 
to  e  sin  u  and  e  cos  u  and  (iii)  to  calculate  the  new  e  and  u. 

Uhen  n  =  -4,  equations  (83)  to  (87)  lead  to  formulae  which  are  listed 

below  (A  =  A  .  ).  The  first  five  of  formulae  (96)  to  (lOl)  agree  with  results 

-4 

obtained  by  Allan  (end  of  section  4.2  in  Ref. 24),  who  uses  a  notation  by  which 
they  can  be  expressed  more  compactly.  Allan,  like  Cook,  does  not  give  an 
expression  for  Ac.  Our  results  are:- 


Aa  =  0 


(96) 


Ae  = 


15*  a(l-e^)2 

32  R4  n2 


[(4+3e2)(4-5A2-5B2)(A  sinu-  B  cos  u) 

.2  _2 


+  35e2  [B  cos  3u  OA^-B*)  +  A  sin  3u  (3B2-A2)]]  ,  (97) 


Ai  =  - 


15*  ae ( 1 -e2 )  2  „  ?  ? 

- r~~~2 - -  [(4+3e  )[10AB  sinu  +  (15A  +5B  -4)  cosuj 

32  Rtf  n  2 

+  35e  { 2  A  B  sin  3<*>  +  (A  -B  )  cos  3<*> I]  ,  (98) 


2  — 1 

15x  p.  ae ( 1 -e  )  2  2  ?  2  , 

AO  =  -  - ^—5 - 14  [(4+3e  ){10AB  cosu  +  (5A  +15B  -4)  sinu] 

32  K*~  n  sin  i  2  2  2 

-  3>s  j  2  A  B  00s  3^  -  (A  -B  )  sin  3<*>]] 


(99) 
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15*  ae*"\i-o2)2 

Au  +  cos  i  AO  s  - , — - - 

32  R4  n2 

x  [ (4+9e2)(4-5A2-5B2)(A  cosu  +  B  sinco) 

-  35e2  {B  sin3u  (3A2-B2)  -  A  cos  3w(3£2-A2)}] 

....  (100) 

and 

,  i  15^  ^  ae 

Atf  +  (l-e  )2  (Aw  +  cos  i  AO)  =  - ; - 5— 

16  R4  n 

x  [3(4+3e2)(5A2+5B2-4)(A  cosw  +  B  sinw) 

+  35e2  [B  sin3u  (3A2-B2)  -  A  cos  3w(3B2-A2)]] 

....  (101) 


When  n  =  -5,  equations  (83)  to  (87)  lead  to  quite  complicated  expressions. 

If  these  expressions  are  expanded  in  powers  of  c,  however,  the  leading  terms  are 

25 

fairly  simple  and  may  be  compared  with  results  obtained  by  Smith  .  Ref. 25, 
like  Refs. 2  and  24,  does  not  give  A  o',  but  we  get  this  from  equation  (87). 
Following  the  point  made  earlier,  in  connection  with  equation  (93)  and  Ref.2, 
wc  note  that  since  Ref.25  omits  e  terms  for  Ae,  all  terms  should  be  omitted 

for  Aw,  the  effective  error  being  then  only  0(e  ).  On  thi3  basis  we  list  the 

A  c  formulae  for  all  six  elements  as  follows:- 

-0 


Aa 


Ae 


Ai 


AO 


Aw  +  cos  i  AQ 


0  , 

0(e)  , 


15*  Pd  a 
2R5  n2 


AC 


15^  Pd  a 
5  2 

2  R3  n  sin  i 


+  0(e2)  , 

2)j  «■  0(e2) 


0(1) 


(102) 

(103) 

(104) 

(105) 

(106) 


and 


A<T  +  (l-e2)2  (Aw  +  cos  iAO)  =  £l  -5(A2+B2)  +  ^  (A2+B2)2J  +  0(e2) 


•  ••• 


(107) 
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Equations  (102)  to  ( 107)  agree  with  Smith's  results  though  there  is  an 
important  difference  in  notation:  the  direction  cosines  A,B,C  are  denoted  by 
C,A,B  respectively  in  Ref, 25*  It  is  important  to  note  that  there  is  no 
contradiction  in  giving,  in  equation  (107),  terms  of  an  order  which  are  neglected 
in  equation  (106)  :  Ao*  +  Aw,  when  e  is  small,  has  to  be  known  more  accurately 
than  Aw,  The  leading  terms  for  Aw,  which  we  have  argued  should  not  be  listed 
in  equation  (106),  are  given  by  Smith,  However  his  expression  contains  some 
errors  which  he  has  confirmed  in  a  private  communication,  Cu  working  out  the 
leading  terms  for  both  e  and  w,  we  get,  instead  of  (103)  and  (106), 

Ae  _  - —  >,  _X(  A  +B  )  [2  AB  cos  2w+(B  -A  )  sin  2wj  +  O(e^) 

8R5  n  l  b  J 


and 


Aw  +  cos  i  Afi  = 


15*  Pd  a 
2  R5  n2 


x|Jl-5(A2*B2)  +  &  (A2+B2)2j 

-  ~  I"1  -|(A2+B2)  {2  AB  sin  2w  +  (A2-B2)cos  2wJ  ^j+  0( 


e2)  . 


The  results  given  in  this  section,  for  n  =  -3,  -4  and  -5»  could  be  obtained 
directly  by  applying  equations  (73)  to  (77)  to  the  explicit  expressions  listed 
in  section  6.4,  instead  of  by  appeal  to  the  general  equations  (83)  to  (87). 

8, 4  Perturbations  caused  by  a  disturbing  body  in  a  Kepler  orbit 

In  allowing  for  the  motion  of  the  disturbing  body  -  with  orbit  assumed 
to  be  a  fixed  Kepler  ellipse  around  the  earth  -  we  are  for  the  first  time 
considering  a  potential  field  which  varies  with  time.  Over  an  integral  number 
of  revolutions  of  the  disturbing  body  the  difficulty  may  be  overcome  by 
averaging  with  respect  to  the  mean  anomaly  of  the  body,  i,e.  with  respect  to 
time.  If  the  Kepler  ellipse  has  zero  eccentricity  tho  average  field  is  then 
axi-symmetric.  In  general,  however,  it  is  not  axi- symmetric,  though  the  axis 
of  the  orbital  plane  of  the  disturbing  body  is  still  a  special  direction.  The 
average  field  may  be  thought  of  as  due  to  an  'ellipse  of  mass’  in  the  orbital 
plane  of  the  body. 

The  analysis  is  carried  out  in  two  parts.  In  the  first,  advantage  is 
taken  of  the  special  axis  even  when  it  is  not  an  axis  of  symmetry  -  the 


47 


disturbing  body  is  not  restricted  to  a  circular  orbit.  As  in  section  8.3, 
special  elements  may  be  introduced  -  this  time  we  use  double  dashejs  -  such  that 
i",  0"  and  w"  relate  to  the  orbit  of  the  disturbing  body  as  'equator'  while 
the  special  axis  points  to  the  'north  pole'. 

Now  general  formulae  for  all  the  elements  have  been  given  by  equations  (83) 
to  (87)  of  section  8.3*  It  would  be  perfectly  reasonable  to  use  these, 
introducing  the  double  dash  natation  and  interpreting  A,B  and  C  appropriately. 
However,  an  averaging  procedure  has  to  be  carried  out  and  it  is  simplest  just 
to  perform  this  on  the  potential  U^av)  given  by  equation  (88),  Lagrange's 
equations  may  then  be  set  up  in  terms  of  the  new  ^^(av,  av)  and  the  R.H.S, 
of  each  is  constant,  as  in  section 

In  the  second  stage  of  the  analysis  it  is  convenient  to  use  a  new 
notation.  Although  the  A£  notation,  for  each  element  £,  was  convenient  for  the 
change  in  £  when  d£/dt  was  integrated  over  a  revolution  of  the  satellite,  now 
that  a  second  integration  is  involved  it  is  preferable  to  give  expressions  for 
the  average  rates  of  change  of  the  £'s  relative  to  a  complete  revolution  of 
both  the  satellite  and  the  disturbing  body.  The  derivative  notation  is 
appropriate,  even  though  the  original  equations  have  been  doubly  integrated;  so 
the  average  rate  of  change  of  £  is  denoted  by  d£/dt. 

Writing  dZ/dt  =  T.  ^  d  ^^^/dt,  the  Lagrangian  equations  formed  in  the 
first  stage  of  the  analysis  give  d  2^/dt,  with  £  =  a,  e,  i" ,  fi",  u"  and  c r,  in 
terms  of  a,  e,  i",  fi"  and  u".  The  second  stage  consists  in  the  replacement, 
by  means  of  axis  transformations,  of  i",  0"  and  w"  by  the  usual  elements  i, 

0  and  u. 


We  start,  then,  by  averaging  U^av)  with  respect  to  the  mean  anomaly 
of  the  disturbing  body. 

Let  the  axis  system  Qx"y"z''be  similar  to  the  system  Oxyz  of  section  8.3, 
except  that  Ox”  i3  now  towards  the  pseudo-node  given  by  the  orbit  of  the 
disturbing  body.  Define  the  system  0X"Y"Z  by  this  orbit,  such  that  OX"  coincides 
with  Ox",  Let  a^,  c^,  and  refer  to  the  disturbing  body  and  its  orbit, 

and  let  0^  bo  the  angle  to  the  perigee  measured  from  OX",  It  is  clear,  from 
Fig.4  for  example,  that  0^  =  -  0",  measuring  0"  from  the  (fixed)  node  of 

the  disturbing  body.  Then  if  P  is  the  instantaneous  position  of  the  disturbing 
body,  the  direction  cosines  of  OP,  in  the  system  0x"y"^',  are  given  by 


(A,B,C)  =  (cos  Uj,  cos  i"  sin  u^,  -  sin  i"  sin  u^) 


> 


(108) 


as  may  be  seen  from  Pig* 5,  u,  being  (0,  +  v, )  and  i"  being  the  angle  by  which 

d.  a.  cL 

the  x"y"  plane  is  inclined  to  the  X"Y"  plane. 

Hence  from  equation  (88),  with  w"  in  place  of  u 


Unk(av) 


Jc  -1-n  _n  _k/  v  k  ,k  /  2  2  .  2  -411  \ 

=  -  nd  a  R  Bv(e)  e  kp  (cos  u,  +  co3  i"  sin  u,,  -  sm  i"  sin  u,) 


cos  k  (w"-^) 


2q  k-2q  k-2q 

cos  u,  cos  i"  sin  u, 
a  a 


-  sin 


q 

k  £  (-1  )5  (2sk+l) 


2  sfl  k-2s-1  .  k-2s- 

cos  u,  cos  1”  sin 

a 


’•«] 


....  (109) 


=  -  |i^  a~n"”^  Rn  B^(e)  P  (v^),  saY»  where  it  is  only  necessary  to  show 

one  argument  of  this  P,  since  u^  =  0^  +  and  other  quantities  are  constant. 
Then,  writing 


R  = 


ad  ^“*ed^ 
1+e^  cos 


and,  as  in  section  4,2, 

d  M, 
a 

it  follows  at  once  that 

_k  -1-n  k  _k/  \  n  2Nn+l|- 

e  ad  (W 

2% 

*  2n  f  F  ^vd^1+e  cos  vd^"n"2  dvd  *  (110) 

o 

Equation  (llO)  can  be  developed  analytically,  on  the  linos  of  section  4*2. 
However,  the  general  development  becomes  complex  and  we  do  not  embark  upon  it 
here.  Instead,  we  use  equation  (llO)  to  give  expressions  for  particular  cases 
of  n  and  k,  performing  each  integration  separately.  The  algebra  is  still 
somewhat  tedious,  but  the  majority  of  the  terms  disappear  on  integration.  The 
particular  cases  chosen  are  essentially  those  which  received  special  attention 


Unk(aV-  av)  =  - 
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in  section  8,3,  viz  n  =  -3,  for  both  k  =  0  and  k  =  2;  n  =  -4,  for  both  k  =  1 
and  k  =  3,  and  n  =  -5,  for  k  =  0  only  (v;o  diffor  horc  from  section  8.3  in 
giving  all  terms  ins  toad  of  leading  tcrmo  only). 

For  n  =  -3  and  k  =  0, 


P(vd) 


3  2  2  2 

=  1  -  yr  (cos  u,  +  cos  i"  sin  u,),  whore  u. 


ea  +  Ta  ' 


leading  to 


U-3,0^aV>  av) 


-  <-4rVi 

8  *1 


i2  (l+|  e2)  (1-3  c°s2i")  . 

....  (ill) 


?or  n  =  -3  and  k  =  2, 


2  2  2 

F(vd)  =  cos  2cj"  (cos  -  cos  i”  sin  u^)  +  2  sin  2  w"  cos  cosi"  sin 


leading  to 


U.3^2(av>  av) 


^-ea)'3/2 

16  ai 


a^  sin^i"  cos  2  w"  . 


(112) 


For  n  =  -if  and  k  =  1, 


P(vd)  = 


f .  5  ,  2  2  .2  0 

1-f  (cos  +cos  i"  sin  u^) 

x  (cos  u,  cos  u"  +  cos  i"  sin  u,  sin  w")  , 


leading  to 
U-4,1(av>  av>  = 


15nd  ead-e|)-5/2  3  ,  ,  2- 

- ; -  a  e  1  +  *•  e 

16  a4  V  4 


d 

in^i"^  cos  0^  cos  w"  +  ^1  -  ^  sin^i"^ 


sin^i"J  cos  i"  sin  0^ 

sin  . 
(113) 


•  ••• 


50 


018 


For  n  =  -4  ani  k  =  3> 


F^v^)  =  sin  3  w”  (3  cos*1  cos  i"  sin  -  cos^i"  sin^  u^) 


z  2  2 

+  cos  3  «"  (cos^  u,  -  3  cos  u  cos  i"  sin  u, )  , 


leading  to 


U-4,3(av’  av) 

525Ma  *d0"«i)"5/2 


256  a‘ 


3  3  2 

sin  i"  (cos  0^  cos  3w"  +  C03  i"  3in  0^  sin  3w") 

....  (114) 


For  n  =  -3  and  k  =  0, 


F(vd)  =  1-5  (cos2  +  cos2i"  sin2  u^)  +  (cos2  +  cos2!1'  sin2  u^)2 


leading  to 


U-5,0(av>  aT> 


n..  p2\-7/2 

d  1  d'  (l  ♦  5  e2 ♦  f  A  g  tI.|Vl-5  sin2i"  ♦  f  sinV) 


a 


64  at 


+  ^  e2  sin2i"  cos  20^  ~~6  s^n^"YI  • 


....  (115) 


For  each  of  these  cases,  the  first  part  of  the  analysis  may  be  completed 
by  writing  down  Lagrange's  planetary  equations  in  terms  of  U^^av,  av)*  These 
are,  in  effect,  already  solved,  since  every  term  on  each  R.H.S.  is  constant. 

But  in  the  second  part  of  the  analysis  the  usual  elements  are  to  be  restored. 
Since  we  aim  at  listing  formulae  for  d  ^^^/dt  at  the  erd  of  the  section,  it  is 
convenient  at  this  point  to  obtain  d  i"^/dt  and  d  fi^^dt  only.  For  the 
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former,  <3U(av,  av)/<D0"  is  involved,  so  that  it  is  necessary  to  replace  0^  by 
before  the  partial  differentiation* 

For  n  =  -3  and  k  =  0,  equations  (4),  (5)  and  (ill)  lead  to 


d  i" 
-3.0 

dt 


=  0 


(116) 


and 


d  n^.o 

dt 


3^a  (1  -4r5/2  (i-e2)_i  (1  +  li  e2) 


4  a 


3 


n 


cos  x"  • 


(117) 


For  n  =  -3  and  k  =  2,  equations  (i*.),  (5)  and  (112)  lead  to 


d  i" 

-3.2  _ 

dt 


15.d  d-e|)-3/2  eWri 


8  at 


n 


cos  i"  sin  i"  sin  2  w" 


(118) 


and 


d  ii»  15nd  (l-ed)'3/2  e2(l.e2)4 

— -  =  - - - 1 ‘ —  cos  i"  cos  2  u"  . 

dt  n  3  n 

o  a. 


(119) 


For  n  =  -4  and  k  =  1,  equations  (4),  (5)  and  ( 1 1 3 )  lead  to 


dt 


and 


15l*d  ed(l-e|r5/2ae(l-e2)~"  (l  +  Je2) 


32  a" 


n 


sin  l' 


5  cos  i"  cos  0^  sin  u"  +  sin2  iM^  sin  0^  cos 


(120) 


d  *-4.1 
dt 


ed(1-ed)"5/2ae(l-e2r"  (l  +  |e2) 


d  dv  A 

32  a* 


n 


x  ^3  cos  i"  oos  6^  cos  w"  +  ^17-4^  sin2  i"^  sin  0^  sin  •  ( 1 21 ) 


n 
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For  n  =  -4  and  k  =  3,  equations  (4),  (5)  and  (114)  lead  to 


- 

dt 


52^d  «d0 -ej)'5/2  ae3(l.e2^  3in  g 


256  a. 


n 


x  [(2-3  sin2i")  sin  0^  cos  3  w"  -  2  cos  i"  cos  0^  sin  3  w"j 


....  (122) 


and 


d  fl".  , 

_ 'rAil  _ 

dt 


525^d  ed(l-e^)-^2  ae3h_e2^ 

n 


256  a. 


2  1 

[2  cos  i"  cos  0^  cos  3  w"  +  (2-3  sin  i" )  sin  0^  sin  3  w"J  • 

....  (123) 


For  n  =  -5  and  k  =  0,  equations  (4),  (5)  and  ( 1 1 5 )  lead  to 


ii=5»£  _ 

dt 


135t*d  4(1-ea)"7/2a2(l-e2ri  (l  +  5e2  +  ¥ 


Ja  a1''  °a 

128  a' 


n 


x  sin  i"  sin  20. 


'  z  .  2 

1  -  T  sm  1" 


....  (124) 


and 


11^0  .  . 

dt 


45^  (l-el*7/2  a2(l-e2)_5? 


32  a; 


(l+5a2  +  ^  e4) 


cos  i" 


20. 


.7.5 
1  -  sin 


(125) 


\7e  now  give  the  relations  by  means  of  which  i",  Q" 
in  favour  of  the  usual  i,  n  and  w.  Let  i^  and  Q^,  like 
of  the  disturbing  body's  orbit;  they  are  referred,  like 


and  w"  may  be  eliminated 
etc,  denote  elements 
the  elements  i  and  0 
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of  the  satellite  orbit,  to  the  standard  (equator-equinox  based)  axes.  Let  a 
be,  as  in  section  8.3,  the  angle  in  the  orbital  plane  of  the  satollite,  from 
the  true  node  to  the  pseudo-node,  i.e, ,  to  the  intersection  with  the  orbital 
plane  of  the  disturbing  body. 

Then  it  is  shown  in  Appendix  G  that  di/dt,  dD/dt  and  dw/dt  are  given  in 
terms  of  the  double  dashed  derivatives  by 


d«3  .  df)  dw"  ...  dfi" 

dt  +  C0S1dt  =  dT  +  C031  IT  • 


C 128) 


Appendix  G  also  establishes  the  important  transformation  equations  which 
relate  the  axis  system  Oxyz  of  section  8.3  to  the  corresponding  system  Ox,y^z^ 
for  the  disturbing  body,  where  Ox,  points  towards  its  node.  In  matrix  form 
the  equations  may  be  combined  as: 


=  cos  (0-0^)  ,  =  cos  i^  sin  =  -sin  i^  sin  , 

=  -cos  i  sin  (Q-0^)  t  ^  =  sin  i^  sin  i  +  cos  i^  cos  i  cos  , 

N2  =  cos  i^  sin  i-  sin  i^  cos  i  cos  (Q-0^)  ,  =  sin  i  sin  , 

=  sin  i^  cos  i  -  cos  i^  sin  i  cos 

and 
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Some  of  the  elements  of  the  transformation  matrix  were  introduced  by 
25 

Smith  .  Indeed  L^,  1^,  L y  M^,  and  are  respectively  Smith's  Y y  Y^, 
Yg,  X^,  and  X^,  the  permutation  of  suffixes  being  in  conformity  with  the 
permutation  of  A,B  and  C  as  used  in  Ref. 25. 

Prom  these  transformation  equations  Appendix  G-  derives,  finally,  the 
relations  by  means  of  which  substitutions  may  be  made  for  i",  w",  a  and  0^. 
These  are 


u>" 

=  u  -  a  , 

(130) 

sin  i" 

sin  a 

=  “N1  * 

(131) 

sin  i" 

cos  a 

•h 

C\J 

S 

II 

(132) 

cos  i" 

II 

VX 

(153) 

sin  i" 

sin  0, 
a 

=  -  L,  cos  u,  -  M_  sin  ui, 

5  d  3  d 

(134) 

and 

sin  i"  cos  0^  =  sin  cos  .  ( 1 35 ) 

The  way  is  now  clear  to  giving  the  final  expressions  for  the  d  S^^/dt,  for 
the  particular  cases  chosen,  introducing  the  remaining  Lagrangian  equations  (2), 
(3),  (6)  and  (7).  Every  d  a ^^/dt  is  zero  and 

and  d  0 ^j^/dt,  di^/dt  and  dfi^/dt  in  equations  (ll6)  to  (l25)  are  used  through 
equations  (l26)  and  (127).  Equation  (128)  may  be  used  in  the  standard  way,  each 
side  being  d  \|r/dt,  in  the  notation  of  section  5»  The  quantities  L^,  ...., 
of  equation  (l29)  are  used  in  the  final  expressions  since  these  quantities 
depend  only  on  tne  orbital  parameters  of  the  satellite  and  the  disturbing  body. 

For  u  =  -3  and  k  =  0, 


may  be  forgotten.  For  d  i  .  /dt 


d  e 


dt 


=  0 


(136) 


d  i 


dt 


3Ka  (1-eJ 


2^3/2 

d- 


4  a|  n 


(l-e2)  2  +  £  e2 


)Ni  N; 


•  •  •  • 


(137) 
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d  (7 


=hO  =  _ 


3(id  (l-ej)'3/2 


dt 


4  at 


1^—  , 

n  sin  i  \  / 


....  (138) 


d  (J 


d  ff 


at 


cos  i 


dt 


(l-e2)~3/2  2  J-  2 

— - r -  (l‘e  )2  (1-3N2,) 

8  a3n  3 


(139) 


and 


d  a 


dt 


K,  (^e?)"372  2  2 

— - A -  (7+ 3  e2)  (1  -  5  N?)  .  (140) 

8  a3  n  3 


For  u  =  —3  and  k  =  2, 


d  e 


15(4,  (l-e2)"3/2  2  , 


-3.2  _  d  w  d 

8  ad  n 


dt 


e(l-e2)2  {(N2-!^)  sin  2o> 


+  2  N_j  Ng  cos  2<4  ,  (141) 


d  i  15n.(l-e2)  3/2  2  2  _i 

— -  - - - - =2 -  e  (1-e  )  2  N,  [N,  003  2  01+14,3^20)}, 

dt  8  a3  n  31  d 

....  ( 142 ) 


d  n 


=Jbl  _ 


15^  (l-e2)_3/2  2  2  1 


dt 


— - - - e^(i-e^)  2  N,  (N2cob2o>-N  sin2  oj)  , 

8  a^  n  sin  i 

....  (143) 


d  W-V.2  ^  ^-3  2 

~dt  *  <=»= 1  -sf* 


and 


15Ma  (i-4)-j/2  .  i 


8  a3  n 
d 


(1-e2)2  {(M^-hJ)  cos  2  ui 


-  2  ^  N2  sin  2  o)j  (144) 


d  o' 


Zl*2  =  - 


dt 


15nd  (l-e2)_3/2  2  f  2  2 

(l+e  )  {(Ng-tfj)  cos  2  o) 


8  a3  n 


-  2  Ng  sin  2  01}  . 


(145) 


a  i 
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For  n  s  -4  and  k  =  1, 


--.4,1 

dt 


15tid  ed(l-^)-5/2 


64  n 

a 


L(l-e2)2  £  e2 


x  jsinfiJ^  sinw  (itM^-5L^N2+1^^N^)  +cos  u^sinu  Ng+ISL^N^N^) 

-  sina^cosu  (ltl^+^L^fJ^+15^^2^^)  ~  cos  cos u>  (4i^-5^^+15iyt>Nj) } 

....  (146) 


d  T-4.1 
dt 


15k,  e,(l-e2r5/2 


d  dv'  d' 
32  n 


ae(1-e2)“i  A  +  £  e2' 


x  5  sin  sin  to  (f  t^N^N^+M^N^N^+M^N^N^)  +  5cosuj^  sinui  (L^N^N^+L^N^N^ 

♦  l3n,n2) 


sinud  cos  oj  l'^?2  +  ^  M3N2  -  TMy-IOL^  j 
cos  wd  cos  w  +  |p-  L^-TL^+IOi-^N^j 


(147) 


l=ital 

dt 


15kd  =dd-4)‘5/2 


32 


-  aed-e2)"*'  (l  +  £  e2") 

n  sin  i  \  / 


x  sinwd  sin 60  ^  ^  J^-TI^+IOL^IgJ 

+  cos  cod  sin  w  ^L3N2  +  ^  L3N^-7L3-ia'tlN2^ 
+  5  sin  u>d  cos  60  (l^N^+MgN^yMy^Nj,) 

+  5cos<^d  cos  w  (L1N2N3+L2N1N3+L3N1N2)  , 

m* 


(146) 
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-4.1 


+  COS  X 


-4,1 

dt 


15Md  d^1  V  -1,.  2^  .  9  2\ 

= - - - l - ae  (1-c  y  ll  +  £°  ) 

64  a^  n  \  / 

x  {sinw^sinw  (41^+5^^ +  15M^N^N^) 

+  cos<*>d  sin  u 

+  sin  cos  w  (41^-54^2+15^^^) 

+  COSW^COS  u  (41^+5^^2+154^^^)] 


(149' 


I--4,1 

dt 


15^ed('-4>'5/2  -1 


64  n 
d 


1  +  ^e2  +  f  e* 

4  4 


x  [sinw^  sinw  ( 41t> + 5L^ + 1 5Myi2N3 )  +  cosu^sinw  ( 442-5Myi^ + 1 ) 

+  sin cos  w  (4^^-51^2+151^^^)  +  cos  u>d  cos  w  (4L1  +5M^N2+1 54^ N^) ] . 

....  (150) 


For  n  =  -4  and  k  =  3, 


1575lid  ed(l-e2)-5/2 

- - -  ae  (1-e  )d 

256  aj  n 


x  [sinu^  sin 3  u  ( 31^^^  ^-4^1^)  -  cos  sin 3  w  ^y^+Ly^Ny^^N2) 
-  sinw^  cos  3»  OL^M^Nj+yy^)  +  cos  cos  3  w  (3M3N1-L3N2N3 


....(151) 


t 

U 
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d  1 


dfca 1 


dt 


525^,  «,d^)-5/2  j 


d  dx '  wd' 
256  ad  n 


;0-e2r* 


x  (2ainwdain3«  -  2oos  wdsin  3 « 

+  8lnudcos3w  (2M^3M^N^-3MjN^-2*L2N,j) 

-  cos  «d  cos  3  w  (2L3-3L^+3L3N^-A-M1N2))  ,  (152) 


d  ff 


-4^5 


dt 


„w>-» 


256  Sj  n  sin  1 
a 


x  [sin«d  sin  3  w  ( 2Td3+5M3N^-3M3N|-i4L2N1 ) 

-  cos  wd  sin  3  <*>  (2L3-3L3N^+3I<3N^*-4MjN2  ) 
-  2siiiudcos  3w  (L^-L^^ILN^I^) 

+  2  cosw^  cos  3  <*>  (M1N1-M2N2-3L3N1N2)) 


053) 


d  “-4.3  d  ^-4.3 

_t  +  cos  i  — 


256  ad  n 

x  |  sin  cod  sin  3  w  ( 31*^  +M3N2N'3+^2N^ ) 

-  cos  <*>d  sin  3  «  (3M3N1-L3N2N3-4L2N^) 

+  sinud  cos  3  w;. * ( 3L3N2-M3N^ Ntj ) 

-  cos  wd  cos  3  «  (3M3N2+L3N1N3+4L1N2)! 


(154) 


and 


d  <r 


=kil  . 

dt 


1575Md  ea(l-e2)-5/2 

- ; - ae(1+e  ) 

256  ad  n 


x  [  sin wd  sin  3  u  ( 3L^  +M-3N2N3+P^N^  ) 

-  cos  ud  sin  3  u  ( 3MjN^  -L^Ny-lfL^ ) 
+  sinwd  cos  3  «  (3L3N2-M3N1N3-pi1N2) 

(3M3N2+L3N1N3+2(L1I^)i 


-  COS  w  cos  3  W 
CL 


(155) 


( 
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For  n  a  -5  and  k  =  0  , 


d  *Ag  s  o 

dt  u 


(156) 


d  i 


-5.0 


dt 


135Md  ('-4)‘7/2_i2(i_e2)-i  ^5e2  +  ^e^ 


128  a^  n 


)  "lN3  (•  -H)+  ®a  (L1L3-M1M3  +  3  L3M3N2 

+  i  (Lj-^)  j  008  2  «d 

+  °d  [LiM3+L3M1+  3tL3W3N1N3+  I  N2(M3"L3}]  sin  2wd]  *  (l57) 


d  ff 


_ 


135H*  (l-e?)“7/2 


dt 


d_l_'  ~d; 

'4 


128  a ^  n  sin  i 


a2(l-o2H  ^1  +  5e2  +  ^e4 


) 


*  [0  + 1  ed)  Vi  i1  -  H)  +  ed  [w^M3  -  3  Wl 

♦  J  NjNj  (Lj-^)J  cos  2 

jL2iyL^2+^tL3M^2N3+  l  N/L^)!  8i"  2“d]  >  (158) 


+  e. 


d  w 


=5*0 


d  0 


dt 


+  cos  i 


_ 


dt 


135**d  (1_ed^  ?/,Z  2,  2,i/  3  2\ 

“1^ - ^  ’  T***) 


*  2  ed)  (1-10!1  +  f-^)  +  f  *d 

X  ( (L2-M2)  cos  2  o»d+  2Ly^  sin  2  «di^j  (159) 


and 


dt 


243m,  (1  -e^)'7/2 


256  a^  n 


2(1  +  if  e2+^et) 

[J  + 1  ed)  (1-1®23  *  f  "3)  +  3  ed  ^ 

x  RLj-M^)  oos2  «d+ 2LjMj  sin2u^!j  , 


(160) 


60 


Olfi 


25 

These  expressions  have  been  compared  with  those  given  by  3raith  as  far 
as  oomp&rison  is  possible  -  Ref* 25  neglects  oubio  and  higher  powers  of  a/a^,  0 
and  e^,  Most  of  the  terms  compared  agree  but  there  are  some  discrepancies  for 
du/dt,  as  there  were  for  Au  in  seotion  8.3,  and  one  discrepancy  for  each  of 
de/dt  and  dT/dt,  In  a  private  communication  Smith  has  agrood  that  the  errors 
are  in  Ref*  25. 

The  details  of  the  algebra  required  to  roach  tho  abovo  expressions  have 
boen  omitted.  The  analysis  is  fairly  straightforward,  but  one  point  of 
difficulty  is  worth  mentioning.  In  applying  equations  (l3l),  ( 1 32 ) ,  (l 34)  and 
( 1 35 )  to  eliminate  a  and  0^,  negative  powers  of  sin  i"  are  temporarily 
introduced.  They  must  be  cancelled  by  factorisation  of  tho  polynomial  expression 
in  L y  VLy  N^,  Ng  and  which  also  appears.  That  this  factorisation  is  possible 
is  not  always  evident  and  an  exariple  is  now  given  of  the  factorisation  for  one 
of  the  polynomial  expressions.  Tho  method  employed  is  to  eliminate  N^,  N2  and 
first,  using  relations  of  the  typo  =  -  (L^L^+M^M^)  and 

Ng  =  Ly^-LjM^,  then  to  regroup  terms,  to  factorise  and  fiaally  to  simplify 
by  re-introducing  ,  N2  and  N^. 

2  2  2 

As  our  example  we  extract  the  factor  sin  i"  (=  L,+M  )  from  the  expression 

3  3  ^ 

L3Wl*3  which  occurs,  multiplied  by  sinw^  sin  3  w/sin  i",  in  the  evaluation 

of  de  ,  ,/dt.  XJq  have 

=  L3 -l1M3 ) ( 1  -LI-,4 )  - 

=  -L1L3?.13(l|+m|)  ♦  M1  )  -  l.^(Lj+^)  -  I.l| ( 1  )  j 

=  -  M^'L^)  +  K, 

=  (Lj+M*) 

9  APPLICATION  TO  SOLAR  RADIATION  PRESSURE  (NO  SHADOtf) 

It  is  interesting  to  note  that,  of  the  results  for  the  complete  range  of 
values  of  n  (-co  to  -wo)  covered  by  the  general  theory  of  sections  4  to  6, 
interpretations  have  been  found  for  every  case  but  one;  for  n  =  2  to  oo  the 
results  apply  to  the  gravitational  field  outside  the  earth,  for  n  =  1,0  and  -1 


L3K2+M3Nl"3 


the  rosults  are  trivial  and  have  been  interpreted  in  eootion  6.1,  and  for 
n  s  -3  to  -co  the  results  apply  to  luni-solar  (gravitational)  perturbations. 

This  leaves  the  results  for  n  s  -2  unusod  and  it  so  happens  that  it  is  Just 
these  results  which  are  relevant  to  the  study  of  the  perturbations  due  to 
solar  radiation  pressure. 

It  is  neoessary  to  assume  that  the  satollito  remains  in  sunlight,  since 
in  shadow  the  radiation  pressure  disappears  abruptly.  Equivalently  wo  may 
think  of  tho  earth  as  transparent,  so  that  it  is  possible  to  have  a  conservative 
field  ovorywhero.  Then,  negleoting  any  asymmetry  in  the  shape  or  surface 
finish  of  the  satellite,  nrd  taking  the  sun  to  bo  at  an  infinite  distance,  the 
radiation  pressure  is  a  constant  force  -  directed  always  away  from  the  sun.  The 
direction  of  the  sun  is  an  axis  of  symmotry  and  the  potential  function  is  given 
by 


U  =  Fr  cos  9  ,  (l6l ) 

whore  9  is  tho  angle,  subtended  at  the  centre  of  the  earth,  between  the 

satellite  and  tho  sun;  F  is  the  (constant)  force  per  unit  mass  towards  the  sun  - 

2 

using  the  notation  of  Cook  -  and  is  always  negative. 

Thus 


U  = 


2 

on  taking  J_g  *  -  F  R  /p;  hero  J  g  is  positive  and  the  interpretation  of  the 
distance  II  is  unnecessary  since  it  will  be  cancelled  out  in  obtaining  tho 
perturbations  A£. 

Tho  treatment  of  section  8  is  now  immediately  applicable  (with  n  =  -2)  and 
there  are  again  two  cases,  according  to  whether  or  not  v/e  average  over  a  year. 

If  the  sun  is  assumod  fixed  wc  use  the  A,B,C  of  section  8.3*  Instead  of  working 
from  the  general  formulae  we  toko  advantage  of  the  formulae  for  ^  in 
section  6.4*  It  is  necessary  to  bear  in  mind  that  these  equations  must  be  taken 
a3  referring  to  the  dashed  elements  of  section  8,3.  To  obtain  normal  elements 
wo  use  equations  (74)  and  (79)  to  (82),  Since 


AS  = 


2x  F  a2 


9 


we  get,  very  easily,  the  following  results: 


Aa 

8 

0  , 

(162) 

Ae 

S 

^2J?a2(l-o2)^  (Boos  w- Asinw)  # 

(163) 

Ai 

S 

-  a2e(l-o2)^  Cooo  w  , 

(164) 

sin  i  AO 

a 

-  ^■a2e(l-o2)7  Csinw  , 

(165) 

Aw  +  oos  i  AO 

s 

•  ^a2e“\l-e2)^  (Acos  w+ Bsin  w) 

(166) 

and 

Aor  =  ^“a2(e"*^+e)  (Aooa  w  +  Bain  w)  •  ( 1 67) 

Tho  results  of  equations  (l63)  to  ( 1 66)  agree  with  those  given,  in  a 
different  form,  by  Allan  -  equations  (14)  of  Ref.  24. 

Passing  to  the  consideration  of  the  moan  effects  of  solar  radiation  over 

ono  year,  it  is  advisablo  to  interpret  R  as  tho  (variable)  distance  of  the  sun 

2  2 

from  the  earth.  If  we  then  assume  that  PR  =  const  =  P  a,,  where  a,  is  the 

0  d7  d 

semi-ma,jor  axis  of  the  sun's  orbit,  we  may  pursue  the  analysis  of  section  8. 4* 
Double-dashed  elements  are  based  on  the  ecliptic  as  reference  plane. 

2 

With  n  =  -2,  k  takes  the  value  1  only.  On  replacing  by  pQad  equation  (88) 

gives 


U-2(av)  .  -#P 


1 

2 


ae 


(sinw"  cos  i"  sin  (0^+v^)  +  cosw"  cos  (0^+v^)) 


ar &' 


(168) 


Hence 

U-2(av,  av) 


2% 

“  2  Po  ae(1"ed^”2  2x  J  lainu"  008  sin  (0^+Vj) 

0 

+  cos  «"  cos  (0^+v^)]  d  Vj 


=  0  . 


Thus 


(169) 


for  every  element. 
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10  SECULAR  AND  I.0NC-PKRI0D  PERTURBATIONS 
10.1  Introductory  remarks 

The  bulk  of  this  paper  has  been  devoted  to  the  first-order  perturbations 
of  a  satellite  orbit  over  a  period  of  one  revolution  of  the  satellite.  The 
perturbations  have  been  'first-order'  in  that  eaoh  term  of  the  perturbing 
potential  U,  given  by  equation  (ll),  is  treated  quite  separately.  The 
integration  of  Lagrange's  equations  (2)  to  (7)  has  been  carried  out  with  the 
elements  on  the  right  hand  sides  held  constant.  Perturbations  over  many 
revolutions  are  assumed  to  be  given  by  multiples  of  the  perturbations  over  one 
revolution. 


This  treatment  would  be  adequate  for  dealing  with  the  potential  field 

/  \  —10 

of  the  earth  if  it  were  the  case  that  each  J„  (n  >  2)  were  of  the  order  10 

n  -3 

say.  The  actual  position,  howover,  is  that  J0  is  0(10  )  and  the  other  J  , 

-6  ^  n 
with  n  >  2,  are  0(10"  ),  Thus  it  is  in  some  ways  more  logical  to  regard 

Jg  terms,  only,  as  first-order  and  the  terms  for  later  J^,  together  with 

terms,  as  second-order.  Up  to  the  present  we  have  preferred,  for  the  sake 

of  generality,  not  to  do  this,  but  now  it  is  desirable  to  consider  matters  in 

the  light  of  the  actual  geopotential  field. 


The  key  effects  of  the  dominant  Jg  term  euro  on  Q  and  w,  these  elements 
being  given  secular,  i.e.  non-periodic,  changes.  If  the  rates  of  secular 
motion  are  donoted  by  ft  and  u  respectively  -  the  dot  notation  is  used  nowhere 
else  in  this  Report  -  then 


ft  =  Jn 


'2  n(|)  °20  ani  "  =  J2  n  (f)  "20  * 

whore  n  is  the  mean  motion  and  w20  arG  2^VGn  ^y  section  6,3;  thus 

■  -  2  J2  n  (I)  003  1 


(170) 


and 


1  =  >J2  n(p)2(1-tf)  (17,: 1 

2  • 
where  f  =  sin  i.  The  clement  <r  also  has  a  secular  component,  c r,  but  this  has 

not  the  importance  of  ft  and  u. 
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In  section  10*2  we  oonsider  long-period  effects  arising  from  the  inter¬ 
action  between  u  and  the  (previously-called)  first-order  perturbations  associated 
with  Jy  eto;  perturbations  are  introduced  and  their  intoraotion  with  u 
considered  also.  In  seotion  10,3  we  consider,  similarly,  the  interaction  of  ft 
and  w  with  luni-solar  perturbations* 


These  introductory  remarks  are  concluded  by  a  brief  mention  of  short- 
period  perturbations.  Until  now  these  have  not  arisen  because  Lagrange’s 
planetary  equations  have  always  been  integrated  over  a  complete  revolution  of 
the  satellite.  Short-period  terms  are  those  which  arise  during  a  revolution 
but  vanish  at  its  end.  The  only  short-period  terms  which  are  significant  are 
those  associated  with  J A  good  way  to  derive  them  is  to  integrate  the  Lagrangian 
equations  using  true  anomaly  v  as  argument  instead  of  mean  anomaly  M, 


Merson^  has  applied  this  method  and  his  results  (for  the  standard  osculating 

elements)  are  given  for  integration  from  the  ascervling  node  of  the  satellite  orbit 

29 

as  initial  point.  Sterne  gives  results  which  apply  for  integration  from  any 
point  of  the  orbit,  the  node  being  one  obvious  choice  and  perigee  another,  ,/e 
follow  Sterne  and  list  here,  for  reference,  the  complete  J ^  perturbations  - 
integrated  to  a  general  point  of  the  orbit  -  for  each  of  our  3ix  standard  elements, 
using  v  as  argument.  The  elements  on  the  right  hand  sides  of  equations  (l 72 ) 
to  (177)  should  be  taken  to  refer  to  conditions  at  the  chosen  initial  point  of  the 
orbit  (when  t  =  0)  and  the  functions  P(e,i,w)  must  be  chosen  to  make  the 
expressions  vanish  at  this  point;  e.  g.  for  integration  from  the  node, 

PQ(e,i,w)  =  6w  +8  e  sinw.  The  expressions  agreo  with  those  given  by  Merson 
and  Sterne,  To  obtain  purely  short-period  terms,  it  is  only  necessary  -  for 
Q,  w  and  <r  -  to  replace  the  initial  secular  v  by  v-K, 

2  ' 

J0  R  a 

6a  =  — - r- 

16  p-5  2  2 

+  (4+ej  f  (cos  (2u+v)+cos  (2w+3v)j  +  12e  (2-3f)cos2v 

+  180  fcos2(«+2v)  +  3c',f  [cos  (2w-v)+cos  (2u+5v)] 

+  2e^  (2-3f)  cos  3  v  -  P  (e,i,w)]  ,  (172) 

a 


[l2  (2+  3e2)f  cos2  (w+v)  +  60  (4+e  ) (2-3 f)  cos  v 


80 


— — j  [6  (4+e^)  (2-3f )  cos  v  + 
32  P 


3(4+11e^)  f  cos  (2«+v) 

+  (28+17©2)  f  cos  (2u+3v) 


+  12e(2-3f )  cos  2  v  +  6ef  (lOcos  2  (w+v)  +  3cos2  (u+2v)| 

0  o 

+  3  e  f  [cos  (2w-v)  +  cos  (2u+5v)j  +  2e  (2-3f)cos3v  -  Fe(o,i,«)]  , 

....  (173) 


* 


Jg  R  sin  2  i 


[3  oos  2  (w+v)  +  3«oos  (2w+v)  +  ecos  (2w+3v) 


(174) 


J2  R  oos  i 

60  »  — — — —  [«6v  +  3oin2  (w+v)  -  6e  sinv  +  3esin  (2w+v) 

^ p  +  esin(2w+3v)  -  Ffi(o,i,w)]  , 


(175) 


j 

-2— r[24  (V-5f)v  +  6s’1  {4  (2-3f )  +  e2  (l4-17f))  sinv 
32 

-  3e-1(4f  +  e2(8-15f ))  sin  (2w+v)  +  e"1  (28f-e2(8-19f)l 

x  sin  (2w+3v) 

+  12  (2-3f)  sin2v  -  12(2-5f )  sin  2  (w+v)  +  I8f  sin2  («+2v) 

-  3ef  sin  (2w-v)  +  2e  (2-3f)sin3v  +  3ef  sin  (2«+5v) 


-  Pu(o,i,«)] 


(176) 


J  R2(l-e2)^  2 

— - 5 - [24  (2-5f)v  -  6e“1  (4-5®)  (2-3f )  sinv 


32  p‘ 


+  3e-1  (4+17e2)fsin(2w+v)  -  o”1  (28-1 3e2)f  sin  (2w+3v) 

-  12(2-3f)  sin 2  v 

+  36fsin2(w+v)  -  l8f  sin2  (w+2v)  +  sin  (2w-v) 

-  2e(2-3f )  sin  3v  -  3ef  sin  (2w+5v)  -  F^e,!^)]  . 

....  (177) 


It  is  to  be  noted  that  Sterne  does  not  give,  as  his  sixth  expression,  a 
formula  like  equation  (177).  Ho  gives,  instead,  a  formula  which  is  more 
convenient  for  the  derivation  of  M,  the  quantity  actually  required*  From 
equation  (8), 


6M  =  6cr  + 


■0 

j  nt  ; 


using  the  constant  n'  introduced  in  section  7  and  equations  (63)  and  (65 ),  with 
the  indicial  suffix  n  *  2,  we  then  have 


6M  a  n’t  +  6c  —  6p  • 


(178) 
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o 


1 

With  tho  expression  for  6p  available  (it  is  -  6<p  of  Merson  )  we  got 


J.R'O-e2)* 

a  nit  . - - - [&,  1  (2^-e^) (2-3f ) sin v  -  3e  n(4+5e  )fsin  (2w+v) 


6M 


32  p 


+  (28-e^) f  sin  (2w+3v)  +  12(2-3f)  sin2  v 

+  I8f  sin  2(aH-2v) 

-  3efsin(2«-v)  +  2e(2-3f ) sin  3v 

+  3 ef  ain  (2w+5v)  -  PM(e,i,w)]  , 


....  (179) 


For  the  pure  short-period  behaviour  the  n't  term  must,  of  course,  be  dropped. 


The  6p  term  in  equation  ( 1 78)  does  not  arise  if  Merson' s  smoothed 
elements  are  used  instead  of  osculating  elements.  The  short-period  expressions 
of  equations  (172)  to  ( 177)  then  become  very  much  simpler;  in  particular,  terms 
in  e  disappear,  though  we  have  already  -  in  section  8,3  -  made  the  point  that 
such  terms  aro  not  of  undue  significance. 


10,2  Secular  and  long-period  perturbations  associated  with  J^,  etc  and  with 


There  is  a  conceptually  straightforward  method,  associated  with  the  names 
of  Poisson  and  Poicare,  which  it  is  possible  to  apply  5.n  order  to  obtain 
perturbations  associated  with  successively  highor  powers  of  the  J  coefficients. 
The  method  has  been  applied  by  Merscn1  to  obtain  the  perturbations,  but 
already  at  the  expense  of  much  laborious  calculation.  For  higher  powers,  e.g, 

J 2  J y  it  becomes  quite  unmanageable,  but  it  is  nevertheless  worth  mentioning 
briefly. 


The  principle  of  the  method  is  an  iterative  one  based  on  the  assunption 
that  the  complete  perturbation  of  each  element,  £,  may  be  expressed  -  in  terms 
of  the  elements  aQ,  e^,  i  and  at  epoch  (t  =  0)  and  of  the  true  anomaly  v  - 
by  a  Taylor-type  series  of  the  form 


•  •  •  • 


S2,S3' 


(180) 


where  v  is  not  confined  to  values  within  one  revolution  of  t  =  0  (we  ignore 
convergence  questions,  of  which  there  is  a  short  discussion  at  the  end  of 
section  4.6  of  Ref,29)#  For  =  1  and  =  ....  =0,  the  functions  f 

are  precisely  those  implied  by  equations  (172)  to  (177).  For  =  1  and 
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S2  *  . .  tho  full  first-ondcr  perturbation  could  be  obtained;  and 

so  on,  Lagrange's  equations  for  these  first-order  perturbations  being  integrated 
v/ith  tho  elements  on  the  right  hand  sides  held  fixed  at  a  ,  otc.  Socond- 
order  perturbations  (with  £Sj  =2)  could  then  bo  obtained  by  substituting,  on 
tho  right  hand  sides  of  the  Lagrangian  equations,  functions  of  tho  elements 
v/ith  first-ordor  perturbations  includod.  Proceeding  in  this  way  perturbations 
of  any  ordor  aro  in  principle  determinate. 

Morson^  has  used  this  method  to  obtain  J?  terms  for  osculating  elements, 

*  30 

Although  ho  does  not  dorive  short-poriod  terms  (for  r/hioh  a  paper  of  Kozai 

2 

may  be  consulted)  the  method  is  already  extremely  laborious  at  the  J2  level, 

Tho  Mcrson  expressions  are  givon  for  the  variation  from  node  to  node;  from 
porigee  to  perigee  the  variations  are  different  (unlike  the  situation  for 
etc  variations)  due  to  tho  J2  variation  in  u.  The  perigee-to-perigee  expressions 
are  used  later  in  the  section  and  are  listed  below;  is  the  second-order 
change  associated  v/ith  the  fir3t-oidor  cliange  given  for  J2  by  equation  (}6)  with 
elements  in  equations  (38)  to  (45)  taken  at  perigee. 


Aaw  «  -  4  r4&2p“5  O+e)3  f(A-5f)  sin  2«  ,  (l8l) 

Aow  =  ’  ll  J2  (^f  (l+e)f  t^5f)  (5+7e)+o(l-o)(l4-15f)l  sin  2w  ,  (182) 

Aiu  s  sin2i  l(^15t)e2  -  4(Jp-5f)  (3+4®)  j  sin 2  u  ,  (183) 

Kh  3  “  if  A  (p)4  003  1  [^(9-10f)+e2(4+5f)-!8(3+4e)(2-5f)  -  2e2(7-15f)l 

oos  2u]  , 

••••  (184) 


Au  s 
u 


[2(l92-412f+215f2)  +  e2  (56-36f-45f2) 

+  [32c’1  f  (4-5f)  -  4(48-322f+3l5f2)  -  64e(4-26f+25f2) 
+  2e2(28-158f+135f*2)J  cos2«] 


(185) 


and 


t 
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A<r, 


w 


3§  ^  (j)  (l-e2)^r2(48-1l6f+67f2)+  I2e(2-3f)2-  2e2(l6-52f  +59f2) 

+  4e3(2-3f)2  -  e^35-41f  -  ^  f2^ 

-f  |l6o‘1(4^5f)-2(l70-2l3f)-4e(l46-l87f) 

+  8e2(49-57f)+4e3(l06-13lf)-e4(l06-111f)l  cos2w 

.*]  •  (186) 


-  ^  e^f2  cos  4 


Some  general,  but  lengthy,  expressions  for  JZ  perturbations  have  been 

31  ^ 

given  by  Claus  and  Lubowe  .  These  expressions  refer  to  a  complete  revolution 

of  the  satellite  from  v+w  =  0  to  0  +2x,  using  the  notation  of  Ref,  31.  For 

o  o 

0Q  =  0,  Merson's  expressions  are  confirmed.  For  0Q  =  uq,  however,  v  changes 
from  0  to  2x-Aw  and  so  we  do  not  get  our  equations  (l8l)  to  (186),  It  is 
interesting  to  observe  that  for  this  case  (0Q  =  wq)  Ref,  31  leads  to  a  J 2 
perturbation  of  semi-major  axis  equal  to  zero. 


The  general  method  outlined  above  is  impractical  then*  The  following 
point  about  it  should  be  noted,  however;  its  significance  will  appear  later. 
Just  as  (secular)  terms  in  v  appear  in  the  perturbations  associated  with 

Jo,Jt ,  ••••  and  there  similarly  appear  -  as  dominant  terms  over  a  long 

c  J  2  3 

interval  -  term3  in  v  associated  with  J y  J 2J^,  ••••••  and  J 2;  for  still 

higher  powers  of  the  J  coefficients  there  appear  higher  powers  of  v, 

A  suitable  treatment,  which  takes  account  of  the  main  higher-order  terms, 

2 

is  in  fact  rather  obvious  and  will  now  be  described*  Leaving  ^  terms  until 
the  end,  we  consider  the  perturbations  due  to  J^, J^,  etc;  we  work  generally 
uoing  J^,  where  4  is  appropriate  as  the  general  suffix  since  only  positive 
values  are  required  and  wo  want  n  to  denote  the  mean  motion.  Since  there  is 
no  long-period  or  secular  variation  of  semi-major  axis  we  start  by  considering 
eccentricity,  e. 


We  eliminate  short-period  terms  by  using  'barred’  elements  and  write,  from 
equations  (36)  and  (37), 


de 

dt 


4=2  k=0 


e 


Ik 


(187) 


Hence,  abbreviating  equation  (39)  to 


e4k  =  G4k^°,i^sin  f 
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whoro  tho  function  G^Oji)  is  thereby  defined, 

at  =  * 

e,k 

Now  our  interest  is  in  the  interaction  of  ( t  >  2)  with  J ^  Since 

a,  e  and  i  aro  constant  to  order  the  variation  in  e  is  given  by 

t  t 

6o  =  n^ (p)  G£k^e,i^  j  sink(w“2,0<it  * 
l,k  ° 

Taking  the  variation  of  u  to  be  given  by  equation  (l7l)»  we  may  write  (since 
only  terms  with  k  >  0  occur) 

6°  =  "  ~kJT  (^)  Gek(e,i)  (cosk  (w^-oosk  (w0-^)J  .  (l88) 

l,k 

Here  t  does  not  appear  so  that  6e  is  not  secular  but  long-periodic  (by  implied 
definition  of  'long-periodic’). 

Two  points  may  be  observed  about  equation  (l88).  They  apply  also  to  the 
perturbations  for  the  other  elements. 

First,  if  cosk  (u-£k)  ia  expanded  by  Taylor's  theorem  (valid  for  «  of 

any  magnitude)  with  w-w  =  —  (v-v  ),  short-period  terms  being  neglected,  then 

o  n  o 

-  -V  [cosk  (u--^)-cosk  («Q-^)|  =  (v**vo)  sin  k(wQ«*^) 

+  i  ~  (v"v0)2  °osk  (wQ-^)  . . 

With  the  first  term  here  we  have  recovered  tho  appropriate  contribution  to  65 
in  equation  (l80)  with  s  1  and  the  other  =  0.  With  the  second  term  we 
havo  precisely  the  dominant  (quadratic)  term  of  those  associated  with 
of  which  tho  significance  was  forecast.  Similarly,  successive  terms  are 
associated  with  dominant  terms  in  v^,v^  etc  associated  with  J^,  eto. 

The  second  point  to  be  observed  about  equation  (l88)  is  that  the  lower 

limit  uq,  corresponding  to  t  =  0,  appears  in  it;  The  term  in  «Q  is  often 

omitted  for  simplicity,  for  example  in  the  work  of  Kozai^.  The  variation  in 

e  is  then  given  from  its  value  at  sane  time  for  which  w  is  such  that 

cosk  (w—'g’x)  =  0.  The  present  author  feels  that  one  may  then  be  only  too  easily 

•  2 

confused  by  what  happens  when  u  is  small.  This  occurs  when  sin  i  is  close  to 
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4/5;  i,©#  when  i  is  about  63°, 4  or  11 6°, 6,  the  ao-called  critical  inclinations. 
If  the  term  is  dropped  from  equation  ( 1 88 ) ,  then  when  «  is  small  the  value 

of  6©^  i©  in  general  going  to  be  large  (rc1  itive  to  J^)  and  there  appears  to 
be  a  singularity,  The  trouble  is,  of  course,  that  when  w  is  small  the  perigee 
is  moving  so  slowly  that  the  time  at  which  cosk  (w-'^c)  =  0,  if  it  exists  at 
all,  may  be  a  long  way  from  the  time  at  which  the  perturbation  is  of  interest. 
Many  papers  have  been  devoted  to  the  subject  of  the  critical  inclinations, 
e,g.  Refs, 32,  33  and  34>and  some  rather  advanced  mathematics  is  involved.  It 
is  perhaps  sometimes  overlooked,  however,  that  the  situation  is  in  some  ways 
simpler  than  for  non-critical  inclinations.  For  if  w  is  zero,  w  is  no  longer 
a  different  type  of  element  from  e  and  i  in  its  Jg  behaviour.  The  integration 
of  equation  (l87)  then  leads  simply  to 


6,k 


to  which  equation  (188)  reduces  as  w  -*■  0. 

The  long-period  perturbation  to  the  element  i,  corresponding  to 
equation  (188),  is  given  by 

=  e(1“e2)*"1  coti  G«ek^e,i^  fcoslc  ("HW -cosk  («0-4*)!.(l89) 

l,  k 


For  the  elements  Q,  w  and  o’  there  is  a  complication  arising  from  their  main 
secular  perturbations.  For  n,  we  start  by  writing  (cf,  equation  (41 )) 

ft  =  nX/4(p)  H«k(s,i)  oosk  » 

C,  k 

where  ,i)  cosk(«-,§^)  =  n^k. 

For  each  C  >  2  and  k  =  0,  we  get  a  secular  term  contributing 

nJe(p)  H«o(e»i)t  (190) 

to  6 0,  For  k  >  0,  we  get  the  long-period  contribution,  corresponding  to 
equation  (l88),  given  by 

1We,i^  fsink  -sink  (og-in)l 


(191) 
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but  it  will  now  be  seen  that  this  ter*«  may  be  oombined  with  part  of  the  main 
secular  (ft)  term  for  l  e  2  and  k  =  0.  The  latter  is  given  by  equation  (170), 
which  it  is  neoessary  to  write  more  precisely,  with  n'  in  place  of  n.  We  have 


n' J2  (p )  «20<M)  -  A  = 


-  2rfJ2  (Fl  d-®2r2oosi  , 


so  that 


Jj  =  -  6  n!  Jg  (“)  e(l-e2)_100Bi 


IT  •  2n'J2^)  8lni  • 

Hence,  using  equation  (40)  and  writing  n'  only  where  necessary, 

2  2 

ft  =  -  |n'J2  (^)  008  ^  "  ¥n  J2  (|)  eO-o2)-1  oo»  i  65 


and  so 


t  2  2 

J  ft  dt  =  -  2"n,J2  oosiQt  -  ^nJ2  ©(l-©2)"1  cos  i  ^  8e  dt 


But,  from  equation  (l88). 


|  Seat  =  S4k(e,i)^V  (sink  (u-^)  -  sink  («o-|«)| 

3  «,k  .  ,_n 

-tcosk  (wq-^) 


Hence  the  combination  of  the  Jg  secular  term  and  the  terms  of  type  given  by 
equation  (191)  -  those  with  k  >  0  -  may  be  represented  by  a  secular  term 

-  fn'J2  oos  iQ  +  SeCl-e2)"1  G^k(e,i)oos  k  («,-£*)  jt  (192) 


together  with  long-period  terms 


i 
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z^d)  [H^(e<i)+2  aIs^SrA  ^(e>i)] ieink  (i^)-ink  • 

«,k  \ 4  / 

....  (193) 

The  oomplete  perturbation  6 is  made  up  from  the  contributions  (190), 

(192)  and  (l93)«  The  presence  of  the  second  term  in  the  square  brackets  of 

(193)  explains,  for  the  case  1*5  and  k  =  1,  an  apparent  discrepancy  between 

1  12 

formulae  of  Merson  aid  Kozai  ;  Refs,  1  and  12  differ  in  having  faotors 
(l  -l5f/4)  and  (l-5f/4)  respectively  but  the  discrepancy,  noted  by  A,  H,  Cook^, 
is  now  seen  not  to  be  a  real  one.  Pull  agreement  between  Refs.l  and  12  has  been 
demonstrated  in  a  further  paper  by  Merson  . 

It  is  observed  that  one  unfortunate  result  follows  from  the  separation  of 
the  Jg  perturbation  into  secular  and  long-period  c  exponents :  near  the  critical 
inclinations  each  component  becomes  very  large,  though  the  combination  remains 
finite  (and  continuous)  always.  Thus  if  u  is  near  to  zero  the  separation  should 
not  be  made. 

The  perturbation,  6 w,  of  u  may  bo  dealt  with  in  the  same  way.  If  I^k  is 
such  that 


Iek(o,i)°°sk  (u-s«)  =  , 


then  6(3  has 


(194) 


arising  for  oaoh  C  >  2  and  k  =  0,  Terms  vdth  k  >  0  may  be  combined  with  the  main 
Jg  secular  term  and  represented  by  a  secular  term 


3J2“'(0  f(l-K)+  ioO-e2)"1  (13-150  Y$(?f  ®4k(e'i) 


C.k 


(195) 


together  with  long-period  terms 


Z  S6  d/  [I*(,'i)  *  G«k(°> i}] fsink  - sinic  • 

e,k  ....  (196) 
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Similarly,  if  K^k(e,i)  is  such  that 


K^,k(o,i)  oosk  (w-W  s  c T/v  , 


ek 


thon  the  perturbation  6<r  has 


n  Je  (ft  IWe’i^t 


(197) 


arising  for  each  C  >  2  aivi  k  a  0,  The  combination  of  terms  for  which  k  >  0 
with  the  main  Jg  seoular  term  may  be  represented  by  a  secular  terra 


1 

2 


^(0  (I-#*  [('-k)*  2e(i-e2)'1 

i,k 

(198) 


•  ••• 


(with  Jg  instead  of  u  in  the  denominator  of  the  summed  terms  since  a  factor 
(l  -  5f/4)  has  been  cancelled),  together  with  long-period  terms 

EnJ//D\^r 

k5 

e,k 


S^(p)  K^k^e,1^~  —  GCk(®,i)  ^sink  (u-^c)-’Bink  (W0-4^)J 


•  ••• 


(199) 


But,  as  with  the  short-period  terms  listed  in  seotion  10.1,  it  is  the 
perturbation  in  M,  rather  than  that  in  or ,  which  is  important.  To  obtain  this 
wc  follow  the  method  of  section  7,  where  n'  was  introduced  by  equation  (63). 
Thus  it  follows  from  equation  (8)  that 

T 

§  .  n'*f  /  (n-h')dt+ff  ,  (200) 

(0  J 
o 

where 


using  equations  (63)  and  (68),  If  L^k(e,i)  is  defined  suoh  that 


Lek(e,i)  cosk  («-^)  =  - 


3  plk 
2(£+1 ) 
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than,  since  K^e,!)  +  Lg^e,!)  &  0,  the  complication  of  terms  in  G^(e,i)  does 
not  arise*  The  pertuxbation  55  consists  of  the  secular  term 

n’  Jt  (|j  !k4o(o>1)  +  L<o(e'i)!t  (201) 

for  eaoh  l  >  0  and  k  s  0,  the  main  secular  term 

n't  ,  (202) 

where  n'  is  absolutely  constant,  and  long-period  terms 

^"kjf  (^)  +  Lek^e,i^  fsin  k  (““M  "  8in  k  (-o-^)I  •  (2°3) 

l,k 


We  now  outline  the  similar  analysis  which  may  be  carried  out  to  obtain 
secular  and  long-period  perturbations  associated  with  J^,  Here,  as  has  already 

been  mentioned,  the  expressions  obtained  depend  on  the  point  of  the  orbit  with 

which  the  orbital  elements  are  associated.  If  the  smoothed  elements  of 
1  12 

Me r son  or  Kozai  are  used  the  long-period  perturbation  for  semi-major  axis 
turns  out  to  be  zero.  For  osculating  elements  defined  at  the  node  or  perigee 
this  is  not  the  case.  Here  we  give  perturbations  for  elements  at  perigee,  for 
which  the  second-order  changes  over  one  revolution  have  been  given  by 
equations  (l8l)  to  (186), 


From  equations  (l8l)  to  (l83)  the  long-period  perturbations  of  a,  e  and  i 
are  given  by 


5  a 


w 


6  o 


u 


l  n'  4  a2 

”*5~  P5 


(l+e)^f  (l (oos  2  u- oos2  w0)  , 


(204.) 


3  1  j2 

2  (|)  0+«)f  ^1-Jf^(5+7e)+e(l-e)  }£fJJ  (oos 2  w-  oos 2«q) 

....  (205) 


and 


6  i 


w 


8i'2  (|)  ain2i  [(''t4')  (3+4c)‘(8'lff)°2]  (<=oa2u-ooa2uo)  , 


#  •  •  • 


(206) 


c 
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since  (with  ffg  -  £p2  ■  0 

f 
~ta 

the  distinction  becomes  vital  for  the  perturbations  60u>  an^  gl’v®n  below* 
Zero  suffixos  refer  to  initial  perigee  values* 

For  60  the  procedure  is  similar  to  that  for  the  J«  perturbations* 

w  .2  ^ 

Combining  the  terms  with  the  main  Jg  secular  term  it  may  be  shown  that  the 

total  perturbation  is  representable  by  a  secular  term 


respectively.  In  the  above  equations  n’  is  usod 

from  section  7)  o*  ,  o  » 

fU  n'  ♦  0(J?,J  ....) 


-  |  »’  J2  0  °osi0[l  <■  i  J2  [(9-1  Of)  ♦  ie2  (4+5f)  +  5f  «2) 


Tl 

0)  1 
o 


x  cos  2  wit. 


(207) 


together  with  a  long-period  term 


3'l6  ■  2  (|)4 oos  i  £(3+4*)  (4-5f )  -  (sin2  u- sin 2  uq)  .  (206) 


Similarly,  the  combination  of  with  the  main  Jg  secular  term  may  be 
represented  by  a  secular  term 

2  f*/  ,2 

3  J2  w  [0 ' 4f °)  +  128  J2  (I)  [z(l92^2f^5f2)  +  *2(5<S-36f^5f2) 

+  I6f(l3-15f)  (3^- 4^5^  »2)  0082  “oj]  * 

.*..  (209) 


together  with  a  long-period  term 


128 


Tp  (|^4  Tl  6e-1  vit-5f)  (f-2«2+2e2f)- 2(48-1 66f+135f2) 

,  e20l2-74^|gOf2-1I5f3n  (8ln2M  ■  aln2  uQ)  . 


(210) 


For  the  final  element,  o',  the  coiribination  of  60^  with  the  main  secular 
term  may  be  represented  by  a  secular  term 


76 


2 

2 


X  ^2 (48-1  l6f+67f2)+  12e(2-3f)2-2e2(l6-52f+59f2) 

+  4e3(2-3f)2-  o^35-41f-^  f2^ 

+  3f(l-0  ^4(4- 5f)  (3+4©)  -  ©2(l4~15f  cos  ^ojj* 


together  with  a  long-period  term 


••••  (21l) 


-  — — f  (i-e2)"^ni6e“1(4-3f)-  2(98-123f)- 4© (98-1 27f)+ e2 (206-231  f) 

128w  W  L 

+  4«3(58-7lf)-2e4(32-33f)l  (  sin2  w -  sin2  uQ) 

+  ^- e^f  (sin4  w  -  sin4  wq)J  .  (212) 

Again,  however,  it  is  the  perturbation  in  M,  rather  than  that  in  o’,  which  is 
important.  This  may  be  obtained  by  starting  from  equation  (200)  where  now, 
introducing  the  second-order  term  into  equation  (63), 

T  T 

t  W  Z  O  IJ2 

J  (n-n’)dt  =  -■jAp-gj  na  —  dt  . 


Since  ^p^  =  OV,  it  follows  that  we  may  write 


dt  =  n’  +  n'  ^  (p)  (°22~  2  P22"  v22^ 

where  o’ ^ »  P22  v22  0X9  defined  as  follows:  0^  i©  suc^  that 


(213) 


Ao’ s  2k  0’22  for  Ac  given  by  equation  (186);  i©  given  similarly  by 

Ap  where,  extending  the  set  of  equations  (l8l)  to  (186), 

Ap  =  ^  *4  (p)4  (l-e2)"^r2(48-108f+51f2)  + 12e(2-3f)2- e2(72-192f+l89f2) 

+  4©3(2-3f)2- 3e^10-1Af  + J  f2^ 

+  f|l2(26-33f)+  12e(38-49f)-  6e2(62-75f) 

-  12e3(30-37f )  +  72e4(l-f  )J  oos2  «  - 1  eVcos  4«j  ; 

(214) 


•  •  •  • 
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Vgg  is  similar,  with 


Av  s 


_  1 
~  2 


T 

/ 


2 

na  dt 


and  it  follows  from  Merson's  work  that 

+  jig"  ^  008  4  wj  • 

Hence,  after  some  simplification, 


e  cos2w 


(215) 


f§  =  n'  4  (|)4  C1-®2)3^  !(8-8f-5f2)+2f(8e“1(4-5f)-(l4-15f))co5  2w)J, 


Thus  we  obtain,  for  M  ,  a  secular  term 


n'  ['-^.4  (|)4  (1-«2)3/2  (8-8f-5f2)Jt 


(216) 


and  a  long-period  term 


3  n’  4 

f  (l-e2)3/2  |8e"1(4-5f)-  (14-15^)1  (sin2  w- sin2  u>  ) .  (217) 


64  w  \P/ 


It  is  now  desirable  to  summarise  the  results  of  this  sub-section,  showing 
how  the  formulae  lead  to  the  possibility  of  an  algorithm  for  computing  the 
position  in  space,  at  a  given  time,  of  an  earth  satellite  which  is  subject  only 
to  earth  asphericity  perturbations.  We  suppose  that  the  position  and  velocity 
of  the  satellite  at  a  given  time,  t^,  are  known  and  that  its  position  at  time 
tg  is  required,  A  possible  sequence  of  steps,  in  outline  only,  would  be  as 
follows : 

(i)  obtain  the  osculating  elements,  a^,  e^,  i^,  Q^,  and  M,j  (a  0^  a  x^)  at 
the  epoch  t.j ,  by  standard  formulae2^;  also  the  true  anomaly,  v^ ; 

(ii)  from  equation  (10)  obtain  U,  given  r,  and  0.  at  t,,  and  hence  obtain  a’ 

2  3  1  1  1 

from  equation  (62)  and  n'  from  n1  a ’  a  {i; 
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(iii)  obtain  the  initial  purely  ahort-period  terms,  relative  to  elements  at  perigee, 

from  equations  ( 1 72 )  to  (178)  and  ( 1 79)  with  secular  terms  omitted  (v  replaced  by 
v-M  in  6n  and  8w,  and  n’t  dropped  from  6M)  and  the  appropriate  F(e,i,w)  terms 
substituted;  remove  these  short-period  terms,  6^a  etc,  from  a^,  «««.,  to  give 
a0  (=  -  81  a)  etc; 

(iv)  obtain  u  s  3  ^  n*  a^s0  *'unc**ona  ^o*  eo  ar*^  *o 

which  are  required  for  secular  and  long-period  perturbations; 

(v)  obtain  «*,  where  «•  =  +  w  (tg-t,,); 

(vi)  obtain  the  combined  secular  and  long-period  perturbations  6gS,  ....  6gM  ; 
for  example  is  given  by  contributions  from  equations  (190),  (192),  (193), 

(207)  and  (20 8),  with  p,e,i,wQ,w  and  t  replaced  by  pQ,  eQ,  iQ,  w^,  w’  and 
tg-t^  respectively;  the  duplicated  main  term  in  equations  (192)  and  (207)  is# 
of  course,  taken  once  only; 

(vii)  obtain  a^,  •  •«.,  from  a^  =  aQ  +  6^a  etc;  hence  also  v^  (using  Kepler’s 
equation) ; 

(viii)as  in  step  (iii)  obtain-short-period  perturbations  6^a  etc,  based  on  a^  etc; 
hence  obtain  a2,  from  a2  s  a^  +  6^a  etc; 

(ix)  obtain  position  from  the  osculating  elements  a2,  e2,  i2,  02,  w2  and 

Steps  (i)  to  (iv)  are  iidependent  of  t2  and  would  be  performed  once  only  if 

position  were  required  at  several  times  No  difficulties  occur  near  the  two 

critical  inclinations  if  the  formulae  are  handled  correctly,  but  the  usual 
12  21 

problems  arise  *  for  the  singularities  in  the  elements  at  e  =  0,  i  =  0  and 
i  =  7C. 

10,3  Secular  perturbations  associated  with  a  disturbing  body 

It  is  supposed  that  the  orbit  of  the  satellite  is  close  enough  to  the  earth 
for  the  secular  variations  in  0  and  w  due  to  the  earth’s  oblateness  to  be  still 
the  dominant  effects.  Rates  of  change  of  satellite  elements  due  to  the 
principal  terms  of  the  disturbing  potential  of  a  distant  body  have  been  given 
by  equations  (l36)  to  ( 1 60 ) ,  averaged  with  respect  to  the  mean  anomalies  of  both 
the  satellite  and  the  disturbing  body. 

To  obtain  perturbations  over  a  period  which  is  long  in  comparison  with 
the  periods  of  the  satellite  and  the  disturbing  body  it  is  necessary,  as  in 
section  9*3,  to  consider  the  secular  variation  of  w.  Due  to  tho  presence  of  0 
in  tho  equations  it  is  now  also  necessary  to  consider  tho  secular  variation  of 
this  element.  Following  the  example  of  Smith  ,  we  do  not  give  long-period 
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variations,  which  are  liable  to  become  rather  complicated,  but  only  the  secular 
variations  obtained  by  averaging  equations  (136)  to  (160)  with  respeot  to 
0  and  «, 

Sinoe  fl  only  occurs  as  part  of  the  trigonous trioal  argument  (Q-flj)  it  is 
actually  with  respect  to  this  quantity  that  the  first  averaging  takes  place. 

To  avoid,  for  n  =  -5  and  k  =  0,  terms  containing  factors  sin 2  and  oos2«d 

we  also,  in  effect,  average  with  respect  to  w, ;  these  terras,  containing  also 

2  a 
the  factor  e^,  are  in  any  case  negligible. 

Since  averaging  with  respect  to  u  is  being  performed  secular  perturbations 
associated  with  any  k  *  0  must  vanish.  Of  the  (n,k)  values  considered  in 
section  8  this  leaves  us  with  (-3,0  and  (-5,0)  only.  For  the  first  case, 
equations  (l 36)  to  (140)  lead  easily  to 


dt 


dt 


0 


(218) 


dt 


-  2  8inZid) 


(l-e2r(l*fe2) 

cosi  , 


(219) 


d  u  ,  d  0  _  . 

+  cos  i  - 


dt 


d  t 


b  ad  n 


0- 


-e2)7  si*2id^  (l  ”  2  sin2i) 


(220) 


and 


d  ^.O  _ 
dt 


4  a; 


0-4r}/2  (7+3e2)  (l  -  f  sin2id) 
(l-^sin2i)  .  (221) 


26 

The  secular  results  given  by  Kozai  follow  at  once  from  equations  (219)  and 
(220),  to  the  order  in  e  to  which  he  works.  We  recall  that  the  inclination  i^ 
is  relative  to  the  earth’s  equator  and  is  not  the  (almost  constant)  inclination 
relative  to  the  ecliptic. 
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For  the  other  case,  equations  (156)  to  (160)  lead  to 


d  e  d  i 

_=5iQ  -  1  -5.0  .  o 

dt  “  dt 


(222) 


d  ”-5.0 
dt 


-45li, 


1024  a^ 


-  a2(l-4)-7/2  . 1  .J)  (l-e2)-i  (l  +  5e2  +  £  .*) 


(8-40  sin2id  +  35  sinS.^)  (4-7  sin2i)  cos  i  , 


(223) 


d  u-5  0  d  ^-5  0 

_iO  +  ooai__iO 


x  (8-40  sin2i.  +  35  sinh, )  (8-40  ain2i  +  35  BirAi) 


(224) 


and 


d  g-s.o  .  -8^d 


dt 


x  (8-40  siiAi,+  35  sirAi.)  (8-40  sini+  35  aiiAi) 


(225) 


pc 

Smith  gives  expressions  containing  terms  covered  by  our  equations  (222) 
to  (224).  These  terms  in  Ref. 25  contain  a  number  of  errors,  however,  as  has  been 
confirmed  by  Smith  in  the  private  communication  mentioned  in  section  8. 

It  should  be  noted  just  why  the  assumption  had  to  be  made  at  the  beginning 
of  this  section  that  the  perturbations  due  to  oblateness  we  re  still  dominant. 

The  interaction  of  these  perturbations  on  those  arising  from  the  disturbing  body 
has  been  considered,  but  the  interaction  of  the  latter  on  the  former  has  been 
disregarded.  Thus  to  analyse  the  complete  effect  of  the  mutual  interaction  on  i 
and  Q  it  is  necessary,  for  first-order  results,  to  integrate  the  simultaneous 
differential  equations  for  di/dt  and  dfl/dt  given  by  equations  (137)  and  (138) 
with  the  oblatcness  terms  added  on  the  right  hand  sides.  Analysis  on  these 
lines  has  been  carried  out  by  Allan  and  Cook^, 

The  above  remarks  make  the  higher-order  terms  given  by  equations  (222)  to 
(225)  rather  irrelevant  in  the  case  of  luni-solar  perturbations  of  an  earth 


satellite;  for  if  a/a^  were  small  they  would  be  very  small  and  otherwise  the 
treatment  would  be  invalid  anyway.  These  terms  would  only  have  been  useful  if 
the  earth  had  been  much  more  oblate  or  if  the  moon  had  been  lighter  and  in  a 
closer  orbit, 

A  final  remark  which  is  worth  making  in  connection  with  perturbations  of  a 
satellite  due  to  a  disturbing  body  concerns  the  perturbation  of  the  mean  anomaly. 
As  in  section  10.2  it  is  the  perturbation  of  M  rather  than  that  in  o'  which  is 
really  required.  Returning  to  the  case  of  the  disturbing  body  regarded  as  fixed 
during  ono  revolution  of  the  satellite  (section  8.3),  the  mean  rate  of  change 
of  M  is  given,  to  first  order,  by 


where  both  Acr  and  Ap  are  given  by  equation  (94)  of  which  the  left  hand  side 
actually  is  Ap.  However,  the  situation  is  more  complicated  when  allowance  is 
made  for  movement  of  the  disturbing  body,  since  the  equation  of  energy  no 
longor  leads,  as  in  section  7,  to  a  constant  n*. 

11  CONCLUSIONS 

The  disturbing  function  for  an  axi- symmetric  potential  field  has  been 
written  as  the  doubly  infinite  Legendre  series 

oo 

-  V  M  Jn  rV"-1  P^sin  P)  , 

n=-oo 

where  r  and  ^x-0  are  polar  co-ordinates,  6  =  n  if  n  >  0,  and  &  =  -n-1  if  n  <  0, 
For  the  field  of  the  earth,  assumed  axi-symmetric,  the  Jn,  for  n  >  2,  are  the 
now  standard  coefficients,  R  being  the  mean  equatorial  radius  of  the  earth; 

J-|,  J0  an*  simple  meanings  and  may  easily  be  made  to  vanish;  the  series 

with  the  for  n  <  -2  only  may  be  applied  to  the  study  of  luni-solar 
perturbations  of  an  earth  satellite,  if  the  and  R  are  properly  interpreted. 

The  effect  of  the  general  term,  b’n,  of  the  series  on  the  motion  of  a 
satellite  in  the  f?eld  has  been  obtained  by  first  expressing  Un  in  terms  of  the 
usual  orbital  elements  and  mean  anomaly,  and  then  averaging  with  respeot  to 
the  mean  anomaly  to  eliminate  short-period  terms.  The  average  becomes  the 
sum  of  a  series  of  terms  in  cosk  (w-^rc),  with  0  <  k  <  l.  For  the  general  such 
term,  U^av),  the  changes  in  the  orbital  elements  over  a  complete  revolution 
of  the  satellite  have  been  found  by  making  use  of  Lagrango’s  planetary  equations. 
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Formulas  for  the  secular  and  long-period  motion  of  a  satellite  in  the 
gravitational  field  of  the  earth  have  been  developed.  The  long-period 
perturbations  arise  from  the  secular  motion  of  perigee  associated  with 
At  the  two  oritioal  inclinations  this  motion  oeases  to  exist,  but  the  formulae 
for  secular  and  long-period  motion  of  each  element  have  been  developed  in  suoh 
a  way  that  for  the  oombined  perturbation  there  is  no  singularity  near  the 
oritioal  inclinations. 

In  all  sets  of  formulae  for  the  perturbation  of  orbital  elements  the 
perturbation  of  the  sixth  element,  omitted  by  many  authors,  has  been  included. 
By  use  of  the  constant  mean  motion,  n',  of  Merson  the  formula  for  the  long- 
period  perturbation  of  mean  anomaly  has  the  same  accuracy  as  the  formulae  for 
the  other  five  orbital  elements.  It  should  bo  remarked,  however,  that  for  an 
actual  close  earth  satellite  the  long-period  variation  of  mean  anomaly  is 
essentially  governed  -  in  a  largely  unpredictable  fashion  -  by  atmosphoric  drag 
when,  in  addition,  n'  is  no  longer  constant. 

Application  to  the  study  of  luni-solar  perturbations  nas  required 
consideration  of  more  than  one  axis  of  symmetry.  Formulae  for  the  perturbations 
of  elements  defined  iu  terms  of  a  pseudo-equator  have  accordingly  been 
transformed  into  formulae  in  the  normal  elements.  The  plane  perpendicular  to 
the  direction  of  a  disturbing  body  is  one  such  pseudo-equator  and  the  orbital 
plane  of  the  body  another. 


Appendix  A 

LEGENDRE  EXPANSIONS  FOR  GRAVITATIONAL  POTENTIAL 


Consider  first  the  simple  case  of  Pig,  1(a)  in  which  P’  represents  the 
location  of  an  isolated  point  mass  m',  P  is  a  general  point  and  0  is  an 
arbitrary  origin  of  polar  co-oidinates. 

The  potential  at  P  due  to  P*  is  m'/fc'P.  Taking  the  axis  of  polar 
co-ordinates  along  OP'  we  have 


m* 

P'P 


m' 


/  2 

(r  +  r 


,2. 


2zr'cos  0) 


I 


(A1) 


There  are  two  cases  to  consider  according  to  whether  r  >  r'  or  r  <  r*# 

If  r  >  r',  the  expansion  of  m'r”^  |l-2(rf/r) cos  0  +  (r'/r)^J“^  as  an 
infinite  series  in  powers  of  r'/r  is  valid.  As  is  remarked  at  the  beginning 
of  Appendix  B,  this  expansion  may  be  used  to  define  the  Legendre  polynomials 
P^(cos  0)  and  in  fact 

°°i  C 

-  m'  }  -^e+T  pe  (°080)  •  (A2) 


pT 


"E 

U  0 


If  r  <  r',  we  have  similarly 

pV  =  (cos0) 


(A3) 


1=0 


Next  consider  an  extension  to  the  case  of  several  point  masses;  in 
Fig,l(b)  there  are  three,  0  is  still  an  arbitrary  origin  and  the  polar  axis 
OZ  is  now  also  chosen  arbitrarily;  OX  is  an  arbitrary  line  of  zero  azimuth  in 
the  plane  perpendicular  to  OZ, 

The  potential  at  P  due  to  m'  at  P'  may  be  expressed  as  a  series  in  the 
P^  (cos<p).  In  order  to  add  this  potential  to  that  due  to  the  other  masses  it 
is  desirable  to  rewrite  the  series,  expressing  it  in  terns  of  the  angles  0  and 
0',  V/e  note  that,  if  X  and  X’  are  the  azimuths  of  OP  and  OP'  respectively, 


cos <p  =  cos0'  cos0  +  sin0'  sin0  cos  (X-X1) 


(M-) 


Hence  the  addition  theorem  proved  in  Appendix  C  may  be  applied  and,  from 
equation  (Cl),  using  the  associated  Legendre  functions  defined  in  Appendix  B, 
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Appondix  A 


P^OOSf) 


c 

y'Um{e^}f  ^(eo» «’)  1^(005  8  )  008 m(X-X')  ,  (A5) 


msO 


where  =  1  if  m  =  0  but  2  if  m  >  0, 

For  tho  case  r  >  r*  we  now  have,  from  o qua t ions  (A2)  and  (A5), 

00  t 


potential  (due  to  m’) 


II 


)  r 

L _ i 

£>aQ  ITfcrO 


-£-1 


P?(co8  0)  (C^(m’)  cos  mX+S^m')  sin  m\) 


•  •••  (A6) 


where 


C%n*),  S%n»)  =  m»  u  r'C  (<^)i  P^(C03  6')  (cos  mX',  sin  mX'). 

The  double  series  is  convergent  aral  must  remain  so  on  combining  with  equivalent 
series  for  other  point  masses,  m”,  m'"  etc,  so  long  as  r  >  r",  r  >  r'”  etc. 

The  potential  duo  to  a  finite  set  of  masses  is  simply  that  given  by  equation  (A6) 
with  C^(m’)  and  S^(m')  replaced  by  C^(m')  +  C^(m")  +  ....  and  S^(m')  +  S®(m”)+  ... 
respectively. 

When  we  consider  the  potential  at  P  due  to  an  actual  gravitating  body, 
such  as  the  earth,  there  is  no  difficulty  with  the  convergence  on  passing  from 
a  sum  to  an  integral,  so  long  as  r1  <  r  for  all  points  P'  of  the  body  -  in  any 
case  it  is  arguable  whether  a  finite  body  contains  a  finite  or  infinite  nunber 
of  points.  Thus  for  any  given  distribution  of  matter  the  general  expression 
for  the  potential  at  any  point  outside  &  sphere  containing  all  the  matter  is 
given  by 

oo 

U.  .  £  P^(oos6)  ((£  cos  mX  +  S*?  sin  mX)  ,  (A7) 


t= 0  msO 


where 


C™,  Sf  =  J  um  t'C  ^(oos0')  (cos  mX*,  sin  mX^dm* 

(matter) 

The  expression,  by  equation  (9)  of  section  4.1,  for  the  earth’s 
gravitational  potential  now  follows  with  minor  changes  of  notation. 
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By  consideration  of  the  case  r  <  r'  we  conclude,  similarly,  that  if  a 
given  distribution  of  matter  occurs  entirely  outside  a  certain  sphere,  then 
the  potential  at  all  points  inside  this  sphere  is  expressible  by 

oo  C 

U*  s  ^  ^  rv  P^(cos  0)  ($£  cos  mX  +  3^  sin  mX)  ,  (A8) 

1=0  m=0 

where 

S™,  ff  3  J  u^  F^(cos  0')  (cos  mX’,  sin  mX')  dm'  • 

(matter) 


The  results  expressed  by  equations  (A7)  and  (A8)  may  be  combined  into  a 
single  equation,  the  validity  of  which  holds  throughout  any  empty  region 
bounded  by  two  spheres.  Ue  take  n  to  range  over  all  integers,  positive  and 
negative,  defining  C  =  n  if  n  >  0  and  C  =  -n-1  if  n  <  0,  The  equation  is 


U» 


oo  t 


n=-oo  m=0 


P^(cos  0)  (Cf  cos  mX  +  sin  mX) 


(A9) 


where 


Sj  if  n  <  0  . 


Absence  of  matter  from  the  spherically  bounded  region  is  a  sufficient 
condition  for  the  convergence  of  the  series  but  not,  in  general,  a  necessary 
one.  For  evidence  that  it  may  not  be  necessary  section  24*06  of  Ref, 37  should 
bo  consulted.  For  a  case  when  it  ^is.  necessary,  see  Fig.  l(c).  Matter  with 
unit  lino  density  is  assumed  to  be  distributed  between  A  and  B,  points  for 
both  of  which  0=0  and  for  which  r  s  a  and  b  respectively.  For  r  >  b,  the 
potential  is 


tfr1  -  a**1 

(1*1)  rM 


p,(cos  e) 

V 


9 


which  at  a  point  P  for  which  0  =  %  gives 
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For  r  <  a  there  is  a  similar  expression  in  positive  powers  of  r,  giving  the  same 
sum.  For  a  <  r  <  b  the  potential,  still  along  0  =  x,  is  of  course  given  by  the 
same  logarithmic  function,  but  now  it  has  no  expansion  in  powers  of  r. 
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LEGENDRE  FUNCTIONS 

The  basic  functions  are  the  Legendre  polynomials  Pn(x)  which,  as  usual 
in  mathematics,  may  be  defined  in  many  different  ways.  One  definition  of 
P  (x)  is  as  the  coefficient  of  yn  in  the  formal  power-series  expansion  of  the 
function  (l  -  2  xy  +  y  )"2,  Prom  this  definition  it  is  straightforward^  to 
obtain  the  explicit  formula  for  Pn(x)  given  by  Rodrigues'  formula: 

PBW  =  -J—  ^  (x2-1  )n  .  (B1) 

n  2n  nj  dxn 

This  formula,  of  course,  provides  another  possible  definition  of  Pn(x), 

The  first  six  polynomials  are  given  by:- 


Po«  - 


1»  P^x)  =  x,  P2(x)  =  i(3xZ-l), 


Pj(x) 


K5x3-3x),  P4(x)  =  iOSx'-iO*2*}), 


P  (x)  =  J(63x5-70x3+15x)  . 


The  associated  Legendre  function  I^(x)  provides  a  generalization  of  Pn(x) 
to  which  it  reduces  when  m  =  0,  It  may  be  defined  by 


Pm(x)  =  (l-x2)1^2  -^-P  (x)  , 

nv  '  '  '  „  ra  nv  7  9 


d x 

so  that,  in  view  of  Rodrigues'  formula,  (Bl), 

/.  2\m/2  ,n+m  „ 

Pm(x)  «  - (x2-l)n  . 

n  2  nl  dxn+m 


(B2) 


(B3) 


Some  warning  remarks  on  the  above  definition  are  in  order.  For  although 
there  is  universal  agreement  on  the  definition  of  Pn(x),  authorities  differ 
over  F^(x). 

Jeffreys  and  Jeffreys^  use  equation  (B2)  and  this  definition  is 

incorporated  in  the  1961  recommendation  of  the  International  Astronomical 
28  38 

Union  .  MacRoberp  ,  however,  and  some  older  texts  write 

P^(x)  =  (x^-1  )rn’/^  dm  Pn(x)/dxm.  This  definition  is  obviously  more  profitable 

if  |x[  >  1,  but  less  profitable  if  |x|  <  1,  x  boing  real. 
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Groves^  and  llerson^  have  tried  for  the  best  of  both  worlds,  following  MacRobert 
with  F^(x)  and  introducing  T^(x),  in  the  notation  of  Ferrers"^,  as  a  function 
identical  with  our  P^‘(x).  The  present  author  feels  that  this  double  usage  is 
unnecessary,  contrary  to  the  I.A.U,  recommendation  and  confusing.  A  possible 
source  of  danger  is  illustrated  in  Ref. 3,  which  introduces ^T^(x)  first  and  then, 
when  |x|  >  1,  uses  MacRobert's  P^x),  writing  F^‘(x)  =  T^(x);  but  this  last 

relation  is  wrong,  since  when  x  is  real  and  |x|  >  1,  the  correct  principal-value 
relation  is  (l-x2)1"1^2  =  (-l)1^2  (x2-'!)"1’^2,  leading  to  F^(x)  =  (-l^2  T^(x). 

A  further  point  of  divergence  in  the  definition  of  the  associated  Legendre 
functions  docs  not  concern  us  in  this  paper  but  is  mentioned  for  completeness. 

It  may  be  convenient  to  multiply  pm(x)  by  an  arithmetic  factor  and  to  use  the 
resulting  function,  denoted  say  by  p™(x),  instead  of  F^(x).  Jeffreys  and  Jeffrey^ 
take,  as  this  arithmetic  factor,  (n-m)l/n'.  but  the  I. A. u.  recommend 
{ (n-m) ’./(n+m)  l]  2.  The  object  of  the  factor  is  to  achieve  some  measure  of 
normalization  and  Kaula^  points  out  thiat  for  true  normalization  the  factor 
{um(n-m)l  (2n+1  )/(n+m)ij 2,  where  u^  is  1  for  m  =  0  and  2  for  m  *  0,  should  be 
used.  Kaula's  factor  is  becoming  wide1y  used  in  evaluating  the  earth's  tesseral 
hiarmonic  coefficients^*-1,  but  it  has  one  major  weakness:  the  factor  is  not  unity 
when  m  =  0,  so  thiat  the  basic  Legendre  polynomials  themselves  arc  affected. 

If  x  =  cos  9,  where  0  <  9  <  tc,  the  functions  F^(x)  may  be  expressed  rather 
simply.  For  values  of  n  up  to  3  the  functions,  excluding  Legendre  polynomials 
(m  =  0)  already  given,  are  as  follows: 


P]|(x)  =  sin9,  P^(x)  =  3  sin9  COS9, 

P^(x)  -  sin9  (5  cos29-l),  P2  =  15 


P?(x)  =  3  '  2 


sin  9, 


sin29  C0S9,  P^(x)  =  15  sin^9  . 
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TH3  ADDITION  THCORT&*  TOR  LZG3NDRE  FUNCTIONS 

The  addition  thoorem  ia  tho  key  to  tho  expressing  of  Legondre  polynomials 
relating  to  ono  axis  in  terms  of  Legendre  functions  relating  to  another  axis* 
Thus  in  Pig,  1  (b )  wo  suppose  that  OZ  and  OP'  are  particular  (given)  axes  and 
that  OP  is  a  general  direction;  OP  makes  angles  8  and  9  with  OZ  and  OP' 
respectively,  the  angle  between  these  axes  being  8’,  Then  the  theorem  states 
that  Legcndro  polynomials  in  cos<p  (i,e,  relating  to  0Pf)  are  given  in  terms 
of  Legendre  functions  in  cos  8  (i.e,  relative  to  OZ)  by  the  expressions 

n 

Pn(cos  9)  a  ^urn  fcnjf  ffi008  e* )  ffi008  0)  008  »  (C1) 

m=sO 

where  X  is  the  (azimuth)  angle  between  the  projections  of  CP  and  OP*  on  the 
plane  perpendicular  to  OZ,  um  =  1  if  m  a  0  but  2  if  m  >  0,  and  n  >  0. 

77  7  0 

Proof  of  the  addition  theorem  is  usually  given' *  as  part  of  a  complete 
mathematical  treatment  of  spherical  harmonics.  It  is  possible,  however,  to 
prove  it  by  elementary  algebra,  and  we  proceed  to  do  this. 

The  theorem  is  trivial  when  n  a  0;  when  n  a  1  it  reduces  to 

cos 9  =  cos 81  cos 8  +  sin 8*  sin 8  cos X  (C2) 

which  is  true  by  spherical  trigonometry,  7e  can,  in  fact,  use  this  last 
expression  to  give  a  purely  algebraic  statement  of  tho  theorem:  if  9  is  defined 
by  equation  (C2),  then  (Cl )  holds  in  general. 

It  is  necessary  to  prove  the  theorem  for  n  >  2  and  this  may  be  done  by 
induction.  Thus  we  suppose  that,  for  an  arbitrary  C  >  2,  the  theorem  is  true 
for  n  =  C-2  and  n  =  l- 1,  If  its  truth  for  n  =  •£  can  be  established  on  the 
basis  of  this  supposition,  proof  for  all  n  will  be  complete. 

From  equation  (D5)  of  Appendix  D  with  n  =  £-1,  m  a  0  and  x  =  0039,  from 
the  induction  hypothesis  for  P^(cos9)  and  ani  *’n>m  ecIua‘tion  (C2), 

it  follows  that 


f 
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Of  A 

P^(costp)  =  £  (cos 0 *  cos 6  +  3in0*  sind  cosX) 

6-1 

*  Yj  u»  (oos  6  ’ )  ’x-i (o03  6 >  cos  "*■ 

m=0 

6-2 

•  (C3) 

m=0 

Thus, since 


cos  X  cos  mX 


^(cos  m+1  X  +  cos  ii»— 1 X ) 


(C4) 


P^(oosif)  is  certainly  expressible  by 


P^(coscp)  =  ^f(0',  0,  6,  m)  cos  irX 
m=0 


and  it  remains  to  be  shown  that 


(C5) 


f(0',  0,  6,  m) 


u 

m 


^fl5(oOSe')^(cO30) 


(06) 


The  cases  m  =  6  and  m  =  6-1  are  quite  easy  and  are  considered  before  the 
general  case  (0  <  m  <  6-2).  We  need  -  immediate  from  equation  (B3)  of  Appendix  B  - 
the  expressions 


P^(cos0) 


(2n)i  sin1^ 
2n  n*. 


and 


P^OjosO)  = 
n  v 


sin0-^©  cos0 


n: 


Using  these  expressions  it  may  be  seen  very  easily  from  equations  (C3),  (C4)  and 
(C5)  that 

or  .  sin0*  sin0  pf  ]  (cos  0'  )  P^"*?  (cos  0 ) 

U**l  i/—  I 

C  (26-2)1 

2  r^(cos0f)  P^(cos0) 

mi 


f(0',  0,  6,  6) 
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Similarly, 


f(0',  0,  e,  c-i)  = 


2^  j-2  cosO'  cos  3  P^_^(cos0* )  P^“^(co3  0) 
~T~  l  (2C-2 ) 

sin©'  sin0  P^_^(cos  d 1 )  P^,_^(cos  ©Xj 

+  (26-3)1  J 


2  P^— 1  (cos  0  * )  P^”1  (cos  0 ) 

=  (2£-l)>. 

For  the  general  case,  with  0  <  m  <  £-2,  bearing  in  mind  the  special  case 
when  m  =  C, 


f(0’,  0,  £,  m)  = 


1  [c°se'  coae  um  (eXl?l  E^.1(°»s6,)rt.1  (cose) 

line-  sir.6  (u-l)  fT]  (cos  8 ' )  (°os6 ) 


6’  3i"6Umv,  IX^fC!(00s9')P”1l('!Me’] 


+  sin 


fe-I )•  rc-2(=ose,)Pe-2(oo30) 


Now  P^*^  may  be  eliminated  by  use  of  equation  (D3)  of  Appendix  D  expressed  in 
the  appropriate  form  with  P^~]  on  the  K.H.S.  Similarly,  may  be  eliminated 

by  use  of  equation  (P6),  with  n  =  £-1.  This  leads  to 

f(0',  6,  £,  m)  =  g^  cos0*  cos©  P^^cosO'  )  F^‘  ^(cos0) 

-  g2  {cos0'  P^^cosOOP^^cosO)  +  cos 0  P^2(cos 0 1 )  F^^cosO)] 

+  g3  P£-2^cos0^  C2(c°86)  > 
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where 


Sh 


26-1 

C 


fu  ) »  .  Um**^  (6-m) 1 _  i  u  (6-m-1  )2 

I  m  te^TTt  .  N2  (-6+111-2)1  *  m+1  (7+mJl 

v-  (6-m)  ^ 


.(££.-,0. (b3)l  ju  (e2-m2)  +  (u  -l)  (6+m)  (6+m-l)  +  £  u_,(6-m)(6«n*-l)) 

(6-m'l2  fiijnll  *  m  ^ 


6 (6-m)  (6+m)  l 

u  (26-1 ) 2 


m 


(6-m)‘ 


(6-m)  l 

(6+m)  I  * 


«2  = 


26 

6 


-i  rvl_  (6h 

u  (2£-1  )(C+m-1 ) 
m 


(6-m)‘ 


(6-m)  \ 
(6+m)  J 


and 


g; 


26- 


*  [(6-m)2  +  ^  Um+1 

u  (6+m-l)2 


u  _i 

in 


(feaH,. .  a  „  )2  .  £=iu  (*-»?); 

^+m  (6+m)  I  J  6  m  (6+m-2)I 


(6-m)' 


(6-m)! 
(6+m)  l 


Hence 


f(8J,  0,  C,  m)  =  »B^j 


,  (26-1 )  cos0f  (cos  0')-( -6+m-l)  P^_2  (cos  ©  * ) 

6-m 

(26-1 )  cos0  (cos  0)  -  (6  +m-1  )  *SL  O(cos0) 


6-m 


=  um{fi$H”(ooae')re(oosfs) 


on  using  equation  (D5), 


This  completes  the  proof  of  the  addition  theorem. 
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RECURRENCE  RELATIONS 

The  purpose  of  this  Appendix  is  to  prove  the  important  relations,  for  the 
A,  D,  C,  D  and  E  functions,  given  by  equations  (50)  to  (57)  of  the  main  text* 
The  starting  point  is  the  definition  of  Legendre  polynomials  by  means  of  the 
identity 

CO 

(1  -  2  xy  +  y2)“J  s  ^  yn  Pn(x)  , 

0 

from  which  it  follows  at  once  that 

00  2 

(l  -  2  xy  +  y2)  QTV  pn(x)^  a  1  • 

0 

On  differentiating  both  sides  partially  with  respect  to  x  and  y  in  turn  the 
resulting  pair  of  relations  leads,  after  some  manipulation,  to  the  classical 
recurrence  formulae 


>p;,W  -  p;(x)  =  (n+l)  Pn+1(x)  (D1) 

and 

P^+1(x)  -  x  P'(x)  =  (n+l)  Pjx)  (D2) 

for  n  >  0. 

These  last  two  formulae  may  be  differentiated  m  times  with  respect 
leading  to: 

xp^l)(x)  -  p<m+1>(x)  =  (n-m+l)  P^j](x) 

and 

P^l5(x)  -  xpj;m+l)(x)  =  (ntm+1)  P^m)(x).  (Bl) 

We  replace,  in  equation  (D3),  n  by  n-1  and,  in  both  equations,  m  by  m-1. 

After  the  elimination  of  )(x)  and  the  introduction  of  the  Legendre-f unction 
notation  by  equation  (B2),  it  then  follows  that 


to  x, 
(D3) 
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(n-m+l)  *^+1(x)  -  (2n+l)x  rj(x)  +  (n+m)  P^-()(x)  «  0  •  (D5) 

Equation  (D5)  is  the  key  recurrence  relation  and  holds  for  0  <  m  <  n-1,  However, 
a  second  relation,  proved  similarly,  is  also  given,  since  it  is  used  in 
Appendix  C: 

(1-x2)^  I^+1(x)  s  (n+m)  P^(x)  -  (n-m)x  F^(x)  ,  (P6) 


V/e  can  now  proceed  to  derive  the  required  recurrence  relations  for  the 
A,  B,  C,  D  and  E  functions. 

Pram  equations  (17)  and  (l8)  of  section  4*2  of  the  main  text, 

sinki  Ak(i)  =  2k  ^k||  k»,  f£(cos  i)  ,  (D7) 

Hence  from  equation  (D5),  with  n  =  l-\  ani  m  =  k, 

♦  aL>w  -  °  • 


The  required  relation  for  Ak(i)  follows  at  once  and  is  given  by  equation  (50)  of 
section  6,2, 

k 

For  the  B^(e)  function,  it  follows  from  equation  (31 )  of  section  4*2  and 
equation  (El  3)  of  Appeniix  E  that 


Bfc>  ■  ^  (})*^  (-1^4,00 


where  £ 


Hence  from  equation  (D5),  with  n  =  h-2  ani  m 


k. 


(D8) 


(h-k-1 ) 


,h-1 


*  (h+k-2){h^5h'3Bh-2(e)  ■  0  • 

The  required  relation  for  B^(e)  follows  at  once  and  is  given  by  equation  (51 ) 
of  section  6,2, 


Appendix  D 
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For  (with  l  =  n  if  n  >  0  and  0  =  -n-1 
equation  (34)  of  seotion  1*2  and  equation  (B3) 
even, 


if  n  <  0)  we  have,  from 
of  Appendix  B,  that  if  t- k  is 


(-D 


•  •  •  #  (n-k) 

22M  (k,)2  [i(e+k)|l 


(D9) 


while  if  0-k  is  odd. 


0  . 


(DIO) 


Considering  only  the  case  0-k  even,  we  have,  by  dealing  with  the'n  factors' 

and  the  *0  factors'  separately,  from  the  'n  factors'  a  contribution  to 

Ck/Ck  o  of 
n  n-2 

-  (n-1 )( n-2) 
n-k-1  )(n-k-: 


T 


27 


and  from  the  *0  factors'  a  contribution  to  or  (accort*in8 

as  n  =  0  or  n  s  -0-1 )  of 


0+k-1 

0-k 


Hence  if  n  =  0  (i*e,  for  n  >  0), 


ck/ck 

n  n-2 


( n-1 )  (n-2 )  ( n+k-1 ) 
(n-k) (n-k-1 ) (n-k-2 ) 


and  if  n  =  -0-1  (i.e.  for  n  <  0), 


Ck/Ck  9 
n  n+2 


( _n+k-2 )( -n+k-1 )(-n+k) 
( -n) (-n-1) (-n-k-1 ) 


These  are  the  same  as  equations  (32)  and  (54)* 

From  equation  (D9)  it  follows  that,  for  k  >  0, 
k  uk  (2k+2) J 


kt2  "  2*+2(k!)2 


and 


jc  .  ,  ,'k-i  V210!  (2^1)'. 

~k'1  2,k(kJ  )5  (ktl);  ’ 


from  which  equations  (53)  and  (55)  may  be  derived  at  once* 
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For  the  function  D*(i),  there  is  no  three-term  recurrence  relation  when 
k  >  1#  Suoh  a  relation  exists  for  the  case  k  =  0,  however,  and  this  case  is 
the  most  important.  From  equation  (46), 


d  P^(cos  i) 
f  d(cos  i) 


s  -sin  i  P^(cos  i) 


(Dll) 


The  required  relation  follows  at  once,  using  equation  (D5),  and  is  given  by 
equation  (56),  It  is  useful  to  note  that  the  leading  term  of  D°(i)  is 
~y?C(C+ l)f  oosiif  Z  is  even,  and  -46(6+1  )f  if  6  is  odd;  the  leading  term  of 
D^(i),  for  k  >  1,  is  k  cos  i  if  6-k  is  even,  and  k  if  6-k  is  odd. 

For  E^(e),  there  is  again  no  three- term  recurrence  relation  unless 
k  s  0,  But  it  is  shown  in  Appendix  E  that  for  n  >  0 

E°(e)  =  -  (-1)*  e(l-e2)^  P„(«  .  (D12) 

Hence  E^(e)  satisfies  a  relation  similar  to  that  which  holds  for  D°(i);  it  is 
given  by  equation  (57). 

For  n  <  0  there  is  no  need  to  use  a  recurrence  relation;  for  it  follows 
from  Appendix  E  that,  when  n  <  0,  E°(e)  is  identical  with  E0^  ^(e).  It  mist 
bo  noted  that  there  is  no  similar  result  when  k  >  1, 

k  2 

Leading  terms  for  En(e)  are  ^n(n+l)e  for  k  =  0,  n  >  0,  and  k  for  k  >  1 
(and  any  n). 


i 


u 


I 
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PROPERTIES  OP  THE  FUNCTIONS  S  ,  (e)  AND  B^e) 

nkv  nx  / 

The  main  result  of  this  Appendix  is  the  establishment  of  equation  (32)  of 
the  main  text,  relating  S^e),  for  n  <  0,  with  B^^Ce)*  After  this,  simple 
mathematical  expressions  for  S^e)  are  obtained,  covering  all  values  of  n. 
Finally,  two  relations  for  En(e)  are  obtained,  those  given  by  equation  (D12)  of 
Appendix  D,  and  equation  (58)  of  the  main  text. 

Instead  of  proving  equation  (32)  itself,  it  is  convenient,  first,  to 
extend  the  definition  of  B^(e)  to  cover  n  <  0,  In  this  case  equal  to 

(n-1 )(n-2)  ,,,,  (n-k)Al,  is  non-zero  and  wo  define  B^(e)  by  the  first  of 
equations  ( 31  )#  For  the  purpose  of  symmetry  in  expressing  the  result,  we  use  a 
slightly  different  range  of  k,  employing  h(=n  if  n  >  0  and  -n+1  if  n  <  0). 

Then  the  result  to  be  proved,  for  all  n  and  k  <  h-1,  is 

B*(e)  =  (1-e2)"-i  B^n+1(»)  .  (El) 


Since  the  starting  point  is  the  function  S^e),  which  is  defined  only  as 
a  series,  the  proof  is  not  short.  We  divide  it  into  stages  and  prove  four 
lemmas. 


Lemma  1  B^e) 

-  rr  ' 

for  k  <  h-1,  where  F(a,  b; 

F(a,  b;  c;  x) 


-n+k+1  -n+k+2  #  , 

»  Q  i 


2  '  2 

c;  x)  is  the  hypergeometric  series  defined  by 


.  a.b  ..  a(a+l)  b(b+l)  „2 
=  1  +  0X+  o(oil)l.2  *  +  ♦*" 


(E2) 


oo 


=  r_(o)  VXa+q)  r(b+q)  q 

r(a)  r(b)  r(o+q)  q! 

q=0 

and  the  gamma  function  T(x)  is  related  to  the  factorial  function  by 
r(x)  =  (x— 1 ) , 

This  lemma  can  be  proved  by  considering  the  general  term  of  the  series 
S^e)  given  by  equation  (29),  We  have 
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u  S 

q 


(n-1 ) (n-2) 

(k+2q)' 


ln-k-2q) 


(n-l)  ....  (n-k)  -n+k+1  -n+k+2 
k».  2  2 


-n+k+2 q  kl  e^ 

2  (k+q)i  qi 


(E3) 


Taking  the  faotors  £(-n+k+l),  £(-n+k+2)  eto  alternately  in  two  sets  we  may 
now  write 


■  U)' - 


-n+k+1  ’\  „/-n+k+2  ' 

+  q)  -n(— 2 —  +  q 


-n+k+1  ^  p^-n+k+2^j" 


Xflstll  (e2)^  . 

1  (k+1+q)  '  '  ' 


i.e. 


•  (?)<: 


-n+k+1  -n+k+2 

> -  »  — § -  ;  k+1 


(E4) 


from  which  the  result  follows  using  equation  (31  )• 

Lemma  2 

For  an  arbitrary  set  of  m  quantities  v. ,  v^,  ....  v^,  where  0  <  m  <  q-1, 
we  have 


2^  (-0s  Q  (v.j-sJCvg-s) 


(vs) 


=  0 


s=0 


(E5) 


This  lemma  is  easily  proved  by  observing  that,  if  =  (v^-s)(v2~s)  .... 
(Vs),  then  wg  is  an  mth  degree  polynomial  in  s.  Since  q  >  m  the  qth  order 
differences  of  this  polynomial  vanish.  But  these  differences  are  simply 


SbO 


Lemma  3 

If  G(a,  b, 

G(a,  b,  c,  q)  = 


c,  q)  is  defined  by 

^  (s)  a(a+0  •••• 
SbO 


(a+q-1-s)  b(b+l)  •  •  •  •  (b+q-1-s)(c+q-s)  •  ••• 
(c+q-1 )(c-a-b)  f  (c-a-b-1+s)  ,  (E6) 


then 

G(a,  b,  c,  q)  s  (c-a)(c-a+l)  ....  (c-a+q-1  )(c-b)(c-b+1 )  ....  (c-b+q-l).  (E7) 
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This  lemma  may  be  proved  by  checking  that  each  factor  on  the  R.H.S.  of 
equation  (E7)  ia  a  factor  of  G(a,  b,  c,  q).  The  result  will  follow  at  once 
since  there  are  the  same  number  (2q)  of  factors  on  the  R.H.S,  as  in  each  term 
of  G(a,  b,  e,  q)  while  the  term  r.  with  unit  coefficient,  occurs  on  the 
R.H.S,  of  equation  (E7)  and  ii  the  last  term  of  G(a,  b,  c,  q). 

To  check  that  c-a+p,  for  0  <  p  <  q-1,  is  a  factor  of  G(a,  b,  c,  q)  we 
show  that  G(a,  b,  a-p,q)  =  0,  Now  the  q+1  factors  a,  a+1,  •  •••  a+q-1-p, 

b,  b+1,  b+p  occur  in  every  term  of  G(a,  b,  a-p,  q ).  Dividing  them  out, 

there  remains  the  series 


^  (-l)S  (fj  (a+q-p-s) 
s=0 


•  •  •  • 


(a+q-1-s)(b+p+1-s)  ....  (b+q-1-s) 


which  vanishes  by  lemma  2,  with  m  =  q-1.  Thus  c-a+p,  and  similarly  c-b+p,  is 
a  factor  of  G(a,  b,  c,  q). 

Lemma  A- 

P(a,  b;  c;  x)  xF(c-a-b,  1;  1;  x)  =  F(c-b,  c-a;  0;  x)  .  (E8) 


On  multiplying  out  the  two  series  on  the  L.H.S. ,  it  may  be  seen  that  the 
coefficient  of  x^  is 


SrO 


a (a+1 )  ....  (a+q-1-s)  b(b+1 )  ....  (b+q-1-s) 
c(c+l)  ....  (c+q-1-s)  (q-s)i 


x 


(c-a-b)(c-a-b+1 )  «...  (c-a-b-1+s) 
si 


} 


which  is  simply 


c 


b.  c.  a) 

....  (c+q-1) 


q 


f 


or  again,  from  lemma  3, 

(c-a)(c-a+l)  ....  (c-a+q-l)(c-b)(c-b+l)  ....  (c-b+q-1 ) 
c (c+1 )  ....  (c+q-1;  q» 

Hit  this  last  quantity  is  the  coefficient  of  x^  in  F(c-b,  c-a;  c;  x).  This 
completes  the  proof  of  lemma  4. 
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0 


The  proof  of  the  main  result  of  this  Appendix  is  now  quite  straightforward. 
In  lomma  4  wo  take 


-n+k+1 


-n+k+2 


»  =  — r~  >  b  =  2 


,  c  s  k+1  ,  x  =  e  , 


This  gives 


-n+k+1  -n+k+2 


;  k+1  ;  e2^  x  F(i>4,  1?  1;  ©2)  =  ,  ^^Y^  J  k+1  j  e2^ 


....  (E9) 


But  F(n-4»  1;  1;  e2)  =  (l-e2)”n+?;  lemma  1  and  equation  (E9)  then  give 


l£(e)  x  (l-e2)"n+2  =  B^n+1(e)  . 


n 


(E10) 


This  is  the  result  required  . 


Next,  in  this  Appendix,  we  give  explicit  expressions  for  ^^(e)  for  all  n 
and  k  <  l. 

First,  if  n  <  0,  then  from  equation  (E3), 


(-n+k+2q)»  _ 1 _ 

(-n)l  (k+q)l  ql 


,  \k  l(2k+lU  ,(,r.h^+2qU  _ (21a.n.2.q).i,  kJL-  , 

1  (2k+1+2q)J  2^^(2k+l)l  ql(k+q)'. 


2q 


But 


and  so 


(2k+1+2q)j  kl 
22<1  (2k+l)j  ql(k+q)l 


(2k+3)(2k+5)  ....  (2k+1+2q) 
21  qf 


* 
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i.o. 


S  J®) 
nkv  ' 


a-n-k-1 

-n-k-1 


de 


(Ell) 


If  n  >  0,  -n-k-1  will  be  negative  and  the  above  formula  breaks  down*  For  n  =  0, 
k  is  restricted  to  0  and  we  have,  easily, 


Soo(e)  =  (1-e2)'*  .  (E12) 

Finally,  if  n  >  1  we  can  prove  a  result  obtained  by  Groves3,  in  a 
different  notation,  viz. 

Snk<e>  -  T^(ff  S1"1  (‘l)3k/2  £ltt)  '  (E13) 

2  — —  3k/2 

where  Z,  =  (l~e  )  2  (the  factor  (-1  y  '  conforms  .vith  a  remark  in  Appendix  B). 
We  need  a  lemma 


Lemma  5 


( 


-n+1+k 


+  s, 


-n+2+k 


+  *1  l)  »  n^-k-a;(2n~2w}J - T  *  (E14) 

/  2n  1  k  (n+k-1 )  l(n-l-s)  I 


This  follows  fairly  easily  from  Gauss’s  theorem  for  the  hypergeometrio 


38 

function,  viz 


'<*.  **  •*  D  ■  swass}  « 


(E15) 


on  taking  a  = 


-n+1 +k 


-n+2+k 


+  8,  b  s  — 2 —  +  s,  c  =  k+1+s.  For  then 


F(....aain(E14)....)  =  ^^^3^  j, 


/.  x.  (2n-3-2s)  1  2n+lc"1 

_ (k+s)i  (2n-2-2s)l 

2n-1-k"28(n+k-l)’(n-1-s): 

Having  established  tho  lemma  we  proceed  to  obtain  the  required  expression 
for  Prom  equation  (E3), 
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I 

$ 


w 


(n-1 ) L 


*  (n-1-k-2q)i(k+q)i  q'.  22q 


|l  -  (l-e2))q 


and  ao 


Snk(e)  ■  Y - - 2 of  (1)  (-1)S  «1-°2)B 

1—1  (n-1-k-2q)  i(k+ql)ql  2  q  l—i  \) 


(El  6) 


q=0 


8=0 


where  wo  oan  take  the  q  3eries  as  infinite,  even  though  it  terminates,  since 
(n-1-k-2q) i  =  oo  when  n-1-k-2q  <  0« 


But 


oo  q  CXs<q  oa  oo 

It-  l-ll 

q=0  s=0  q,s  s=0  q=s 


ani  hence 


3  (e)  =  V  V  (-l)3(n-l);(l-a2l3 _ 

^  Li  Li  (n-1-k-2q)i(k+q)  i  s'.(q-s)[  22q 


3=0  q=s 


=  \  (-l)S(n-l)j _ (V _ (n-!,-k-23)i,.(kt3j.; - J  (l-e2)3 

l—i  (n-1-k-2s)i(k+s)  i  si  22s  \J—i  22q”2s(n-1-k-2q)  '.(k+q)i(q-s)  'J 


s=0 


q=S 


(El  7) 


The  quantity  inside  the  central  brackets 


-n+1+k 


+  s. 


-n+2+k 


+  s;  k+1+s; 


;1> 


So 


may  easily  be  seen  to  be 
from  lemma  5* 


c 


ii 

ti 


I 
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00 


■«<•> 


■I 


SsO 


2\a 


(.l)a(n-l)l(2n-2-2s)l(l-eg) 
(n+k-l^n-l-k^sMn-l-s)!  s!2 


n-1-k 


■  fes#  S^E('i)8^^(n;1)5 


n-1-k-2s 


8=0 


_  (n«l)l  /eV^I-n  (&e)k 

'  {2j  4  z^Vl)! 


,  n+k-1  * 


(El  8) 


2  2  2 

Since  £  e  =  4-1,  the  required  expression  for  S  ,  (e)  follows  at  once  from 


equation  (B3)  of  Appendix  B,  noting  that  (^-l)1^  _  (l-^2)1^2. 

The  remaining  results  proved  in  this  Appendix  relate  to  the  function 
E^(e)  in  the  particular  case  when  k  =  0.  First  we  obtain  equation  (D12)  of 
Appendix  D. 

From  equations  (47),  (31 )  and  (E13),  when  n  >  0, 


K>) 


=  (2n-l)  e2  B°(e)  +  e(l-e2)  Jj-B°(e) 

=  e(l-92)n+^f(l-e2rnB°(e)i 

-  «(i-«2)n4^fsn  Vl(S)i  • 


(E19) 


3  p  -1 

But  d£/de  =  qw  since  £  =  (l-e  )~?j  also,  from  equation  (D2),  with  n  replaced 
by  n-1, 

Sn  P'  .  (S)  +  n  S1"1  P„  .  (S)  =  S'"1  P!  . 
n— i  n— i  n 


Hence 


E°(e)  =  82(l-e2)*('"l)P^(S)  .  (E20) 

Using  equation  (B2)  to  introduce  P^(£),  the  result  given  by  equation  (D12) 
follows  at  once. 

The  other  result  -  equation  (58)  of  section  6  -  is  best  proved  by 
0 

expressing  E^(e),  for  n  <  0,  in  terms  of  P^^fe),  i«e#  in  terms  of  gO®)* 


P 

!>> 
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Prom  equations  (48)  and  (E13),  if  n  <  0, 


E>>  -  e  air  Bh<e) 


d  (s/1-h 


=  ef  [5 

de 


(E2l) 


But  d£/de  =  e£^  and  from  equation  (Dl),  with  n  =  h-2, 


S1_h  Pj^fe)  +  (l-h)  5"h  PW(S)  =  «'h  P^(S) 


(E22) 


Hence 


E  (e)  = 

nx  '  0 


e2(l-e2)^h_3)P^2(S) 


=  -(-l)*  e(l-e2)i(h_2)  P^_2(C) 


=  »£*<•> 


(E23) 


by  equation  (D12),  just  proved. 


This  is  the  result  required. 
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AXIS  TRANSFORMATION  RELATIONS  REQUIRED  BY  THE  THEORY  OF  PERTURBATIONS 
DUE  TO  A  STATIONARY  DISTURBING  BODY 

In  this  Appendix  we  obtain  the  formulae  required  by  section  8.3*  Four 
systems  of  axes  are  relevant  and  are  illustrated  in  Fig.  3.  The  fixed  system 
Oxyz  has  been  defined  in  section  8*3,  Ox  pointing  to  seme  instantaneous  position 
of  the  node  and  Oz  being  perpendicular  to  the  orbital  plane.  The  system 
OXYZ  relates  to  the  (true)  equator,  OZ  being  towards  the  north  pole  and  OX 
coincident  with  Ox.  Since  it  is  only  the  perturbation  in  0  that  is  of  interest, 
its  origin  is  irrelevant  and  there  is  no  need  to  introduce  the  direction  of 
the  first  point  of  Aries.  The  systems  Ox'y’z*  and  OX’Y’Z*  are  defined  in  an 
exactly  analogous  manner  but  relate  to  the  pseudo-equatorial  plane  which  is 
perpendicular  to  the  direction  of  the  fixed  disturbing  body. 

Then 


where 


~1 


0 

cos  i 
sin  i 


and 


h 


(cos  a  -sin  a  0 

sin  a  cos  a  0 

0  0  1 


/I  o 

0  cos  i* 
\o  -sin  i' 


0 

•  •  • 
sxn  x* 

cos  V 


We  can  at  once  derive  equation  (74)  of  section  8.3.  For  A,  B,  C  are  the 
direction  cosines,  in  the  Oxyz  system,  of  the  Z1  axis,  so  that 


(pi) 
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To  obtain  equations  (75) ,  (76)  and  (77)  we  observe  that  the  (small)  changes 
Ai,  AQ  and  Aw  represent  the  rotation  of  the  orbital  plane  from  one  orientation 
to  smother.  The  total  rotation  is  in  fact  compounded  of  rotations  Ai  about  Ox, 

Afl  about  02  and  Aw  about  Oz  and  hence  has  corrponents,  in  the  0 xyz  system,  given 
by  (Ai,  AO  sin  i,  An  cos  i  +  Aw),  using  R  ,  In  the  Ox’y'z*  system  the  components 
of  the  total  rotation  are,  similarly,  given  by  (Ai’,  Anf  sin  i',  A O'  cos  i’  +  Aw'), 

The  preceding  paragraph  depends  on  the  fact  that  small  rotations  are 
vectors.  Using  this  fact  again  we  get:- 


(F2) 

Equations  (75),  (76)  and  (77)  follow  at  once,  an  introducing  A,  B  and  C  from 
equation  (74), 


■  h 


Ai* 

sin  i  '  AO  ‘ 
lAw*  +  AO'  cos  iV 


c 
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RELATIONS  REQUIRED  BY  THE  THEORY  OF  PERTURBATIONS  DUE  TO  A 
DISTURBING  BODY  IN  A  KEPLER  ORBIT 

The  relations  required  are  those  given  by  equations  ( 126)  to  ( 1 35)  of 
section  8,4,  It  is  convenient  to  define  altogether  six  systems  of  axes  and 
these  are  illustrated  diagrammatical ly  in  Pig, 4*  The  systems  OXYZ  and 
OX^YjZ^  are  based  on  the  equatorial  plane,  intersected  by  the  satellite  orbit 
and  the  disturbing  body  orbit  respectively.  The  systems  Cbcyz  and  Ox"y,,z" 
are  based  on  the  satellite  orbital  plane,  intersected  by  the  equator  and  the 
disturbing  body  orbit  respectively.  The  systems  and  OX'T'Z"  are 

based  on  the  orbital  plane  of  the  disturbing  body,  intersected  by  the  equator 
and  the  satellite  orbit  respectively. 


All  the  relations  derive  from  transformations  between  these  systems  of 
axes.  As  in  Appendix  P,  the  first  three  -  given  by  equations  ( 126)  to  (l28)  - 
follow  from  expressing  the  angular  momentum  vector  for  the  satellite  orbit 
in  each  of  the  systems  Oxyz  and  Ox,,ywzM,  and  transforming.  Thus 


/  dl 

1  dt 


,  .  dfl 
3in  1  dt 


dw  .  d? 5 

\at  +  003  1  at 


/ 


cos  a  -sin  a  0 


sin  a  cos  a  0 


\  0 


dl" 

dt 

sin  i" 


dfl" 

dt 


du"  d?5" 

idT 4  003  1  aT 


The  important  relation  given  by  equation  (129)  expresses  the  transformation 
of  axes  between  systems  Oxyz  and  Ox^y^z^,  carried  out  via  the  auxiliary  systems 
OXYZ  and  OX^Y^Z^.  The  relation  follows  from  the  equations  for  these  auxiliary 
transformations,  viz. 


/I  0 

0  cos  i 

\0  -sin  i 


sin  i  | 
cos  i/ 


/  cos(n-nd) 
-sin(0-0d) 
\  0 


sin(n-nd) 

cos(n-nd) 

o 


o ' 
o 

i/ 


In  \ 

t  1 
d  1 

\z  ! 

v  d' 
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< 


and 


| 

$ 


i: 


o 

cos  i, 
a 

sin  i, 
a 


0  \ 

-sin  i, 
a 

008  tj 


As  in  section  8,4  we  abbreviate  the  full  transformation,  obtainod  by  multiplying 
out  the  above  matrixes,  to 


/y(  _ 

X  \  » 

lL  1 

M1 

V 

(Xi\ 

y 

I1"2 

“2 

N2 

iyd 

• 

(Gl) 

w 

\L3 

m3 

V 

w/ 

liquation  (l30)  of  section  8,4,  viz  w"  =  w-a,  is  immediate  from  the 
definition  of  a  as  the  angle  in  the  orbital  plane  of  the  satellite  between 
Ox  and  0**. 

To  obtain  equations  (l3l)  to  ( 1 33 ) ,  we  express  the  direction  cosines  of  the 
axis  Oz^,  identical  with  the  axis  OZ" ,  in  the  system  Oxyz  by  two  different 
methods:  first,  using  equation  (Gl),  and  second,  transforming  via  the  system 
Ox"y,,z".  Thus 


|  =  /cos  a 

-sin  o 

h 

0 

0 

l°\ 

N2 

I  sin  a 

cos  a 

0 

!° 

cos  i" 

sin 

i" 

1  0 

\N3 1 

1  \  o 

0 

i/ 

\° 

-sin  i" 

cos 

i'V 

\1  / 

-sin  i"  sin  a  ^ 
sin  i"  cos  a  j 
cos  i"  I 


Finally,  to  obtain  equations  ( 1 34)  and  ( 1 35 )  we  express  z,  identical  with  z", 
in  terms  of  x^,  y^,  z^  in  two  different  ways:  first,  using  equation  (Gl ),  and 
second,  transforming  via  the  system  0X"Y"Z".  Then 


(Lj  N3) 


(0  -sin  i"  cos  i" )  /  cos  0" 

| -sin  0" 


sin  fl" 
cos  0" 
0 
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On  taking  components,  we  have  (in  addition  to  a  second  proof  of  equation  (l 33) ) 

=  sin  i"  sin  Q" 


and 


=  -sin  i"  cos  0"  • 


But 


0" 


w. 


and  hence 


L,  cos  w 
3  < 


sm  co. 


-sin  i" 


sin  0, 
a 


and 


L,  sin  u 
3  i 


cos 


sin  i"  cos  0^  . 


These  are  the  same  as  equations  (l34)  and  ( 1 35  )• 


0 


1 

cos  i 


Table  1 


The  function  A„(i) 


COS  1 


cos  i  ^1  -  ^  f 
\  > 


1  -  5  f  +  ^  f 2 


cos 


cos  1 


1-i' 


COS  1 


Table  2 


The  function  B,  (e) 


\ 

h\ 

0 

1 

2 

3 

4 

3 

1 

1 

0 

2 

1 

1 

0 

3 

1  +  2e2 

1 

1 

0 

4 

1  +  |-  e2 

1  +  f  e2 

4 

1 

1 

0 

5 

1  +  3e2  +  |-  e* 

— k 

+ 

© 

IV) 

1  +T  e2 
o 

1 

1 

0 

113 


SYMBOLS 


Symbol 


first  appearance 


Definition 


A,B,C 


a*(d 

=£(•) 


D»(i) 

Ek(e) 


H^e>^ 


3 

P  (  ) 
nv 

E^(  ) 

nv 


9 

10.2 

10.2 

10.2 


10.2 

10.2 


8*w 


direction  cosines  of  disturbing  body  in 
x,y,z  axis  system 

function  of  i  given  by  equation  (l8) 

function  of  e  given  by  equation  (31 )  or  (33) 

constant  given  by  equation  (34) 

function  of  i  given  by  equation  (46) 

function  of  e  given  by  equations  (47)  and  (48) 

solar  radiation  force  per  unit  mass 
given  by  e^^/sin  k  (w-^x) 

given  by  fi^^cos  k  (w-^x) 

given  by  u^^/cos  k  («-^x) 

gravitational  coefficients  of  aspheric ity 

given  by  c^/cos  k  (w-^x) 

given  by  “^^2(^+1)  cosk  (<*>-|x) 

cos  (n-n^) 

cos  id  sin  (Q-fi^) 

-sin  i^  sin  (fi-Q^) 

mean  anaiAaly 

-cos  i  sin  (0-0, ) 
a 

sin  i^  sin  i  +  cos  i^  cos  i  cos  (0-0^) 

cos  i^  sin  i  -  sin  i^  cos  i  cos  (0-0^) 

3 in  i  sin  (0-0,) 
a 

sin  i^  cos  i  -  cos  i^  sin  i  cos  (0-0^) 

cos  i^  cos  i  +  sin  i^  sin  i  cos  (0-0^) 

Legendre  polynomial 

Legendre  function 

some  fixed  value  of  r,  normally  the  mean 
equatorial  radius  of  the  earth  or  the  distance 
to  a  disturbing  body 

function  of  e  given  by  equation  (29) 

orbital  period 


«MM W!w~. 


114 


Symbol 


U 

U 


n 


X,Y,Z 


Soctlon  of 
first  appearance 

3 

4.1 

4.1 

Appendix  P 


X',Y',Z»  Appendix  P 


X”,Y”,Z”  8.4 


a 

a* 

b 

d 

e 

f 

h 

i 

k 

C 

n 


3 

7 

4.2 

8.2 

3 

4*  2 
4.2 
3 

4.2 

4.1 

3 


n  4.1 

n»  7 


P 

r 

t 

\ 

ua 

V 

x,y,z 


42 

3 

3 

42 

8.4 

4.2 

8.3 


x',y ',2*  Appendix  P 
x",y”,z”  8.4 


SYMBOLS  (Contd) 


Definition 


disturbing  function 

1 

term  of  U  containing  the  factor  r 

term  of  containing  the  factor  cosky 

axes  based  on  the  equator,  with  OX  towards  the 
node  of  the  satellite  orbit 


axes  based  on  the  plane  perpendicular  to  the 
direction  of  a  fixed  disturbing  body  as  pseudo¬ 
equator,  with  OX’  towards  the  satellite  pseudo- 
node 

similar  to  X'jY’jZ1  but  with  the  orbit  of  a 
disturbing  body  as  pseudo-equator 
semi-major  axis 

constant  related  to  a  and  U  by  equation  (62) 

a  or  p  according  as  n  <  0  or  n  >  1 

suffix  for  disturbing  body 

eccentricity 

2 

sin  i 

-n+1  or  n  according  as  n  <  0  or  n  >  1 
orbital  inclination 
integer  in  the  range  0  <  k  <  C 
-n-1  or  n  according  as  n  <  0  or  n  >  0 
mean  motion 

index  associated  with  component  U  of  U 

2  3  n 

constant  such  that  n*  a’  =  p 
a(l-e2) 

distance  of  satellite  from  centre  of  force 
time 

1  if  k  =  0,  but  2  if  k  >  0 

9d +  va 

true  anomaly 

axes  based  on  the  satellite  orbital  plane,  with 
Ox  along  OX 

axes  based  on  the  satellite  orbital  plane,  with 
Ox?  along  OX  ’ 

axes  based  on  the  satellite  orbital  plane,  with 
Ox”  along  OX” 


1 


SYMBOLS  (Contd) 


Symbol 


Section  of 
first  appearance 


Definition 


a 

P 

Y 

A 

6 


8.3 

4.1 

4.2 
5 

10.1 


e 


4.2 


5 

5 

5 


8.4 

3 


P 

cr 


5 

3 


X  3 

*  5 

n  3 

ft  10.1 

<j  3 

u  10. 1 

i’,n»,w»  8.3 

i",n",u>"  8.4 


angle  between  x  and  x'  axes 
geocentric  latitude 
(w+v) 

perturbation  for  a  complete  satellite  revolution 
perturbation  other  than  for  a  complete  satellite 
revolution 

0  if  n  <  0  but  j  if  n  >  1 
representative  orbital  element 
given  by  equation  (36) 

term  of  £n  associated  with  k 
a 

gravitational  constant  for  primary  centre  of 
force 

auxiliary  element  introduced  by  equation  (43) 
modified  mean  anomaly  at  the  epoch,  used  instead 
of  x 

mean  anomaly  at  the  epoch  t  =  0 

auxiliary  element  introduced  by  equation  (42) 

right  ascension  of  the  ascending  node 

rate  of  secular  motion  of  n  associated  with  J ^ 

argument  of  perigee 

rate  of  secular  motion  of  w  associated  with  Jg 
pseudo-elements  associated  with  the  axes 

X’jY'jZ* 

pseudo-elements  associated  with  the  axes 
X" ,Y",Z" 


The  above  list  does  not  include  symbols  introduced  temporarily,  with  a 
limited  use,  Zome  of  the  symbols  in  the  list  have  other  uses  of  this  type. 
The  symbol  n  is  used  widely  with  both  its  meanings,  but  no  confusion  should 
arise. 
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FIG. I  (a)  A  SINGLE  POINT  MASS 
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FIG. I  (C)  EXAMPLE  WITH  A  LINE  MASS 

FIG.  I  SIMPLE  MASS  CONFIGURATIONS  OF  WHICH  THE 

POTENTIAL  IS  STUDIED 
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FIG.2  VECTORS  FOR  PERTURBATION  EFFECTS 
DUE  TO  A  DISTURBING  BODY 


OXYZ  BASED  ON  EQUATOR  WITH  OX  ALONG  INTERSECTION 
WITH  SATELLITE  ORBIT  PLANE 

OX'V'Z'  BASED  ON  PLANE  PERPENDICULAR  TO  DIRECTION  OF 
DISTURBING  BODY,  OX1  ALONG  INTERSECTION  WITH  SATELLITE 
ORBIT  PLANE 

OXyZ  BASED  ON  ORBITAL  PLANE  WITH  OX  ALONG  OX 
OX'y'Z1  BASED  ON  ORBITAL  PLANE  WITH  Ox1  ALONG  OX1 

FIG. 3  AXIS  SYSTEMS  FOR  THEORY  OF  PERTURBATIONS 
DUE  TO  A  STATIONARY  DISTURBING  BODY 
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FIG. 4  AXIS  SYSTEMS  ASSOCIATED  WITH  A  DISTURBING  BODY 

IN  A  KEPLER  ORBIT 
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FIG.5  THE  OxYt"  AND  OX,,Y,,Z,,  AXIS  SYSTEMS 
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